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Abstract 



This thesis focuses on link scheduling in wireless mesh networks by taking into account 
physical layer characteristics. The assumption made throughout is that a packet is 
received successfully only if the Signal to Interference and Noise Ratio (SINR) at the 
receiver exceeds a certain threshold, termed as communication threshold. The thesis 
also discusses the complementary problem of flow control. 

First, we consider various problems on centralized link scheduling in Spatial Time 
Division Multiple Access (STDMA) wireless mesh networks. We motivate the use of 
spatial reuse as performance metric and provide an explicit characterization of spatial 
reuse. We propose link scheduling algorithms based on certain graph models (communi- 
cation graph, SINR graph) of the network. Our algorithms achieve higher spatial reuse 
than that of existing algorithms, with only a slight increase in computational complexity. 

Next, we investigate a related scenario involving link scheduling, namely random 
access algorithms in wireless networks. We assume that the receiver is capable of power- 
based capture and propose a splitting algorithm that varies transmission powers of users 
on the basis of quaternary channel feedback. We model the algorithm dynamics by a Dis- 
crete Time Markov Chain and consequently show that its maximum stable throughput 
is 0.5518. Our algorithm achieves higher maximum stable throughput and significantly 
lower delay than the First Come First Serve (FCFS) splitting algorithm with uniform 
transmission power. 

Finally, we consider the complementary problem of flow control in packet networks 
from an information-theoretic perspective. We derive the maximum entropy of a flow 
which conforms to traffic constraints imposed by a generalized token bucket regulator, by 
taking into account the covert information present in the randomness of packet lengths. 
Our results demonstrate that the optimal generalized token bucket regulator has a near 
uniform bucket depth sequence and a decreasing token increment sequence. 
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Chapter 1 



Introduction 

1.1 Link Scheduling in Wireless Networks 

Wireless and mobile communications have revolutionized the way we communicate over 
the past decade. This impact has been felt both in voice communications and wireless 
Internet access. The ever-increasing need for applications like video and images have 
driven the need for technologies like 3 rd Generation Partnership Project Long Term 
Evolution (3GPP LTE), 3 rd Generation Partnership Project 2 (3GPP2), IEEE 802.16 
Worldwide Interoperability for Microwave Access (WiMAX) networks and IEEE 802.11 
Wireless Local Area Networks (WLANs) which promise broadband data rates to wireless 
users. This revolution in wireless communications has had a great impact in India, where 
the number of cellular subscribers is 250 million (as of November 2008) and is growing 
at a rate of approximately 3% per month [2]. 

Wireless networks can be broadly classified into cellular networks and ad hoc net- 
works. A wireless ad hoc network is a collection of wireless nodes that can dynamically 
self-organize into an arbitrary topology to form a network without necessarily using any 
pre-existing infrastructure. Based on their application, ad hoc networks can be further 
classified into Mobile Ad Hoc Networks (MANETs), wireless mesh networks and wireless 
sensor networks. A wireless mesh network can be considered to be an infrastructure- 
based ad hoc network with a mesh backbone carrying most of the traffic. 

Wireless Mesh Networks (WMNs) have been recently advocated to provide connectiv- 
ity and coverage, especially in sparsely populated and rural areas. For example, several 
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Cellular network 

I I Mesh router Wireless link 

I I Mesh router with gateway/bridge Wired link 

Figure 1.1: Wireless mesh network, adapted from [1]. 

Wireless Community Networks (WCNs) are operational in Europe, Australia and USA 
[3]. Peer to peer wireless technology is also being developed by companies such as [4]. 
WMNs are dynamically self-organized and self-configured, with nodes in the network au- 
tomatically establishing an ad hoc network and maintaining mesh connectivity [1]. An 
example of a WMN is shown in Figure 1.1. Typically, a WMN comprises of two types 
of nodes: mesh routers and mesh clients. A mesh router consists of gateway/bridge 
functions and the capability to support mesh networking. Mesh routers have little or 
no mobility and form a wireless backbone for mesh clients. The gateway/bridge func- 
tionalities in mesh routers aid in the integration of WMNs with heterogeneous networks 
such as Ethernet [5], cellular networks, WLANs [6], WiMAX networks [7] and sensor 
networks. WMNs are witnessing commercialization in various applications like broad- 
band home networks, enterprise networks, community networks and metropolitan area 
networks. Moreover, WMNs diversify the functionalities of ad hoc networks, instead of 
just being another type of ad hoc network. These additional functionalities necessitate 
novel design principles and efficient algorithms for the realization of WMNs. 

Significant research efforts are required to realize the full potential of WMNs. Among 
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the many challenging issues in the design of WMNs, the design of the physical as well 
as the Medium Access Control (MAC) layers is important, especially from a perspective 
of achieving high network throughput. At the physical layer, techniques like adaptive 
modulation and coding, Orthogonal Frequency Division Multiplexing (OFDM) [8], [9] 
and Multiple Input Multiple Output (MIMO) techniques [10] can be used to increase 
the capacity of a wireless channel and achieve high data transmission rates. At the 
MAC layer, various solutions like directional antenna based MAC [11], MAC with power 
control [12] and multi-channel MAC [13] have been proposed in the literature. 

In this thesis, we primarily focus on the design of the MAC layer for wireless mesh 
networks. We abstract out essential features of the MAC and physical layers of a WMN 
and propose techniques that deliver high network throughput. We take into account 
wireless channel effects such as propagation path loss, fading and shadowing [14]. To- 
wards the end of the thesis, we provide an information-theoretic perspective on flow 
control. The main body of this thesis, however, focuses on MAC layer design for two 
types of networks: Spatial Time Division Multiple Access (STDMA) networks and ran- 
dom access networks. We next describe these two types of networks along with their 
potential applications in WMNs. 

An STDMA network can be thought of as a mesh network in which multiple trans- 
mitter receiver pairs can communicate at the same time. More specifically, consider 
a WMN comprising of store-and-forward nodes connected by "point to point" wireless 
communication channels (links). A link is an ordered pair (t, r), where t is a transmitter 
and r is a receiver. Time is divided into fixed-length intervals called slots. In STDMA, 
we allow concurrent communications between collections of nodes that are "reasonably 
far" from each other, i.e., we exploit spatial reuse. An STDMA link schedule describes 
the transmission rights for each time slot in such a way that communicating entities 
assigned to the same slot do not "collide" . In this thesis, we design centralized STDMA 
link scheduling algorithms that take into account physical layer characteristics such as 
Signal to Interference and Noise Ratio (SINR) at a receiver. 

STDMA link scheduling algorithms can be implemented at the MAC layer of wireless 
mesh networks, as shown in Figure 1.2. A mesh network can be constructed with mesh 
routers and mesh clients functioning as relay nodes in addition to their sender and 
receiver roles. The link schedule can be computed by a designated mesh router and then 
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Internet 




links scheduled by 
random access algorithm 



Figure 1.2: Potential applications of link scheduling and flow control in wireless networks. 

disseminated to all other nodes. The mesh routers form the mesh backbone to provide 
connectivity to (possibly mobile) mesh clients. 

In a related problem involving link scheduling, we consider a multipoint to point wire- 
less network with random access. When random access algorithms are directly translated 
from a wired network to a wireless network, they yield equal or lower throughput. This 
is because they do not consider the time variation of the wireless channel and interfer- 
ence conditions at the receiver. In this thesis, we design a distributed random access 
algorithm that takes into account wireless channel attributes such as propagation path 
loss and physical layer characteristics such as SINR at the receiver. 

Random access algorithms can be applied to the MAC layer of wireless networks, as 
shown in Figure 1.2. The BS and SSs are organized into a cell-like structure. Both uplink 
(from SS to BS) and downlink (from BS to SS) channels are shared among the SSs. This 
mode requires all SSs to be within the communication range and line of sight of the BS. 
A random access algorithm can be implemented in the SSs to resolve contentions on the 
uplink channel. 

In a complementary problem, we consider a packet level flow from a source to a 
destination over a data network. The packets transmitted by the source are regulated 
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at the ingress of the network, as shown in Figure 1.2. In this thesis, we investigate 
the maximum amount of information that can be transmitted from the source to the 
destination by utilizing the idea of covert information channels. 

To summarize, this thesis deals with the design of MAC layer algorithms (equiv- 
alently, link scheduling algorithms) for mesh networks. The proposed link scheduling 
algorithms take into account physical layer characteristics such as SINR at a receiver. 
Finally, we also consider the problem of flow control. 

Various solutions to the link scheduling problem have been proposed in literature 
depending on the modeling of the wireless network and interference conditions. In the 
next section, we motivate our work by briefly outlining the essential differences between 
our approach and the methodology of existing approaches. 

1.2 Motivation for the Thesis 

Consider the problem of determining a link schedule for an STDMA wireless network. 
STDMA link schedules can be classified into point to point and point to multipoint 
link schedules. In a point to point link schedule, the transmission right in each slot is 
assigned to certain links, while in a point to multipoint link schedule, the transmission 
right in each slot is assigned to certain nodes. An STDMA scheduling algorithm is a 
set of rules that is used to determine a link schedule so as to satisfy certain objectives. 
An STDMA link schedule should be so designed that, in every time slot, all packets 
transmitted by the scheduled transmitters are received successfully at the corresponding 
(intended) receivers. 

Two models have been proposed in literature for specifying the criteria for successful 
packet reception. According to the protocol interference model [15], a packet is received 
successfully at a receiver only if its intended transmitter is within the communication 
range and other unintended transmitters are outside the interference range of the re- 
ceiver. In essence, the protocol interference model mandates a "silence zone" around 
every scheduled receiver in a time slot. On the other hand, according to the physical 
interference model [15], a packet is received successfully at a receiver only if the SINR 
at the receiver is no less than a certain threshold, called communication threshold. 

Throughout this thesis, we assume that a packet is received successfully if the SINR at 
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the receiver is greater than or equal to the communication threshold, i.e., we employ the 
physical interference model. Moreover, we assume that, as long as the SINR threshold 
condition is satisfied at the receiver of a link, a constant rate of data transfer occurs along 
that link. In other words, the existence of a channel coding technique that guarantees a 
fixed data rate is assumed, when the SINR threshold condition is satisfied. 

To maximize the aggregate traffic transported by an STDMA wireless network, most 
link scheduling algorithms employ the protocol interference model and seek to minimize 
the schedule length. These algorithms model the network by a communication graph 
and employ novel techniques to color all the edges of the graph using minimum number 
of colors [16]. Such approaches have three lacunae. First, they transform the link 
scheduling problem to an edge coloring problem in a graph, which is a simplification 
of the true system model. Second, they do not incorporate wireless channel effects 
like propagation path loss, fading and shadowing. Finally, they do not consider SINR 
threshold conditions at a receiver. 

In this thesis, we seek to address these issues by designing polynomial time link 
scheduling algorithms that employ the physical interference model, provide a reasonably 
accurate representation of the wireless network and aim to maximize the number of 
successful packet transmissions per time slot. These algorithms take into account wire- 
less channel effects like propagation path loss, fading and shadowing, as well as SINR 
conditions at a receiver. We design and evaluate algorithms for both point to point and 
point to multipoint link scheduling. Our work falls under the realm of joint PHY-MAC 
design of wireless networks. 

In a related scenario involving link scheduling, consider the problem of designing a 
random access algorithm for a multipoint to point wireless network. When traditional 
random access algorithms like ALOHA [17] and tree-like algorithms [18] are employed 
in a wireless network, they yield equal or lower throughput compared to the wired case. 
This is because such algorithms are incognizant of wireless channel effects and physical 
layer characteristics. Thus, it is important to design a random access algorithm that 
incorporates wireless channel effects and exploits flexibilities provided by the physical 
layer. Towards this step, we assume a receiver that is capable of power-based capture 
[19]. Also, we assume that users can vary their transmission powers to increase the 
chances of successful packet reception under the physical interference model. Conse- 
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quently, we design and analyze a variable-power tree-like algorithm for a random access 
wireless network. 

In the final scenario, we formulate the problem of analyzing flow control in packet 
networks from an information-theoretic perspective. We focus on the problem of analyz- 
ing regulated flows in a point to point network. It is well-known that information (in the 
Shannon sense) can be transmitted from a source to a destination only by encoding it in 
the contents, lengths and timings of data packets from the source to the destination [20], 
[21]. We investigate the maximum amount of information that can be transmitted by a 
source whose flow is linearly bounded. Specifically, we assume that covert information 
is conveyed by randomness in packet lengths and investigate properties of the regulating 
mechanism that leads to maximum information transfer. 

1.3 Overview and Contributions of the Thesis 

In the first part of the thesis (Chapters 2 to 5), we consider various problems on central- 
ized link scheduling in STDMA wireless networks; each problem represents a different 
nuance of the overall link scheduling problem. In the second part of the thesis (Chapters 
6 and 7), we consider a related link scheduling problem, namely, distributed medium 
access control in a random access wireless network. In the third and final part of the 
thesis (Chapter 8), we consider flow control in networks from an information-theoretic 
perspective. 

Chapter 2 presents a generic framework and system model for link scheduling in 
STDMA wireless networks. We describe the system parameters of an STDMA wire- 
less network and explain two prevalent models used to specify the criteria for successful 
packet reception, namely protocol interference model and physical interference model 
[15]. We argue that STDMA link scheduling algorithms can be classified into three 
classes: algorithms based on modeling the network by a two-tier or communication 
graph, "hybrid" algorithms based on modeling the network by a communication graph 
and verifying SINR conditions and algorithms based on modeling the network by an 
SINR graph. We review representative research papers from each of these classes. We 
explain the relative merits and demerits of each class of algorithms in terms of com- 
putational complexity, performance and accuracy of the network model. We discuss 
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limitations of link scheduling algorithms based only on the communication graph model 
by providing illustrative examples. Finally to compare the performance of various link 
scheduling algorithms, we motivate and introduce spatial reuse as a performance metric. 
Various "spinoffs" of the "parent" link scheduling problem constitute the subproblems 
considered in Chapters 3, 4 and 5. 

In Chapter 3, we consider STDMA point to point link scheduling algorithms which 
utilize a communication graph representation of the network. Initially we examine the 
ArboricalLinkSchedule (ALS) algorithm [16], which represents the network by a com- 
munication graph, partitions the graph into minimum number of planar subgraphs and 
colors each subgraph in a greedy manner. We suggest a modification to the ALS al- 
gorithm based on reusing colors from previously colored subgraphs to color the current 
subgraph. We compare the performance of the modified algorithm with the ALS algo- 
rithm and derive its running time complexity. Subsequently, we propose the Conflict- 
FreeLinkSchedule algorithm, which is a hybrid algorithm based on the communication 
graph and verifying SINR conditions. Under various wireless channel conditions, we 
demonstrate that ConflictFreeLinkSchedule achieves higher spatial reuse than existing 
link scheduling algorithms based on the communication graph. However, this improve- 
ment in performance is achieved at a cost of slightly higher computational complexity. 

In Chapter 4, we consider the point to point link scheduling problem under the 
physical interference model. The STDMA network is represented by an SINR graph, in 
which weights of edges correspond to interferences between pairs of nodes and weights 
of vertices correspond to normalized noise powers at receiving nodes. We propose a 
link scheduling algorithm based on the SINR graph representation of the network. We 
prove the correctness of the algorithm and show that it has polynomial running time 
complexity. Finally, we demonstrate that the proposed algorithm achieves higher spatial 
reuse than ConflictFreeLinkSchedule. 

In Chapter 5, we consider point to multipoint link scheduling (broadcast scheduling) 
under the physical interference model. The problem addressed herein can be considered 
as the "dual" of the problem considered in Chapters 3 and 4. We generalize the definition 
of spatial reuse to the point to multipoint link scheduling problem. We propose a 
greedy scheduling algorithm which has demonstrably higher spatial reuse than existing 
algorithms, without any increase in computational complexity. 
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In Chapter 6, we consider another flavor of the link scheduling problem, namely 
random access algorithms for wireless networks. While random access algorithms for 
satellite networks, packet radio networks, multidrop telephone lines and multitap bus 
("traditional random access algorithms") is a well-researched and mature subject, the 
study of random access algorithms for wireless networks that take into account physical 
layer characteristics such as SINR and channel variations has yet to gain momentum. 
This chapter reviews representative research work which investigate such random access 
algorithms, most of them being generalizations of the ALOHA protocol (by adapting 
the retransmission probability) or the tree algorithm (by adapting the set of contending 
users). We motivate the use of variable transmission power to increase the throughput 
in random access wireless networks. 

We consider random access for wireless networks under the physical interference 
model in Chapter 7. We design an algorithm that adapts the set of contending users and 
their corresponding transmission powers based on quaternary (2 bit) channel feedback. 
We model the algorithm dynamics by a Discrete Time Markov Chain and subsequently 
derive its maximum stable throughput. Finally, we demonstrate that the proposed 
algorithm achieves higher throughput and substantially lower delay than the well-known 
First Come First Serve splitting algorithm [22]. 

In Chapter 8, we formulate the problem of analyzing flow control in packet networks 
from a perspective of maximizing mutual information between a source and a destination. 
We focus on the simpler, yet insightful, problem of analyzing regulated flows in a point 
to point network. More specifically, we consider a source whose flow is bounded by 
a "generalized" Token Bucket Regulator (TBR) and analyze the maximum amount of 
information (in the Shannon sense) that the source can convey to its destination by 
encoding information in the randomness of packet lengths. This chapter reveals two 
interesting results. First, under certain "bandwidth" constraints on cumulative tokens 
and cumulative bucket depth, we demonstrate that a generalized TBR can achieve higher 
flow entropy than that of a standard TBR. Second, we provide information-theoretic 
arguments for the observations that the optimal generalized TBR has a decreasing token 
increment sequence and a near-uniform bucket depth sequence. 

In Chapter 9, we summarize the thesis and provide possible directions for future 
work. Specifically, we suggest generalizations of the two-level power control algorithm 
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proposed in Chapter 7. We also provide pointers for deriving the approximation factors 
of the algorithms proposed in Chapters 3 and 4. 



Chapter 2 



A Framework for Link Scheduling 
Algorithms for STDMA Wireless 
Networks 

An STDMA wireless network consists of a finite set of nodes wherein multiple pairs of 
nodes can communicate concurrently, as discussed in Chapter 1. In this chapter, we 
outline a framework for modeling STDMA link scheduling algorithms. We consider a 
general representation of an STDMA wireless network, i.e., this model is not specific to 
any technology or protocol. This abstraction lends simplicity to the network model and 
helps us focus on the design of scheduling algorithms for the network. Since the problem 
of determining an optimal link schedule is NP-hard [16], researchers have proposed 
various heuristics to obtain close-to-optimal solutions. In our view, such heuristics can 
be broadly classified into three categories: algorithms based on modeling the network by 
a two-tier or communication graph, "hybrid" algorithms based on modeling the network 
by a communication graph and verifying SINR conditions and algorithms based on 
modeling the network by an SINR graph. We review representative research papers from 
each of these classes. The relative merits and demerits of each class of algorithms are 
also elucidated in the chapter. Our observations motivate us to propose a performance 
metric that is proportional to aggregate network throughput. 

The rest of this chapter is structured as follows. In Section 2.1, we describe the 
system model of an STDMA wireless network and explain the protocol and physical 



12 Chapter 2. A Framework for Link Scheduling Algorithms for STDMA Wireless Networks 



interference models. In Section 2.2, we elucidate the equivalence between a point to point 
link schedule for an STDMA network and the colors of edges of the communication graph 
model of the network. This is followed by a review of research work on point to point 
link scheduling algorithms based on the protocol interference model. In Section 2.3, we 
describe the limitations of algorithms based on the protocol interference model from a 
perspective of maximizing network throughput in wireless networks. We review research 
work on link scheduling algorithms based on the physical interference model in Sections 
2.4 and 2.5. Specifically, Section 2.4 reviews algorithms based on communication graph 
model of the network and SINR conditions, while Section 2.5 reviews algorithms based 
on an SINR graph model of the network. Finally, in Section 2.6, we propose spatial reuse 
as a performance metric and argue that it corresponds to network throughput from a 
physical layer viewpoint. 



2.1 System Model 

We consider a general model of an STDMA wireless network with N static store-and- 
forward nodes in a two-dimensional plane, where AMs a positive integer. Nodes are 
indexed as 1, 2, . . . , N. In a wireless network, a link is an ordered pair of nodes (t, r), 
where t is a transmitter and r is a receiver. We assume equal length packets. Time 
is divided into slots of equal duration. During a time slot, a node can either transmit, 
receive or remain idle. The slot duration equals the amount of time it takes to transmit 
one packet over the wireless channel. We make the following additional assumptions: 

• Synchronized nodes: All nodes are synchronized to slot boundaries. 

• Homogeneous nodes: Every node has identical receiver sensitivity, transmission 
power and thermal noise characteristics. 

• Backlogged nodes: We assume a node to be continuously backlogged, i.e., a node 
always has a packet to transmit and cannot transmit more than one packet in a 
time slot. 
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Let: 

(xj,yj) = Cartesian coordinates of node j =: r^, 

P = power with which a node transmits its packet, 

N = thermal noise power spectral density, 

D(j, k) = Euclidean distance between nodes j and k. 

The received signal power at a distance D from the transmitter is given by j^, where 
(3 is the path loss exponent 1 . An STDMA link schedule is a mapping from the set 
of links to time slots. We only consider static link schedules, i.e., link schedules that 
repeat periodically throughout the operation of the network. Let C denote the number 
of time slots in a link schedule, i.e., the schedule length. For a given time slot i, j th 
communicating transmitter- receiver pair is denoted by t it j — > r it j, where t it j denotes the 
index of the node which transmits a packet and r it j denotes the index of the node which 
receives the packet. Let Mj denote the number of concurrent transmitter-receiver pairs 
in time slot i. A point to point link schedule for the STDMA network is denoted by 
• • • ,S C ), where 

Si := {t it i — > r^i, • • • , t itMi — > n,Mi} 

= set of transmitter-receiver pairs which can communicate concurrently 
in time slot i. 

Note that a link schedule repeats periodically throughout the operation of the network. 
More specifically, transmitter-receiver pairs that communicate concurrently in time slot i 
also communicate concurrently in time slots i+C, i+2C and so on. Thus, Si = Si ( mo d c)- 
Finally, note that all transmitters and receivers are stationary. 
Every point to point link schedule must satisfy the following: 

• Operational constraint: During a time slot, a node can transmit to exactly one 
node, receive from exactly one node or remain idle, i.e., 

{Us, r i:j } D {t iik , r hk } = cf> V % = 1, . . . , C V 1 ^ j < k ^ M t . (2.1) 



We do not consider fading and shadowing effects. 



14 Chapter 2. A Framework for Link Scheduling Algorithms for STDMA Wireless Networks 



5 = (-75,0) 



l = (-40,5) 



Ay 



2= (0,0) 



3 = (95,0) 



X 

4= (135,0) 



f 6= (0,-75) 

(a) An STDMA wireless network with six nodes, 
point to point link schedule 



transmitter-receiver pairs 




2 -> 6\ 1 -> 5 l->5 
5->l\ 6^2 6^1 5^2 5->l 6^2 

il — ^ 2i3 — ^ 4i3 — > 2i4 — ^ 3i2 — > 5i2 — ^ lil — ^ 6.2 — ^ 3il — ^ 2i3 — ^ 4i3 — > 2i4 — ^ 3i ... 

1 TT~^ 8^ 2~~l 3^ 4 I ••• I time 

time slots 

(b) A point to point link schedule for the network shown in Figure 2.1(a). 
Figure 2.1: Example of STDMA network and point to point link schedule. 
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As an illustration, consider the STDMA wireless network shown in Figure 2.1(a). 
It consists of six nodes whose coordinates (in meters) are 1 = (—40,5), 2 = (0,0), 
3 = (95,0), 4 = (135,0), 5 = (-75,0) and 6 = (0,-75). An example point to point 
link schedule for this STDMA network is shown in Figure 2.1(b). Note that this sched- 
ule is only one of the several possible schedules and is given here only for illustrative 
purposes. The schedule length is C = 8 time slots and the schedule is defined by 
<S 2 , S 3 , S 4 , S 5 , S 6 , S 7 , S 8 ), where 

Si = {t M -> n,i} 

= {1-2}, 

^2 = {^2,1 — ?"2,1, ^2,2 — ► T 2,2, ^2,3 ~~ ¥ r 2,3} 

= {3-4, 5-1, 2-6}, 

$3 = {*3,1 — ^3,l} 

= {3-2}, 

<S 4 = {£4,1 — r 4,l) ^4,2 — r 4,2, ^4,3 — r 4,3} 

= {4^3, 6^2, 1-5}, 

<^5 = {^5,1 — ?"5,1, ^5,2—^ ^5,2} 

= {2^5, 6-1}, 

<^6 — {^6,1 — r 6,l} 

= {2-1}, 

$7 = {*7,1 — ^7,1, *7,2 — ^7,2} 
= {1-6, 5^2}, 

^8 — {^8,1 — r 8,l} 

= {2-3}. 

After 8 time slots, the schedule repeats periodically, as shown in Figure 2.1(b). 

A scheduling algorithm is a set of rules that is used to determine a link schedule 
^(•). Usually, a scheduling algorithm needs to satisfy certain objectives. 

Consider j th receiver in time slot i, i.e., receiver r^. The power received at r it j from 
its intended transmitter t it j (signal power) is D p, t P , - t - Similarly, the power received at 



Tij from its unintended transmitters (interference power) is Ylk=i WTtfk r- ■) • Thus, the 
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Signal to Interference and Noise Ratio (SINR) at receiver r itj is given by 

p 

SINR r . . = gpZi) _ (2.2) 

Without considering the interference power, the Signal to Noise Ratio (SNR) at receiver 
n,j is S iven h Y 

SNR r . . = Ar -. (2.3) 

According to the protocol interference model [15], transmission t it j — > r it j is successful 

if: 

1. the SNR at receiver r i: j is no less than a certain threshold j C: termed as the 
communication threshold. From (2.3), this translates to 

mi, ,.'■;./) ^ (j^y =■ R " ( 2 - 4 ) 

where R c is termed as communication range, and 

2. the signal from any unintended transmitter t^k is received at r^j with an SNR less 
than a certain threshold 7«, termed as the interference threshold. From (2.3), this 
translates to 

D(t itk ,r itj ) > (j^y =-Ri V fc = l,...,M i , k^j, (2.5) 

where Ri is termed as interference range. 

In essence, the transmission on a link is successful if the distance between the nodes is 
less than or equal to the communication range and no other node is transmitting within 
the interference range from the receiver. 

The STDMA network is denoted by $(JV, (n, . . . , r N ), P, j c , 7^ (3, N ). Note that 
< 7j < 7c, thus Ri > R c . The relation Ri = 2R C is widely assumed in literature [23], 
[24], [25], [26]. 

According to the physical interference model [15], the transmission on a link is suc- 
cessful if the SINR at the receiver is greater than or equal to the communication threshold 
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7 C . More specifically, the physical interference model states that transmission — > rjj 
is successful if: 

p 

v I ^ ; i — ^ 7c - (2 - 6) 

Note that the physical interference model is less restrictive but more complex. Usually, 
this representation has been employed to model mesh networks with TDMA like access 
mechanisms [27]. We will discuss this aspect later in the thesis. 

A point to point link schedule is conflict-free if the SINR at every intended 
receiver does not drop below the communication threshold, i.e., 

SINR r .. > 7c V i= 1,...,C, Vj = l,.,M, (2.7) 



2.2 Link Scheduling based on Protocol Interference 
Model 

2.2.1 Equivalence of Link Scheduling and Graph Edge Coloring 

In this section, we describe the communication and two-tier graph representations of an 
STDMA wireless network. We explain the equivalence between a point to point link 
schedule for the STDMA network and the colors of edges of the communication graph 
representation of the network, and illustrate this equivalence with an example. 

The STDMA network $(•) can be modeled by a directed graph C/(V, £), where V 
is the set of vertices and £ is the set of edges. Let V = {t>i, V2, . . . ,vn}, where vertex 
Vj represents node j in $(•). In the graph representation, if node k is within node j's 
communication range, then there is an edge from vj to Vk, denoted by vj A and termed 
as communication edge. Similarly, if node k is outside node j's communication range 
but within its interference range, then there is an edge from Vj to Vk, denoted by Vj A Vk 
and termed as interference edge. Thus, £ = £ c U £j, where £ c and £i denote the set of 
communication and interference edges respectively. The two-tier graph representation 
of the STDMA network $(•) is defined as the graph Q(V,£ C U £j) comprising of all 
vertices and both communication and interference edges. The communication graph 
representation of the STDMA network $(•) is defined as the graph Q C (V,£ C ) comprising 
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of all vertices and communication edges only. We will illustrate these representations 
with an example. 



Parameter 


Symbol 


Value 


transmission power 


P 


10 mW 


path loss exponent 


P 


4 


noise power spectral density 


N 


-90 dBm 


communication threshold 


7c 


20 dB 


interference threshold 


[i 


10 dB 



Table 2.1: System parameters for STDMA networks shown in Figures 2.1(a), 2.5 and 
2.10. 




Figure 2.2: Communication graph model of STDMA network described by Figure 2.1(a) 
and Table 2.1. 

Consider the STDMA wireless network <£>(•) whose deployment is shown in Figure 
2.1(a). The system parameters for this network are given in Table 2.1. From (2.4) and 
(2.5), it can be easily shown that R c = 100 m and Ri = 177.8 m. The corresponding 
communication graph representation Q C (V,S C ) is shown in Figure 2.2. The communi- 
cation graph comprises of 6 vertices and 14 directed communication edges. The vertex 
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Figure 2.3: Two-tier graph model of STDMA network described by Figure 2.1(a) and 
Table 2.1. 

and communication edge sets are given by 

V = {v 1 ,v 2 ,v 3 ,v 4 ,v 5 ,vq}, (2.8) 

c c c c c c c 

S c = {V!^ V 2 ,V 2 ^ Vi,V!^ V 5 ,V 5 ^ Vi,V!^ v e ,v 6 ^ Vi,V 2 ^ v 5 , 

v 5 A t> 2 , v 2 A v 6 , v 6 A v 2 , v 2 A v 3 , v 3 A v 2 , v 3 A v 4 , v 4 A V 3 } . (2.9) 

The two-tier graph model Q(V,£ C U £j) of the STDMA network $(•) is shown in Figure 
2.3. The two-tier graph comprises of 6 vertices, 14 directed communication edges and 
10 directed interference edges. The vertex and communication edge sets are given by 
(2.8) and (2.9) respectively, while the interference edge set is given by 

Si = {^i A v 4 ,v4 A Vl , v 2 A v 4 ,t> 4 A v 2 , v 3 A v 6 , v 6 A v 3 , 

v A A v 6 , v 6 A t> 4 , v 5 A t> 6 ,t; 6 A v 5 , }. (2.10) 

Given the above representations, a point to point link schedule \l/(-) for an STDMA 
wireless network $(•) can be considered as equivalent to assigning a unique color to 
every edge in the communication graph, such that transmitter-receiver pairs with the 
same color transmit simultaneously in a particular time slot. For the example network 
considered, the link schedule shown in Figure 2.1(b) corresponds to the coloring of the 
edges of the communication graph shown in Figure 2.4. Time slots 1, 2, 3, 4, 5, 6, 7 and 
8 in \l/(-) correspond to colors red, blue, green, magenta, yellow, cyan, brown and gold in 
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Figure 2.4: Edge coloring of communication graph shown in Figure 2.2 corresponding 
to the link schedule shown in Figure 2.1(b). 

S c respectively. Note that a coloring algorithm that uses the least number of colors also 
minimizes the schedule length. This aspect is further addressed in subsequent sections. 

2.2.2 Review of Algorithms 

In this section, we provide an overview of past research in the field of STDMA point to 
point link scheduling algorithms based on the protocol interference model. The protocol 
interference model is widely studied in literature because of its simplicity. It has been 
usually employed to model networks such as Carrier Sense Multiple Access with Collision 
Avoidance (CSMA/CA) based WLANs 2 [27], [25]. Centralized algorithms [16], [28], [29], 
[30], [25] as well as distributed algorithms [31] have been proposed for generating link 
schedules based on the protocol interference model. 

A link scheduling algorithm based on the protocol interference model utilizes a com- 
munication or two-tier graph model of the STDMA network to determine a point to 
point link schedule [32], [33]. Algorithms based on the protocol interference model for 
assigning links to time slots (equivalently, colors) require that two communication edges 

2 Consider an IEEE 802.11 based WLAN wherein CSMA with RTS/CTS/ACK is used to protect 
unicast transmissions. Due to carrier sensing, a transmission between nodes j and k may block all 
transmissions that are within a distance of Ri from either j (due to sensing RTS and DATA) or k (due 
to sensing CTS and ACK). 
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Vi — > Vj and — > vi can be colored the same if and only if: 

i. vertices v iy Vj, Vk, v\ are all mutually distinct, i.e., there is no primary edge conflict, 
and 

ii. Vi — > vi Q(-) and Vk — > Uj ^ <?(•)> i- e > there is no secondary edge conflict. 

The first criterion is based on the operational constraint (2.1). The second criterion 
states that a node cannot receive a packet if it lies within the interference range of 
any other transmitting node. A scheduling algorithm utilizes various graph coloring 
methodologies to obtain a non-conflicting link schedule, i.e., a link schedule devoid of 
primary and secondary edge conflicts. 

To maximize the throughput of an STDMA network, algorithms based on the pro- 
tocol interference model 3 seek to minimize the total number of colors used to color all 
the communication edges of (?(•). This will in turn minimize the schedule length. It is 
well known that for an arbitrary communication graph, the problem of determining a 
minimum length schedule (optimal schedule) is NP-hard [16], [29]. Hence, the approach 
followed in the literature is to devise algorithms that produce close to optimal (sub- 
optimal) solutions. The efficiency of a sub-optimal algorithm is typically measured in 
terms of its computational (run time) complexity and performance guarantee (approxi- 
mation factor). 

The concept of STDMA for wireless networks was formalized in [28]. The authors 
assume a multihop packet radio network with fixed node locations and consider the 
problem of assigning an integral number of slots to every link in an STDMA cycle 
(frame). To solve this problem, they model the network by a communication graph, 
determine a set of maximal cliques and then assign a certain number of slots to all the 
links in each maximal clique. Finally, the authors develop a fluid approximation for the 
mean system delay and validate it using simulations. 

In [29] , the authors consider pre-specified link demands in a spread spectrum packet 

radio network. They formulate the problem as a linear optimization problem and use 

3 Link scheduling algorithms based on the protocol interference model are sometimes referred to 
as "graph based algorithms" in literature [32], [33]. This term is slightly confusing since scheduling 
algorithms based on the physical interference model also construct graphs prior to determining a link 
schedule. 
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the ellipsoid algorithm [34] to solve the problem. They assume that the desired link data 
rates are rational numbers and develop a strongly polynomial algorithm 4 that computes 
a minimum length schedule. Finally, they consider the problem of link scheduling to 
satisfy pre-specified end-to-end demands in the network. They formulate this problem as 
a multicommodity flow problem and describe a polynomial time algorithm that computes 
a minimum length schedule. As pointed out by the authors, their algorithm is not 
practical due to its high computational complexity. 

A significant work in link scheduling under protocol interference model is reported in 
[16], in which the authors show that tree networks can be scheduled optimally, oriented 
graphs 5 can be scheduled near-optimally and arbitrary networks can be scheduled such 
that the schedule is bounded by a length proportional to the graph thickness 6 times the 
optimum number of colors. 

In [16], the authors seem to have missed a subtle point that colors from previously 
colored oriented graphs can be used to color the current oriented graph. Instead, they 
use a fresh set of colors to color each successive oriented graph. Consequently, their 
algorithm leads to a higher numbers of colors, especially if the number of oriented graphs 
is large. The authors employ such a heuristic primarily to upper bound the number of 
colors used by the algorithm ([16], Lemma 3.4) and consequently obtain bounds on the 
running time complexity and performance guarantee of the algorithm ([16], Theorem 
3.3). Though the ArboricalLinkSchedule algorithm has nice theoretical properties such 
as low computational complexity, it can be shown that it may yield a higher number of 
colors in practice. This leads to lower network throughput. 

We should point out here that, if we modify the ArboricalLinkSchedule algorithm to 
reuse colors from previously colored oriented graphs to color the current oriented graph, 
then the schedule length will always be lower than the schedule length obtained by the 

4 An algorithm is strongly polynomial if (a) the number of arithmetic operations (addition, multi- 
plication, division or comparison) is polynomially bounded by the dimension of the input, and (b) the 
precision of numbers appearing in the algorithm is bounded by a polynomial in the dimension and 
precision of the input. 

5 An in-oriented graph is a directed graph in which every vertex has at most one outgoing edge. An 

out-oriented graph is a directed graph in which every vertex has at most one incoming edge. 

6 The thickness of a graph (/(•) is the minimum number of planar graphs into which Q(-) can be 

partitioned. 
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ArboricalLinkSchedule algorithm. This can lead to higher network throughput. We 
develop this idea further in Chapter 3. Furthermore, we show that this can be achieved 
with only a slight increase in computational complexity. 

In [26], the authors investigate throughput bounds for a given wireless network and 
traffic workload under the protocol interference model. They use a conflict graph 7 to 
represent interference constraints. The problem of finding maximum throughput for a 
given source-destination pair under the flexibility of multipath routing is formulated 
as a linear program with flow constraints and conflict graph constraints. They show 
that this problem is NP-hard and describe techniques to compute lower and upper 
bounds on throughput. Finally, the authors numerically evaluate throughput bounds and 
computation time of their heuristics for simple network scenarios and IEEE 802.11 MAC 
(bidirectional MAC). Though the authors provide a general framework for joint routing 
and scheduling, they neither derive the computational complexity of their heuristics nor 
describe their link scheduling algorithm explicitly. 

Recently, in [25], the authors investigate joint link scheduling and routing under the 
protocol interference model for a wireless mesh network consisting of static mesh routers 
and mobile client devices. Assuming that l(u) denotes the aggregate traffic demand on 
node u, they consider the problem of maximizing A, such that at least Xl(u) amount 
of traffic can be routed from each node u to a fixed gateway node. Since this problem 
is NP-hard, the authors propose heuristics based on linear programming and re-routing 
flows on the communication graph. They derive the worst case bound of their algorithm 
and evaluate its performance via simulations. Though the authors make a reasonable 
attempt to solve the joint routing and scheduling problem, their algorithm is extremely 
complex 8 and brute force in nature. Furthermore, the authors have not provided intuitive 
arguments for their algorithm. 

Another recent work which jointly investigates link scheduling and routing under 

7 Under the protocol interference model, the conflict graph F(Vf, Ef) is constructed from the com- 
munication graph Q C (V, £ c ) as follows. Let Uj denote the communication edge Vi —> Vj. Vertices of F(-) 
correspond to directed edges Uj in £ c . In F(-), there exists an edge from vertex lij to vertex l pq if any 

of the following is true: (a) D(i, q) < i?j or (b) D(p,j) < iij. 

8 The algorithm in [25] consists of five steps: solve linear program, channel assignment, post process- 
ing, flow scaling and interference free link scheduling. Moreover, the channel assignment step consists 
of three algorithms. 
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protocol interference model is reported in [30]. The authors consider wireless mesh 
networks with half duplex and full duplex orthogonal channels, wherein each node can 
transmit to at most one node and/or receive from at most k nodes (k ^ 1) during 
any time slot. They investigate the joint problem of routing and scheduling to analyze 
the achievability of a given rate vector between multiple source-destination pairs. The 
scheduling algorithm is equivalent to an edge-coloring on a multi-graph representation 9 
and the corresponding necessary conditions lead the routing problem to be formulated as 
a linear optimization problem. The authors describe a polynomial time approximation 
algorithm to obtain an e-optimal solution of the routing problem using the primal dual 
approach. Finally, they evaluate the performance of their algorithms via simulations. 

It has been observed that high data rates are achievable in a wireless mesh network 
by allowing a node to transmit to only one neighboring node at fixed peak power in 
any time slot [30] . We point out here that a similar assumption of uniform transmission 
power has been made in our system model in subsequent chapters of the thesis. 

Algorithms based on the protocol interference model represent the network by a 
communication or two-tier graph and employ a plethora of techniques from graph theory 
[35] and approximation algorithms [36], [37] to devise heuristics which yield a minimum 
length schedule. Consequently, such algorithms have the advantage of low computational 
complexity (in general). However, recent research suggests that these algorithms result 
in low network throughput. This aspect is further illustrated in the following section. 

2.3 Limitations of Algorithms based on Protocol In- 
terference Model 

Due to its inherent simplicity, the protocol interference model has been traditionally 
employed to represent a wide variety of wireless networks. However, it leads to low 
network throughput in wireless mesh networks. To emphasize this point, we provide 
examples to demonstrate that algorithms based on the protocol interference model can 
result in schedules that yield low network throughput. 

9 A multi-graph is a directed graph in which multiple edges can emanate from a vertex u, and 
terminate at another vertex Vj (vj ^ Vi). 
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Intuitively, the protocol interference model divides the deployment region of the 
STDMA wireless network into "communication zones" and "interference zones". This 
transforms the scheduling problem to an edge coloring problem for the communication 
graph representation of the network. However, this simplification can result in schedules 
that do not satisfy the SINR threshold condition (2.7). 

Specifically, algorithms based on the protocol interference model do not necessarily 
maximize the throughput of an STDMA wireless network because: 



1. They can lead to high cumulative interference at a receiver, due to hard-thresholding 
based on communication and interference radii [32], [33]. This is because the SINR 
at receiver r^j decreases with an increase in the number of concurrent transmis- 
sions Mj, while the communication radius R c and the interference radius Ri have 
been defined for a single transmission only. 

Y 



1 = (-360,0) 3 =(90,0) 6 = (450,0) x 

2 =(-450,0) [4 =(0,0) 5 = (360,0) 



Figure 2.5: An STDMA wireless network with six nodes. 






Figure 2.6: Two-tier graph model of the STDMA wireless network described by Figure 
2.5 and Table 2.1. 
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Figure 2.7: Subgraph of two-tier graph shown in Figure 2.6. 



26 Chapter 2. A Framework for Link Scheduling Algorithms for STDMA Wireless Networks 



V2 



Vi 




(§)J5) 




V5 



v 6. 



Figure 2.8: Coloring of subgraph shown in Figure 2.7. 
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Figure 2.9: Point to point link scheduling algorithms based on protocol interference 
model can lead to high interference. 



For example, consider the STDMA wireless network whose deployment is shown 
in Figure 2.5. The network consists of six labeled nodes whose coordinates (in 
meters) are 1 = (-360,0), 2 = (-450,0), 3 = (90,0), 4 = (0,0), 5 = (360,0) and 
6 = (450, 0). The system parameters are shown in Table 2.1, which yield R c = 100 
m and Ri = 177.8 m. The two-tier graph model of the STDMA network is shown 
in Figure 2.6; note that interference edges are absent. Consider the transmission 
requests 1 — > 2, 3 — > 4 and 5^6, which correspond to communication edges of 
the subgraph shown in Figure 2.7. The communication edges V\ A v 2 , v s — > t>4 
and A vq shown in Figure 2.7 do not have primary or secondary edge conflicts. 
To minimize the number of colors, such an algorithm will color these edges with 
the same color, as shown in Figure 2.8. Equivalently, transmissions 1 — >• 2, 3 — > 4 
and 5 — > 6 will be scheduled in the same time slot, say time slot i. However, 
our computations show that the SINRs at receivers r^i, and are 21.26 dB, 
18.42 dB and 19.74 dB respectively. Figure 2.9 shows the nodes of the network 
along with the labeled transmitter-receiver pairs, receiver-centric communication 
and interference zones and the SINRs at the receivers. From the SINR threshold 
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condition (2.6), transmission t^i — > r^i is successful, while transmissions t i)2 — > r i)2 
and t it3 — "> r i 3 are unsuccessful. This leads to low network throughput. 

2. Moreover, these algorithms can be extremely conservative and result in higher 
number of colors. 
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Figure 2.10: An STDMA wireless network with four nodes. 





Figure 2.11: Two-tier graph model of STDMA wireless network described by Figure 2.10 
and Table 2.1. 



For example, consider the STDMA wireless network whose deployment is shown 
in Figure 2.10. The network consists of four labeled nodes whose coordinates (in 
meters) are 1 = (0,0), 2 = (50,0), 3 = (220,0) and 4 = (170,0). The system 
parameters are shown in Table 2.1, which lead to R c = 100 m and Ri = 177.8 m. 
The two-tier graph model of the STDMA network is shown in Figure 2.11. Consider 
the transmission requests 1 — > 2 and 3 — > 4, which correspond to communication 
edges of the subgraph shown in Figure 2.12. The communication edges vi A v 2 
and i>3 A t> 4 shown in Figure 2.12 have secondary edge conflicts. Hence, such an 
algorithm will typically color these edges with different colors, as shown in Figure 
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Figure 2.12: Subgraph of two-tier graph shown in Figure 2.11. 




Figure 2.13: Coloring of subgraph shown in Figure 2.12. 

2.13. Equivalently, a link scheduling algorithm based on the protocol interference 
model will schedule transmissions 1 — ► 2 and 3 — > 4 in different time slots, say 
time slots i and j respectively, where i ^ j. Our computations show that the 
resulting SINRs at receivers r^i and r^i are both equal to 32.04 dB. Figure 2.14 
shows the nodes of the network along with the labeled transmitter-receiver pairs, 
receiver-centric communication and interference zones and SINRs at the receivers. 
Observe that, with an algorithm based on the protocol interference model, the 
SINRs at both receivers are well above the communication threshold of 20 dB. 
Alternatively, consider an algorithm (perhaps based on the physical interference 
model) that schedules transmissions 1 — > 2 and 3 — ► 4 in the same time slot, 
say time slot i. The corresponding edge coloring is shown in Figure 2.15. Our 
computations show that the resulting SINRs at receivers r^i and are both 
equal to 20.91 dB, which are also above the communication threshold. Figure 
2.16 shows the nodes of the network along with the labeled transmitter-receiver 
pairs and SINRs at the receivers. In essence, with the alternate algorithm, both 
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SINR,, 1 = 32.04 dB SINR^ = 32.04 dB 

Figure 2.14: Point to point link scheduling algorithms based on protocol interference 
model can lead to higher number of colors. 




Figure 2.15: Alternative coloring of subgraph shown in Figure 2.12. 
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Figure 2.16: A point to point link schedule corresponding to Figure 2.15 that yields 
lower number of colors. 

transmissions t^i — > r^i and £^2 — > r i 2 are successful, since signals levels are so 
high at the receivers that strong interferences can be tolerated. In summary, a 
point to point link scheduling algorithm based on the protocol interference model 
will typically schedule the above transmissions in different slots and yield lower 
network throughput compared to the alternate algorithm. 

3. Lastly, these algorithms are not aware of the topology of the network, i.e., they 
determine a link schedule without being cognizant of the exact positions of the 
transmitters and receivers. 

The above examples demonstrate that scheduling algorithms based on the protocol 
interference model can result in low network throughput. Observe that algorithms that 
construct an approximate model of the STDMA network (two tier graph or communi- 
cation graph) and focus on minimizing the schedule length do not necessarily maximize 
network throughput. This observation is developed into a proposal for an appropriate 
performance metric in Section 2.6. 

Since link scheduling algorithms based on the protocol interference model yield low 
throughput, researchers have propounded algorithms based on the physical interference 
model to improve the throughput of STDMA wireless networks. To achieve higher 
throughput, one possible technique is to model the STDMA network by a communication 
graph and check SINR threshold conditions during assignment of links to time slots; 
this is the approach most commonly employed, for example in [27], [32], [38]. The other 
technique is to incorporate SINR threshold conditions into a special graph model of 
the network; this approach is more challenging and (to the best of our knowledge) is 
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considered only in research work such as [39] , [40] , [41] . Research papers which employ 
the former approach are reviewed in Section 2.4, while research papers which employ 
the latter approach are reviewed in Section 2.5. 

2.4 Link Scheduling based on Communication Graph 
Model and SINR Conditions 

In this section, we examine recent research in link scheduling based on modeling the 
STDMA network by a communication graph and verifying SINR conditions at the re- 
ceivers. Though algorithms based on this model [24], [42], yield higher throughput, they 
usually result in higher computational complexity than algorithms based on the protocol 
interference model. 

In [27], the authors investigate throughput improvement in an IEEE 802.11 like 
wireless mesh network with CSMA/CA channel access scheme replaced by STDMA. 
For a successful packet transmission, they mandate that two-way communication be 
successful, i.e., a packet transmission is defined to be successful if and only if both data 
and acknowledgement packets are received successfully. Under this "extended physical 
interference model", they present a greedy algorithm which computes a point to point 
link transmission schedule in a centralized manner. Assuming uniform random node 
distribution and using results from occupancy theory [43] , they derive an approximation 
factor for the length of this schedule relative to the shortest schedule. Though the 
analysis presented in [27] is novel, their model is restrictive because it is only applicable 
to wireless networks using link-layer reliability protocols. 

The throughput performance of link scheduling algorithms based on two-tier graph 
model G(V, £ c U £*) has been analyzed under physical interference conditions in [32]. The 
authors determine the optimal number of simultaneous transmissions by maximizing 
a lower bound on the throughput and subsequently propose Truncated Graph-Based 
Scheduling Algorithm (TGSA), an algorithm that provides probabilistic guarantees for 
network throughput. Though the analysis presented in [32] is mathematically elegant 
and based on the Edmundson-Madansky bound [44] , [45] , their algorithm does not yield 
high network throughput. This is because the partitioning of a maximal independent 
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set of communication edges into multiple subsets (time slots) is arbitrary and not based 
on network topology, which can lead to significant interference in certain regions of the 
network. This is further elucidated by the simulation results in Chapter 3. 

The performance of algorithms based on the protocol interference model versus those 
based on communication graph model and SINR conditions is evaluated and compared 
in [33]. To generate a non-conflicting link schedule based on the protocol interference 
model, the authors use a two-tier graph model with certain SINR threshold values chosen 
based on heuristics and examples. To generate a conflict-free point to point link sched- 
ule based on the physical interference model, the authors employ a method suggested 
in [46] which describes heuristics based on two path loss models, namely terrain-data 
based ground wave propagation model and Vogler's five knife-edge model. Their sim- 
ulations results indicate that, under a Poisson arrival process, algorithms based on the 
protocol interference model result in higher average packet delay than algorithms based 
on communication graph model and SINR conditions. 

In [42] , the authors investigate the tradeoff between the average number of concur- 
rent transmissions (spatial reuse) and sustained data rate per node for an IEEE 802.11 
wireless network. They show that spatial reuse depends only on the ratio of transmit 
power to carrier sense threshold [6] . Keeping the carrier sense threshold fixed, they pro- 
pose a distributed power and rate control algorithm based on interference measurement 
and evaluate its performance via simulations. 

In [24], the authors investigate mitigation of inter-flow interference in an IEEE 
802. lie wireless mesh network from a temporal-spatial diversity perspective. Measure- 
ments of received signal strengths are used to construct a virtual coordinate system 
to identify concurrent transmissions with minimum inter-flow interference. Based on 
this new coordinate system, one of the nodes, designated as gateway node, determines 
the scheduling order for downlink frames of different connections. Through extensive 
simulations with real-life measurement traces, the authors demonstrate throughput im- 
provement with their algorithm. 

Algorithms based on representing the network by a communication graph and veri- 
fying SINR threshold conditions yield higher network throughput than algorithms based 
on the protocol interference model. However, this is achieved at the cost of higher com- 
putational complexity. Furthermore, the gains in throughput may not be significant 
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enough to justify the increase in computational complexity. This has prompted few 
researchers to solve the link scheduling problem in a more fundamental manner. These 
researchers have proposed an altogether different model of the network, termed as SINR 
graph model, and developed heuristics. Such algorithms are reviewed in the following 
section. 

2.5 Link Scheduling based on SINR Graph Model 

In literature, many authors refer to algorithms based on communication graph model and 
checking SINR conditions as "algorithms based on physical interference model" . In this 
thesis, only algorithms that embed SINR threshold conditions into an appropriate graph 
model of the network are referred to as "algorithms based on the physical interference 
model". Though the physical interference model is more realistic, algorithms based 
on this model [39], [40], [41] have, in general, higher computational complexity than 
algorithms based on the protocol interference model. 

Point to point link scheduling for power-controlled STDMA networks under the phys- 
ical interference model is analyzed in [39]. The authors define scheduling complexity as 
the minimum number of time slots required for strong connectivity of the graph 10 con- 
structed from the point to point link schedule. They develop an algorithm employing 
non-linear power assignment 11 and show that its scheduling complexity is polylogarith- 
mic in the number of nodes. In a related work [40], the authors investigate the time 
complexity of scheduling a set of communication requests in an arbitrary network. They 
consider a "generalized physical model" wherein the actual received power of a signal 
can deviate from the theoretically received power by a multiplicative factor. Their algo- 
rithm successfully schedules all links in time proportional to the squared logarithm of the 
number of nodes times the static interference measure [47]. Though the authors of [39], 
[40] allow non-uniform transmission power at all nodes and develop novel algorithms, 

10 A directed graph Q(-) is strongly connected if there exists a directed path from every vertex to 
every other vertex. 

n In uniform power assignment, all nodes transmit with the same transmission power. In linear power 
assignment [39] , a node transmits with minimum power required to satisfy the SINR threshold condition 
at the receiver, i.e., transmission power equals Noj c D^ . Non- linear power assignment refers to a power 
assignment scheme that is neither uniform nor linear. 
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their algorithms are impractical. This is because wireless devices have constraints on 
maximum transmission power, while the algorithms in [39], [40] can result in arbitrarily 
high transmission power at some nodes. 

In [26], the authors provide a general framework for computation of throughput 
bounds for a given wireless network and traffic workload. Though their work primarily 
focuses on the protocol interference model, they briefly allude to the physical interference 
model too. Specifically, they describe a technique to construct a weighted conflict graph 
to represent interference constraints. They briefly describe methods to compute lower 
and upper bounds on throughput and the issues involved therein. However, the authors 
do not describe simulation results under the physical interference model, perhaps due 
to the tremendous complexity incurred in solving linear programs for representative 
network scenarios. 

Remark 2.5.1. Under physical interference model, the weighted conflict graph F(Vp, E F ) 
[26] is constructed from the network as follows. Let Sij := jy£ ^ denote the received 
signal power at node j due to the transmission from node i. In F(-), a vertex corresponds 

S' 

to a directed link kj (equivalently, node pair (i,j)) provided ^ j c . F(-) is a perfect 
graph wherein the weight wf- of the directed edge from vertex l vq to vertex Uj is given by 



We should point out here that, analogous to a conflict graph, an SINR graph rep- 
resentation of an STDMA wireless network has been proposed by us in Chapter 4. 
Furthermore, the authors of [26] do not propose any specific link scheduling algorithm 
and use the weighted conflict graph only to compute bounds on network throughput. On 
the other hand, we use an SINR graph representation of the network under the physical 
interference model and develop a link scheduling algorithm with lower time complexity 
and demonstrably superior performance. 

More specifically, in Chapter 4, we investigate link scheduling for STDMA wireless 
networks under the physical interference model. Unlike [39], [40], we assume that a 
node transmits at fixed power, i.e., we assume uniform power assignment. Moreover, 
unlike [39], [40], we do not assume a minimum distance of unity between any two nodes. 
Consequently, our system model is more practical than those of [39], [40]. Under these 
realistic assumptions, we propose a link scheduling algorithm based on an SINR graph 
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representation of the network. In the SINR graph 12 , weights of the edges correspond 
to interferences between pairs of nodes. We prove the correctness of the algorithm 
and derive its computational complexity. We demonstrate that the proposed algorithm 
achieves higher throughput than existing algorithms, without any increase in computa- 
tional complexity. 

So far, we have provided a brief glimpse into three classes of link scheduling algo- 
rithms, each with its relative merits and demerits. For example, algorithms based on 
the protocol interference model have low computational complexity and are simple to 
implement, but yield low network throughput. On the other hand, algorithms based on 
SINR graph representation have higher computational complexity and are more cumber- 
some to implement, but achieve higher network throughput. Also, there exist algorithms 
based on communication graph and SINR conditions whose performance characteristics 
lie between these two classes. Hence, in general, these three classes of algorithms exhibit 
a tradeoff between complexity and performance. Finally, algorithms based on the pro- 
tocol interference model are better suited to model WLANs, while the latter two classes 
of algorithms are better suited to model wireless mesh networks. For these reasons, we 
investigate and develop algorithms from each of these classes in this thesis. 

Prior to proposing efficient algorithms in each of these classes, we seek to address 
the following question: Is schedule length an appropriate performance metric for an al- 
gorithm that considers the SINR threshold condition (2.6) as the criterion for successful 
packet reception? In other words, should algorithms based on communication graph and 
SINR conditions and algorithms based on SINR graph representation focus on minimiz- 
ing the schedule length? We answer this important question in detail in the following 
section. 



2.6 Spatial Reuse as Performance Metric 

In literature, link scheduling algorithms have only focused on minimizing the sched- 
ule length. However, algorithms that minimize the schedule length do not necessarily 
maximize network throughput, as explained in Section 2.3. Thus, from a perspective 

12 Thc SINR graph is analogous to a line graph [35] constructed from the communication graph 
representation of the network. 
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of maximizing network throughput observed by the physical layer, it is imperative to 
consider a performance metric that takes into account SINR threshold condition (2.6) as 
the criterion for successful packet reception, i.e., a metric also suitable for the physical 
interference model. We propose such a performance metric, spatial reuse, in this sec- 
tion. We show that maximizing spatial reuse directly translates to maximizing network 
throughput. 

Consider an STDMA wireless network that operates over (k 2 — ki + 1) time slots 
ki, ki + 1, . . . , k 2 — l,k 2 . The total number of successfully scheduled links from slot ki 
to slot k 2 is 

k 2 Mi 

T[h,k 2 ] = ^^/(SINR rij > 7c ). (2.11) 

i=k\ j=l 

So, the number of successfully scheduled links per time slot from slot k x to slot k 2 is 

= fc2 _ fcl + 1 • ( 2 - 12 ) 

We define spatial reuse a as the limiting value of rj[ki, k 2 ] (assuming that the limit 
exists). In other words, spatial reuse is the limiting value of r}[ki, k 2 ] as the duration of 
the time interval becomes very large. Mathematically, 

Spatial Reuse := lim f][ki,k 2 ], 

;.a = lim ^ =fclZ ^r 13 1 ' . (2.13) 

|fe 2 -fei|-»oo k 2 — ki + 1 

Assuming a constant data rate of R bits per second on each successful link and a 
slot duration of t s seconds, the (aggregate) network throughput is given by <tRt s bits 
per second. Thus, spatial reuse is directly proportional to network throughput. Note 
that spatial reuse is cognizant of the physical interference model, thereby making it an 
appropriate performance metric for the comparison of various link scheduling algorithms. 

The fact that the interference at a receiver is an increasing function of the number 
of concurrent transmissions in a time slot limits the value of spatial reuse (for a given 
STDMA network). More specifically, if too many transmissions are scheduled in a single 
time slot, the interference at some receivers will be high enough to drive the SINRs 
below the communication threshold, leading to lower spatial reuse. Therefore, for a 
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given STDMA network, there are certain fundamental limits (upper bounds) on the 
spatial reuse. 

In our system model, we only consider static link schedules, i.e., the same fixed 
pattern of slots repeats cyclically. Hence, for our system model, the equation for spatial 
reuse, (2.13), can be simplified to 



The essence of STDMA is to have a reasonably large number of concurrent and 
successful transmissions. For a network which is operational for a long period of time, 
say L time slots, the total number of successfully received packets is La. Thus, a high 
value of spatial reuse directly translates to higher network throughput and the number 
of colors C is relatively unimportant. Hence, spatial reuse 13 turns out to be a crucial 
metric for the comparison of different STDMA algorithms in Chapters 3, 4 and 5. 



Spatial Reuse — a — 



Eti EfJi /(SINR rij > 7c ) 



(2.14) 



c 



13 Note that spatial reuse in our network model is analogous to spectral efficiency in digital commu- 
nication systems. Both performance metrics correspond to the "rate of data transfer" and are upper 
bounded by their respective system parameters. 



Chapter 3 



Point to Point Link Scheduling 
based on Communication Graph 
Model 



We begin our investigation in link scheduling by critically examining the ArboricalL- 
inkSchedule algorithm proposed in [16]. The algorithm is based only on the communi- 
cation graph (protocol interference model) and seeks to minimize the schedule length. 
Though ArboricalLinkSchedule has good properties such as low computational com- 
plexity, it can yield higher schedule length in practice. Towards this end, we propose a 
novel modification to ArboricalLinkSchedule that results in lower schedule length. We 
compare the performance of the modified algorithm with the ArboricalLinkSchedule al- 
gorithm and derive its run time (computational) complexity in Section 3.1. We then 
propose the ConflictFreeLinkSchedule point to point link scheduling algorithm, which is 
based on communication graph model and SINR conditions, in Section 3.2. The perfor- 
mance of the proposed algorithm is compared with existing link scheduling algorithms 
under various wireless channel conditions. We show that the proposed algorithm has 
polynomial run time complexity. Finally, we summarize the implications of our work. 
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3.1 ArboricalLinkSchedule Algorithm Revisited 

In this section, we propose a modification to the ArboricalLinkSchedule point to point 
link scheduling algorithm. Since both the original algorithm and the proposed modifi- 
cation are based on the protocol interference model, we compare their performance in 
terms of average schedule length. Finally, we also derive the run time complexity of the 
modified algorithm. 

Our system model and notation are same as described in Section 2.2. We seek an 
algorithm that determines a minimum length point to point link schedule for an STDMA 
wireless network under the protocol interference model. For consistency with the graph 
model described in [16], we assume that the STDMA wireless network $(•) is modeled 
by the communication graph <? C (V,£ C ) only, i.e., interference edges are absent (£ = £ c ). 

It is well known that, under the protocol interference model, the problem of de- 
termining an optimal schedule, i.e., a minimum length schedule, is NP-hard [48]. As 
pointed out in Section 2.2.1, this is closely related to the problem of coloring all edges 
of the communication graph with minimum number of colors, which is also known to 
be NP-hard [16] . Consequently, the only recourse is to devise approximation algorithms 
(heuristics) and show their efficiency theoretically and experimentally. 

One such algorithm, ArboricalLinkSchedule, has been described in [16]. First, the 
algorithm uses the labeler function to label all the vertices of the communication graph. 
Next, it partitions the communication graph into edge-disjoint subgraphs, which are 
termed as "oriented graphs" . Finally, the oriented graphs are colored in sequence. Specif- 
ically, the vertices in each oriented graph are scanned in increasing order of label and 
the unique edge associated with each vertex is colored using the NonConflictingEdge 
function [16]. The labeler function and the partitioning technique are described later in 
the section. 

In [16], the authors appear to have missed a delicate point that colors from previously 
colored oriented graphs can be used to color the present oriented graph. Specifically, 
they use a fresh set of colors to color each successive oriented graph. In our opinion, the 
authors employ this method to upper bound the number of colors used by the algorithm 
([16], Lemma 3.4) and thus derive the running time complexity of the algorithm ([16], 
Theorem 3.3). However, such a heuristic can potentially lead to a higher number of 
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colors (and higher schedule length) in practice. 

Therefore, we propound a modification to the ArboricalLinkSchedule algorithm that 
reuses colors from previously colored oriented graphs to colors the current oriented graph. 
The resulting schedule length will always be lower than that of ArboricalLinkSchedule, 
leading to potentially higher throughput. Our proposed link scheduling algorithm is 
ALSReuseColors, which considers the communication graph Q C (V, E) and is described 
in Algorithm 1. 

In Phase 1, we label all the vertices using the labeler function [16]. The labeler 
function is reproduced in Algorithm 2 for convenience. It is a recursive function that 
assigns a unique label (from 1 to N) to every vertex of the communication graph. Let 
L{w) denote the label assigned to vertex w. The notation Q r \{u} denotes the graph that 
results when vertex u and all its incident edges are removed from graph Q r (-). At every 
step in the recursion, it chooses the minimum degree vertex u in the residual graph Q r {-) 
and assigns it the highest label that has not been assigned so far. Note that vertices 
with lower degree tend to be assigned higher labels. The labeler function ensures that, 
for any given node, the number of neighbors with lower labels is much lower than the 
number of vertices in Q c {-). 

In Phase 2, the communication graph Q c (-) is decomposed into what are called as 
out-oriented and in-oriented graphs 71, T 2 , . . . , Tk, similar to the technique employed in 
[16]. Recall that an in-oriented graph is a directed graph in which every vertex has at 
most one outgoing edge, while an out-oriented graph is a directed graph in which every 
vertex has at most one incoming edge. Each Tj is a forest 1 and every edge of G c {-) is 
in exactly one of the T;'s. This decomposition is achieved by partitioning graph G c {-), 
the undirected equivalent of G c {')-> into undirected forests. The number of forests can be 
minimized by using techniques from Matroid theory ([49], A;-forest problem). However, 
this optimal decomposition requires extensive computation. Hence, we adopt a faster 
albeit non-optimal approach of using successive breadth first searches [50] to decompose 
G c (-) into undirected forests. Each undirected forest is further mapped to two directed 
forests. In one forest, the edges in every connected graph point away from the root and 
every vertex has at most one incoming edge, thus producing an out-oriented graph. In 
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the other forest, the edges in every connected graph point toward the root and every 
vertex has at most one outgoing edge, thus producing an in-oriented graph. 

In Phase 3, the oriented graphs are considered sequentially. For each oriented graph, 
the vertices are considered in increasing order of label and the unique edge associated 
with each vertex is colored using the NCEReuseColors function. The NCEReuseColors 
function is explained in Algorithm 3. For the edge x under consideration, it discards 
any color from any oriented graph that has an edge with a primary or secondary conflict 
with x. It returns the least color among the residual set of non-conflicting colors from 
all oriented graphs colored so far. If no non-conflicting color from any oriented graph is 
found, it returns a new color. 

Algorithm 1 ALSReuseColors 
1: input: Directed communication graph £? C (V, £) 

2: output: A coloring C : £ — > {1, 2, . . .} 

3: n <— labeler(^ c ) {Phase 1} 

4: use successive breadth first searches to partition Q c (-) into oriented graphs T h 1 ^ 

% ^ k {Phase 2} 
5: for i <— 1 to k do {Phase 3 begins} 
6: for j <— 1 to n do 
7: if Tj is out-oriented then 
8: let x = (s, d) be such that L(d) = j 

9: else 

10: let x = (s, d) be such that L{s) = j 

11: end if 

12: C(x) <- NCEReuseColors (x) 

13: end for 

14: end for{Phase 3 ends} 



3.1.1 Performance Results 

In the simulation experiment, every node location is generated randomly, using a uniform 
distribution for its X and Y coordinates in the deployment area. We assume that 
the deployment region is a square of length L. Thus, if (Xj,Yj) are the Cartesian 
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Algorithm 2 integer labeler(£ r ) 
1: if Q r { ) is not empty then 

2: let u be a vertex of G r {-) of minimum degree 

3: L(u) <- 1 + labeler(£ r \ {u}) 

4: else 

5: return 

6: end if 



Algorithm 3 integer NCEReuseColors(a;) 
1: input: Directed communication graph £? C (V, £) 

2: output: A non-conflicting color 

3: C <— set of existing colors 

4: C\ <— {C(h) : h is colored and x and h have a primary edge conflict} 

5: C 2 <— {C*(/i) : h is colored and x and /i have a secondary edge conflict} 

6: C nc = C\ {& U C 2 } 

7: if C nc 7^ then 

8: return the least color G C nc 

9: else 

10: return \C\ + 1 
11: end if 
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coordinates of j th node, then Xj ~ U[0,L] and Yj ~ U[0,L]. The values chosen for 
system parameters P, 7 C , (3 and N , are prototypical values of system parameters in 
wireless networks [42] . After generating random positions for iV nodes, we have complete 
information of <&(•). Using (2.4), we compute the communication range, and then map 
the STDMA network $(•) to the communication graph Q c {-)- Once the schedule *&(•) 
is computed by every algorithm, we know its schedule length \C\. For a given set of 
parameters (N,L,R C ), we calculate the average schedule length by averaging \C\ over 
1000 randomly generated networks. Keeping all other parameters fixed, we observe 
the effect of increasing the number of nodes on the average schedule length. In our 
experiments, we compare the performance of the following algorithms: 

• ArboricalLinkSchedule [16], 

• Proposed ALSReuseColors. 



L= 750.0 m, R = 100.0 m 

c 




Q| 1 1 1 1 1 1 1 1 1 

1 00 110 120 130 140 150 160 170 180 190 200 

number of nodes 

Figure 3.1: Schedule length vs. number of nodes. 

We assume that P = 10 mW, (3 = A, N = -90 dBm and 7 C = 20 dB. From (2.4), 
we obtain R c = 100 m. We assume that L = 750 m, and vary the number of nodes from 
100 to 200 in steps of 10. Figure 3.1 plots the average schedule length vs. number of 
nodes for both the algorithms. 
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For both the algorithms, we observe that average schedule length increases almost 
linearly with the number of nodes. The average schedule length of ALSReuseColors is 
about 23% lower than that of ArboricalLinkSchedule. 

Note that an increase in the number of nodes in a given geographical area leads to 
an increase in the number of edges incident on a vertex and a subsequent increase in 
the number of oriented graphs. ArboricalLinkSchedule, which is based on using a fresh 
set of colors for each oriented graph, requires increasingly higher number of colors to 
color the communication graph compared to ALSReuseColors. Consequently, the gap 
between the average schedule lengths increases with number of nodes in Figure 3.1. 

3.1.2 Analytical Results 

We now derive upper bounds on the running time (computational) complexity of the 
ALSReuseColors algorithm. With respect to the communication graph C? C (V,£), let: 

e = number of edges, 

v = number of vertices, 

p = maximum degree of any vertex, 

9 = thickness of the graph 

:= minimum number of planar graphs into which the undirected equivalent of Q c {-) 
can be partitioned, 

to = maximum number of neighbors with lower labels (for any vertex) . 

Recall that the modified algorithm partitions the communication graph Q c { ) into 
oriented graphs 7\, T 2 , . . . , T k , and colors the oriented graphs in that order. 7\ is termed 
as the first oriented graph, while any oriented graph Tj, where 2 ^ j ^ k, is termed as 
a subsequent oriented graph. 

Lemma 3.1.1. Suppose that each vertex of the first oriented graph T\ has at most uj 
neighbors with lower labels. Then, T x may be colored using no more than 0{ujp) colors. 

Proof. This is similar to the proof of Lemma 3.1 in [16]. ■ 

Lemma 3.1.2. Any subsequent oriented graph Tj, where 2 ^ j ^ k, can be colored using 
no more than 0(p 2 ) colors. 
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Figure 3.2: Potential conflicting edges when coloring edge (u,v). 

Proof. We prove the lemma for an out-oriented graph. A similar proof holds for an 
in-oriented graph. Let Q c be partitioned into edge-disjoint oriented graphs T l5 . . . ,TV 
Consider the coloring of edge (u, v) in j th oriented graph Tj, where 2 ^ j $C k, as shown 
in Figure 3.2. Now, edges of previously colored oriented graphs Ti, . . . , Tj-\ must also 
be considered for potential edge conflicts with edge (u,v) of Tj. Define 

j 



S\ := : (v,x) G [_J?i and (v,x) is coloredj, 

i=i 
j 

5 2 := x) : (u,x) G [J Tj and (w, x) is coloredj, 

i=i 
j 

5 3 ■= \{y,Xi) ■ (y,Xi) G [Jr, and (y,Xi) is colored and (u,^) G £ c j>, 

i=i 
i 

54 := : (zj,y) G (Jt,; and (xj,y) is colored and G £ c j. 

i=i 

Any edge which can cause a primary edge conflict with (u,v) must belong to Si or S 2 - 
Also, any edge which can cause a secondary edge conflict with (it, t>) must belong to £3 
or S4. Let rii = \Si\ for i = 1,2,3,4. The lemma reduces to proving that ni + n 2 + n 3 + n 4 
is 0(p 2 ). 

By definition of maximum vertex degree, rii ^ p — 1 and n 2 ^ p — 1. Thus, rii + n 2 
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is O(p). For the computation of 113, we must also consider secondary edge conflicts with 
edges of previously colored oriented graphs, as shown in Figure 3.2. The worst-case 
value of 113 is (p — l)(p — 1). Thus, 113 is 0(p 2 ). Similarly, by considering secondary edge 
conflicts with edges of previously colored oriented graphs, it follows that n 4 is 0(p 2 ). 
Finally, ri\ + n 2 + n 3 + is 0(p 2 ). ■ 

Lemma 3.1.3. For the first oriented graph T l7 the running time of Phase 3 of AL- 
SReuseColors is 0(yujp). 

Proof. This is similar to the proof of Lemma 3.2 in [16]. ■ 

Lemma 3.1.4. For any subsequent oriented graph Tj, where 2 ^ j ^ k, the running 
time of Phase 3 of ALSReuseColors is 0{vp 2 ). 

Proof. From Lemma 3.1.2, for any subsequent oriented graph Tj, the size of the set of 
conflicting colors (C\ U C2) of function NCEReuseColors is 0(p 2 ). Thus, determining a 
new color for an edge in Phase 3 of ALSReuseColors takes 0(p 2 ) steps. Since this is 
done for every label and hence for every vertex, it follows that the overall running time 
of Phase 3 of ALSReuseColors is 0(vp 2 ). ■ 

Theorem 3.1.5. For an arbitrary graph of thickness 9 and maximum degree p, AL- 
SReuseColors has a running time of 0(evlogv + vOp 2 ). 

Proof. The running time of the labeler function is 0(e + v logv) using a Fibonacci Heap 
[51]. The partitioning method of [49] results in a decomposition of a graph of thickness 
9 into at most 69 oriented graphs in time O(evlogv). Thus, k ^ 69. From Lemma 3.2 
in [16], the first oriented graph Ti can be colored in time 0{vujp). However, consider 
the coloring of j th oriented graph Tj, where 2 ^ j ^ k. From Lemma 3.1.4, Tj can be 
colored in time 0{vp 2 ). Hence, the for loop of ALSReuseColors runs in time 0(v9p 2 ). 
Therefore, the overall running time of ALSReuseColors is 0(e + v \ogv + ev log v + v9p 2 ). 
Since e + v logv < ev logv holds for any directed graph Q c {-) that models a wireless mesh 
network, the overall running time of ALSReuseColors simplifies to 0(ev log v+v9p 2 ). ■ 

3.1.3 Discussion 

In this section, we have considered an STDMA wireless network with uniform transmis- 
sion power at all nodes and presented an algorithm for point to point link scheduling 
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under the protocol interference model. The proposed algorithm, which is a modification 
of the ArboricalLinkSchedule algorithm in [16], models the network by a communication 
graph, partitions the communication graph into edge-disjoint oriented graphs and colors 
each oriented graph successively. However, unlike [16], we reuse colors from previously 
colored oriented graphs to color the current oriented graph. The proposed algorithm re- 
sults in around 26% lower schedule length than that of [16], albeit at the cost of slightly 
higher computational complexity 2 . Since schedules are constructed only once offline and 
then used by the network for a long period of time, our approach has the potential of 
providing higher long-term network throughput. 

For the rest of this chapter, we consider point to point link scheduling under the 
physical interference model. The algorithm developed in this section will be further 
refined to design a link scheduling algorithm in the next section. 

3.2 A High Spatial Reuse Link Scheduling Algo- 
rithm 

In this section, we propose a point to point link scheduling algorithm based on the com- 
munication graph model of an STDMA wireless network as well as SINR computations. 
We adopt spatial reuse as the performance metric, which has been motivated in Sec- 
tion 2.6. We compare the performance of the proposed algorithm with link scheduling 
algorithms which utilize a communication graph model of the network. We show that 
the proposed algorithm achieves higher spatial reuse compared to existing algorithms, 
without any increase in computational complexity. 

3.2.1 Problem Formulation 

Our system model and notation are exactly as described in Section 2.2. A link schedule 
is feasible if it satisfies the following conditions: 

1. Operational constraint (2.1). 

2 The computational complexity of ArboricalLinkSchedule is 0(ev\ogv + v6 2 p) [16]. 
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2. Range constraint: Every receiver is within the communication range of its intended 
transmitter, i.e., 

D(Uj, Uj) ^ R c V i = l,...,C V j = l,...,M t . (3.1) 

A link schedule is exhaustive if every pair of nodes which are within commu- 
nication range occur exactly twice in the link schedule, once with one node being the 
transmitter and the other node being the receiver, and during another time slot with 
the transmitter-receiver roles interchanged. Mathematically, 

c c 
D(j, k) < R c ^ j -> k e[jSi and k -> j e [J Si V 1 < j < k < N. (3.2) 

i=l i=l 

Our aim is to design a low complexity conflict-free STDMA point to point link 
scheduling algorithm that achieves high spatial reuse, where spatial reuse is given by 
(2.14). We only consider STDMA link schedules which are feasible and exhaustive 3 . 
Thus, our schedules satisfy (2.1), (2.7), (3.1) and (3.2). 

3.2.2 Motivation 

We briefly describe the essential features of STDMA link scheduling algorithms. An 
STDMA link scheduling algorithm is equivalent to assigning a unique color to every 
edge in the communication graph, such that transmitter-receiver pairs corresponding 
to communication edges with the same color are simultaneously active in a particular 
time slot, as described in Section 2.2.1. The core of a typical link scheduling algorithm 
consists of the following functions: 

1. An order in which communication edges are considered for coloring. 

2. A function which determines the set of all existing colors which can be assigned 
to the edge under consideration without violating the problem constraints. 

3. A BestColor rule to determine which color to assign to the edge under considera- 
tion. 

3 The set of edges in Q c {-) to be scheduled is determined by a routing algorithm. For simplicity, we 
only consider exhaustive schedules, i.e., schedules which assign exactly one time slot to every directed 
edge in Q c {-). 
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The second function considers only operational and range constraints in link scheduling 
algorithms based on the protocol interference model (equivalently, based on the com- 
munication graph). However, in the link scheduling algorithm that we propose, SINR 
constraints are also taken into account. 

Algorithms based on the protocol interference model are inadequate to design ef- 
ficient link schedules. This is because the communication graph Q C {V,£ C ) is a crude 
approximation of $(•). Even the two-tier graph Q(V, £ c U £j), which is a better approxi- 
mation of $(•), leads to low network throughput, as argued in Section 2.3. On the other 
hand, from $(•) and £ c (-), one can exhaustively determine the STDMA schedule which 
yields the highest spatial reuse. However, this is a combinatorial optimization problem 
of prohibitive complexity (0(\£ c \^)) and is thus computationally infeasible. 

To overcome these problems, we propose a new algorithm for STDMA link scheduling 
under the realistic physical interference model. Our algorithm is based on the commu- 
nication graph model G C (V,£ C ) as well as SINR computations. Motivated by techniques 
from matroid theory [52] , we develop a computationally feasible algorithm with demon- 
strably high spatial reuse. The essence of our algorithm is to partition the set of com- 
munication edges into subsets (forests) and color the edges in each subset sequentially. 
The edges in each forest are considered in a random order for coloring, since random- 
ized algorithms are known to outperform deterministic algorithms, especially when the 
characteristics of the input are not known a priori [53]. 

A similar matroid-based network partitioning technique is used in [54] to gener- 
ate high capacity subnetworks for a distributed throughput maximization problem in 
wireless mesh networks. Techniques from matroid theory have also been employed to 
develop efficient heuristics for NP-hard combinatorial optimization problems in fields 
such as distributed computer systems [55] and linear network theory [56]. 

3.2.3 ConflictFreeLinkSchedule Algorithm 

We call the proposed point to point link scheduling algorithm as ConflictFreeLinkSched- 
ule (CFLS). The algorithm considers the communication graph Q C {V,£ C ) and SINR con- 
ditions and is explained in Algorithm 4. 

In Phase 1, we label all the vertices randomly. Specifically, if Q c (-) has v vertices, we 
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perform a random permutation of the sequence (1,2, ... ,v) and assign these labels to 
vertices with indices 1,2, ... ,v respectively. L(u) denotes the label assigned to vertex 
u. 

In Phase 2, the communication graph Q c {-) is decomposed into what are called out- 
oriented and in-oriented graphs T\, T 2 , . . . ,T k [16]. Each Tj is a forest and every edge of 
Gc(-) is in exactly one of the Tj's. This decomposition is achieved by partitioning graph 
G c (-), the undirected equivalent of £/ c (-)> i n ^° undirected forests. The number of forests 
can be minimized by using techniques from Matroid theory ([49], A;- forest problem). 
However, this optimal decomposition requires extensive computation. Hence, we adopt 
the faster albeit non-optimal approach of using successive breadth first searches [50] 
to decompose G c (-) into undirected forests. Each undirected forest is further mapped 
to two directed forests. In one forest, the edges in every connected component point 
away from the root and every vertex has at most one incoming edge, thus producing an 
out-oriented graph. In the other forest, the edges in every connected component point 
toward the root and every vertex has at most one outgoing edge, thus producing an 
in-oriented graph. 

In Phase 3, the oriented graphs are considered sequentially. For each oriented graph, 
vertices are considered in increasing order by label and the unique edge associated with 
each vertex is colored using the FirstConflictFreeColor (FCFC) function. 

The FCFC function is explained in Algorithm 5. For the edge under consideration 
x, it discards any color that has an edge with a primary conflict with x. Among the 
residual set of colors, we choose the first color such that the resulting SINRs at the 
receiver of x and the receivers of all co-colored edges are no less than the communication 
threshold j c . If no such color is found, we assign a new color to x. Hence, this function 
guarantees that the ensuing schedule is conflict-free. 

3.2.4 Performance Results 
Simulation Model 

In the simulation experiments, the location of every node is generated randomly, using 
a uniform distribution for its X and Y coordinates, in the deployment area. For a fair 
comparison of our algorithm with the Truncated Graph-Based Scheduling Algorithm 
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Algorithm 4 ConflictFreeLinkSchedule (CFLS) 



1: input: STDMA network $(•), communication graph Q c {-) 

2: output: A coloring C : £ c — > {1, 2, . . .} 

3: label the vertices of Q c randomly {Phase 1} 

4: use successive breadth first searches to partition Q c into oriented graphs T i; 1 ^ i ^ k 

{Phase 2} 

5: for % <— 1 to A; do {Phase 3 begins} 

6: for j <— 1 to n do 

7: if Tj is out-oriented then 

8: let x = (s, d) be such that L(d) = j 

9: else 

10: let x = (s, d) be such that L(s) = j 

11: end if 

12: C(x) <— FirstConflictFreeColor(a;) 

13: end for 

14: end for{Phase 3 ends} 



Algorithm 5 integer FirstConflictFreeColor(x) 
1: input: STDMA network $(•), communication graph Q c { ) 

2: output: A conflict-free color 

3: C <— set of existing colors 

4: C c <— {C{h) : h G S c , h is colored, x and ft, have a primary edge conflict} 

5: C c f = C\C C 

6: for % 1 to |C C /| do 

7: r <— color in C c / 

8: {h:he £ c , C(h) = r} 

9: C(x) <— r 

10: if SINR at all receivers of Ei U {x} exceed 7 C then 

11: return r 

12: end if 

13: end for 

14: return \C\ + 1 
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(TGSA) [32], we assume that the deployment region is a circular region of radius R. 
Thus, if (Xj,Yj) are the Cartesian coordinates of j th node, j = l,...,iV, then Xj ~ 
U[-R, R] and Yj ~ U[-R, R] subject to X? + Yf < i? 2 . Equivalently, if (i2 i? 9-,) are the 
polar coordinates of j th node, then i? 2 ~ U[0,R 2 ] and Gj ~ U[0, 2ir]. After generating 
random positions for N nodes, we have complete information of $(•)• Using (2.4) and 
(2.5), we compute the communication and interference radii, and then map the network 
$(•) to the two-tier graph Q(V,£ C U Once the link schedule is computed by an 
algorithm, a is computed using (2.14). System parameters are chosen based on their 
prototypical values in wireless mesh networks [42] . For a given set of system parameters, 
we calculate the average spatial reuse by averaging a over 1000 randomly generated 
networks. Keeping all other parameters fixed, we observe the effect of increasing the 
number of nodes N on the average spatial reuse. 

In our experiments, we compare the performance of the following algorithms: 

• ArboricalLinkSchedule (ALS) [16], 

• Truncated Graph-Based Scheduling Algorithm 4 (TGSA) [32], 

• GreedyPhysical (GP) [27], 

• Proposed ConflictFreeLinkSchedule (CFLS). 

Performance Comparison under Path Loss Model 

In the first set of experiments (Experiment 1), we assume that R = 500 m, P — 10 

mW, (3 = 4, N = -90 dBm, lc = 20 dB and 7i = 10 dB [42]. Thus, R c = 100 m and 

Ri = 177.8 m. We vary the number of nodes from 30 to 110 in steps of 5. Figure 3.3 

plots the average spatial reuse vs. number of nodes for all the algorithms. 

In the second set of experiments (Experiment 2), we assume that R = 700 m, P — 15 

mW, p = A, N = -85 dBm, 7c = 15 dB and 7i = 7 dB. Thus, R c = 110.7 m and 

4 In Truncated Graph-Based Scheduling Algorithm, for the computation of optimal number of trans- 
missions M*, we follow the method described in [32]. Since < £ < we assume that £ = 0.9999-^ 
and compute successive Edmundson-Madansky (EM) upper bounds [44] , [45] , till the difference between 
successive EM bounds is less than 0.3%. We have experimentally verified that only high values of £ lead 
to reasonable values for M* , whereas low values of £, say £ = 0.1-^r, lead to the extremely conservative 
value of M* = 1 in most cases. 
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Figure 3.3: Spatial reuse vs. number of nodes for Experiment 1. 



Ri = 175 A m. We vary the number of nodes from 70 to 150 in steps of 5. Figure 3.4 
plots the average spatial reuse vs. number of nodes for all the algorithms. 

For the ALS algorithm, we observe that spatial reuse increases very slowly with 
increasing number of nodes. 

For the TGSA algorithm, we observe that spatial reuse is 18-27% lower than that 
of ALS and 30-55% lower than that of GP. A plausible explanation for this behavior 
is as follows. The basis for TGSA is the computation of M*, the optimal number of 
transmissions in every slot [32]. M* is determined by maximizing a lower bound on the 
expected number of successful transmissions in a time slot. Since the partitioning of a 
maximal independent set of communication arcs into subsets of cardinality at most M* 
is arbitrary and not geography-based, there could be scenarios where the transmissions 
scheduled in a subset are in the vicinity of each other, resulting in moderate to high 
interference. In essence, maximizing this lower bound does not necessarily translate to 
maximizing the number of successful transmissions in a time slot. Also, due to its design, 
the TGSA algorithm yields higher number of colors compared to ALS and GP. 

Though the GP algorithm is based on communication graph and SINR conditions, it 
yields slightly lower spatial reuse than CFLS. A possible reason for this observation is as 
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Figure 3.4: Spatial reuse vs. number of nodes for Experiment 2. 

follows. The GP algorithm colors edges of the communication graph in the decreasing 
order of interference number. The interference number of edge e is the number of edges 
Ci such that, if (e, e^) are scheduled simultaneously, then the SINR threshold condition 
(2.7) is violated along one or both links. Edges with higher interference number tend to 
be located towards the center of the deployment region. Since these edges are colored 
first, a large number of colors are utilized in the initial stages of the algorithm, lead to 
potentially higher schedule length and lower spatial reuse. A better technique would be 
successively examine edges at the centre and the periphery, which is achieved by the 
partition technique employed by CFLS. 

For the proposed CFLS algorithm, we observe that spatial reuse increases steadily 
with increasing number of nodes and is about 15% higher than the spatial reuse of ALS, 
TGSA and GP. 



Performance Comparison under Realistic Conditions 

In a realistic wireless environment, channel impairments like multipath fading and shad- 
owing affect the received SINR at a receiver [14]. In this section, we compare the 
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performance of the ALS, TGSA, GP and CFLS algorithms in a wireless channel which 
experiences Rayleigh fading and lognormal shadowing. 

In the absence of fading and shadowing, the SINR at receiver r^j is given by (2.2). 
We assume that every algorithm (ALS, TGSA, GP and CFLS) considers only path loss 
in the channel prior to constructing the two-tier graph Q(V, £ c U £») and computing the 
link schedule. 

However, for computing the average spatial reuse of each algorithm, we take into ac- 
count fading and shadowing channel gains between each pair of nodes. More specifically, 
for computing the spatial reuse using (2.14), the (actual) SINR at receiver r it j is given 
by 

SINR r . = 7? ^^^ fej,n ' j)1 ° H/(tlJ ' r ' j) , (3.3) 

,J ^ + E^^^n^,n,)io^^)' 

where random variables V(-) and W(-) correspond to channel gains due to Rayleigh 
fading and lognormal shadowing respectively. We assume that {V(k, ^ k, I ^ N, k ^ 
1} are independent and identically distributed (i.i.d.) random variables with probability 
density function (pdf) [10] 

fv(v) = \e4 u {v), (3.4) 
°V 

where u(-) is the unit step function. Also, {W(k, l)\l ^ k,l ^ N,k ^ /} are assumed to 
be i.i.d. zero mean Gaussian random variables with pdf [57] 

1 -TO 2 

fw(w) = —^ e *%. (3.5) 
v 



Random variables V(-) and W(-) are independent of each other and also independent of 
the node locations. 

The simulation model and experiments are exactly as described before. In the simu- 
lations, we assume ay = = 1. For Experiment 1, Figure 3.5 plots the average spatial 
reuse vs. number of nodes for all the algorithms. For Experiment 2, Figure 3.6 plots 
the average spatial reuse vs. number of nodes for all the algorithms. 

From Figures 3.3, 3.4, 3.5 and 3.6, we observe that spatial reuse decreases by 20- 
40% in a channel experiencing multipath fading and shadowing effects. A plausible 
explanation for this observation is as follows. Since the channel gains between every pair 
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Figure 3.5: Spatial reuse vs. number of nodes for Experiment 1 under multipath fading 
and shadowing channel conditions. 
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Figure 3.6: Spatial reuse vs. number of nodes for Experiment 2 under multipath fading 
and shadowing channel conditions. 
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of nodes are independent of each other, it is reasonable to assume that the interference 
power at a typical receiver remains almost the same as in the non-fading case. This 
is because, even if the power received from few unintended transmitters is low, the 
power received from other unintended transmitters will be high (on average); thus the 
interference power remains constant. Consequently, the change in SINR is determined by 
the change in received signal power only. If the received signal power is higher compared 
to the non-fading case, the transmission is anyway successful and spatial reuse remains 
unchanged (see (2.14)). However, if the received signal power is lower, the transmission 
is now unsuccessful and spatial reuse decreases. Hence, on average, the spatial reuse 
decreases. 

Finally, from Figures 3.5 and 3.6, we observe that the proposed CFLS algorithm 
achieves 5-17% higher spatial reuse than the ALS and GP algorithms and 40-80% higher 
spatial reuse than the TGSA algorithm, under realistic wireless channel conditions. 

3.2.5 Analytical Results 

In this section, we derive upper bounds on the running time (computational) complexity 
of ConflictFreeLinkSchedule algorithm. We use the following notation with respect to 
the communication graph Q C (V,S C ): 

e = number of communication edges, 
v = number of vertices, 
9 = thickness of the graph 

:= minimum number of planar graphs into which the undirected equivalent of Q c {-) 
can be partitioned. 

Before we prove our results, it is instructive to observe Figure 3.7, which shows the 
variation of 9 and e with v for the two experiments described in Section 3.2.4. Since 
determining the thickness of a graph is NP-hard [58], each value of 9 in Figure 3.7 is 
an upper bound on the actual thickness based on the number of forests into which the 
undirected equivalent of the communication graph has been decomposed using successive 
breadth first searches. We observe that the graph thickness increases very slowly with 
the number of vertices (^Co), while the number of edges increases super-linearly with 
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R = 500 m, P = 1 mW, p = 4, N Q = -90 dBm, y c = 20 dB, y. = 1 dB 




number of nodes 



(a) Experiment 1 



R = 700 m, P = 15 mW, R = 4, N = -85 dBm, y = 15 dB, v = 7 dB 




number of nodes 



(b) Experiment 2 

Figure 3.7: Comparison of thickness and number of edges with number of vertices. 
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the number of vertices. 

Lemma 3.2.1. An oriented graph T can be colored using no more than 0{v) colors 
using ConflictFreeLinkSchedule. 

Proof. Since an oriented graph with v vertices has at most v edges, the edges of T can 
be colored with at most v colors. ■ 

Lemma 3.2.2. For an oriented graph T , the running time of ConflictFreeLinkSchedule 
is 0(v 2 ). 

Proof. Assuming that an element can be chosen randomly and uniformly from a finite 
set in unit time ([53], Chapter 1), the running time of Phase 1 can be shown to be 
0(v). Since there is only one oriented graph, Phase 2 runs in time 0(v). In Phase 3, 
the unique edge associated with the vertex under consideration is assigned a color using 
FirstConflictFreeColor. From Lemma 3.2.1, the size of the set of colors to be examined 
\C C U C c f\ is 0(v). In FirstConflictFreeColor, the SINR is checked only once for every 

\C I 

colored edge in the set (Ji=i an d at most v times for the edge under consideration 
x. With a careful implementation, FirstConflictFreeColor runs in time 0(v). So, the 
running time of Phase 3 is 0(v 2 ). Thus, the total running time is 0(v 2 ). ■ 

Theorem 3.2.3. For an arbitrary graph Q, the running time of ConflictFreeLinkSched- 
ule is 0(ev\ogv + ev9). 

Proof. Assuming that an element can be chosen randomly and uniformly from a finite 
set in unit time [53], the running time of Phase 1 can be shown to be 0(v). For Phase 
2, the optimal partitioning technique of [49] based on Matroids can be used to partition 
the communication graph Q c into at most 69 oriented graphs in time O(evlogv). Thus, 
k ^ 69 holds for Phase 3. From Lemma 3.2.2, it follows that the first oriented graph Ti 
can be colored in time 0(v 2 ). However, consider the coloring of j th oriented graph Tj, 
where 2 ^ j ^ k. When coloring edge x from Tj using FirstConflictFreeColor, conflicts 
can occur not only with the colored edges of Tj, but also with the edges of the previously 
colored oriented graphs Ti,T 2 , . . . , Tj_i. Hence, the worst-case size of the set of colors 
to be examined \C C U C c /| is 0(e). Note that in FirstConflictFreeColor, the SINR is 
checked only once for every colored edge in the set Ui=i &i an d at most e times for the 
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edge under consideration x. With a careful implementation, FirstConflictFreeColor runs 
in time 0(e). Hence, any subsequent oriented graph Tj can be colored in time 0(ev). 
Thus, the running time of Phase 3 is 0(ev9). Therefore, the overall running time of 
ConflictFreeLinkSchedule is 0(ev\ogv + evO). ■ 

3.2.6 Discussion 

In this section, we have developed ConflictFreeLinkSchedule, a point to point link 
scheduling algorithm for an STDMA wireless mesh network under the physical interfer- 
ence model. The performance of the proposed algorithm is superior to those of existing 
link scheduling algorithms for STDMA wireless networks with uniform power assign- 
ment. A practical experimental modeling shows that, on average, the proposed algorithm 
achieves 20% higher spatial reuse than the ArboricalLinkSchedule [16], GreedyPhysical 
[27] and Truncated Graph-Based Scheduling [32] algorithms. Since link schedules are 
constructed offline only once and then used by the network for a long period of time, these 
improvements in performance directly translate to higher long-term network throughput. 

The computational complexity of ConflictFreeLinkSchedule is comparable to the 
computational complexity of ArboricalLinkSchedule and is much lower than the com- 
putational complexity of GreedyPhysical and Truncated Graph-Based Scheduling al- 
gorithms. Thus, in cognizance of spatial reuse as well as computational complexity, 
ConflictFreeLinkSchedule appears to be a good candidate for efficient STDMA link 
scheduling algorithms. 



Chapter 4 

Point to Point Link Scheduling 
based on SINR Graph Model 

In this chapter, we propound a somewhat different approach for point to point link 
scheduling in an STDMA wireless network under the physical interference model. This 
approach is based on SINR graph representation of the network wherein weights of edges 
correspond to interferences between pairs of nodes and weights of vertices correspond to 
normalized noise powers at receiving nodes. We develop a novel link scheduling algorithm 
with polynomial time complexity and improved performance in terms of spatial reuse. 

The rest of the chapter is organized as follows. We motivate our SINR graph approach 
in Section 4.1. We describe the proposed link scheduling algorithm and provide an 
illustrative example in Section 4.2. We prove the correctness of the algorithm and 
derive its computational complexity in Section 4.3. The performance of the proposed 
algorithm is compared with existing link scheduling algorithms in Section 4.4. We discuss 
the implications of our work in Section 4.5. 

4.1 Motivation 

The system model, notation and problem formulation are exactly as described in Section 
3.2.1. Specifically, we seek a low complexity conflict-free point to point link scheduling 
algorithm that achieves high spatial reuse. 

In general, for the STDMA wireless network $(•), the set of links to be scheduled 
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is determined by a routing algorithm. For simplicity, we only consider exhaustive link 
schedules, i.e., we consider uniform load on all links. 

Note that for point to point link schedules that are conflict-free, i.e., for link schedules 
that satisfy (2.7), the equation for spatial reuse (2.14) reduces to 

g 

Spatial Reuse = o = — , (4.1) 

where e denotes the number of directed edges in the communication graph Q C (V,S C ) 
and C denotes the number of slots in the link schedule. Therefore, for conflict-free link 
schedules, maximizing spatial reuse is equivalent to minimizing the number of colors, 
i.e., minimizing the schedule length. 

To the best of our knowledge, there is no known polynomial time algorithm that 
determines a provably optimal schedule (minimum length schedule) for an STDMA 
wireless network with constrained transmission power. Hence, the only recourse is to 
devise heuristics and show their efficiency theoretically and experimentally. Towards 
this end, we propose a heuristic based on an SINR graph representation of the network. 

Consider any directed graph G(V, E), where V is the set of vertices and E is the set 
of edges. The line graph of G(V, E) is the graph G'(V, E') whose vertices are the edges 
of G(-), i.e., V — E [35]. The SINR graph that we consider in this chapter is analogous 
to the concept of line graph in [35]. However, unlike the line graph, we assume that the 
SINR graph is a complete graph, i.e., for any two distinct vertices v •, v'j 6 V, there is a 
directed edge from v[ to v'j in E' . 

The crux of the proposed link scheduling algorithm can be understood by revisit- 
ing the condition for successful packet reception under the physical interference model 
(Equation 2.6), i.e., 

p 

Rearranging the terms in (4.2), we obtain 

—D (^ >rw ) + ^ 7e __ 5 < 1. (4.3) 

Dropping time slot index i for clarity, we obtain the "equivalent" SINR threshold con- 
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dition 

^hpP( t . r ) + f /^'^ < 1 (44) 

where tj, rj and M can be interpreted as j th transmitter, j th receiver and number of 
concurrent transmissions, respectively, in a given time slot. The terms appearing in (4.4) 
correspond to vertex and edge weights in a special graph representation of the STDMA 
network, termed as SINR graph. This idea will be elucidated further in Section 4.2.1. 

4.2 SINRGraphLinkSchedule Algorithm 

In this section, we explain the proposed link scheduling algorithm based on SINR graph 
representation of the STDMA network. We provide an illustrative example to elucidate 
the intricacies of the proposed algorithm. 

4.2.1 Description 

The proposed link scheduling algorithm under the physical interference model is 
SINRGraphLinkSchedule (SGLS), which considers the communication graph £/ c (V, £ c ). 

First, we construct a directed complete SINR graph G'(V',£') that has the edges of 
Gc(-) as its vertices, i.e., V' = £ c . Let the edges of Gc(-) and the corresponding vertices of 
£/'(•) be labeled 1,2, ... ,e. Let U and r\ denote the transmitter and receiver respectively 
of edge i in Gc(-)- For any two edges i and j in graph G c {-), the interference weight 
function Wij is defined as: 

{1 if i and j have a common vertex, 

^m^W otherwise - 

The interference weight function indicates the interference energy at r,j due to trans- 
mission from ti to rj scaled with respect to the signal energy of tj at rj. Note that the 
interference weight function appears as a summand in the equivalent SINR threshold 
condition (4.4). 

We then compute the co-schedulability weight function w'. For any two edges i and 
j in Gc(-), the weight of edge e\, in G'(-) is given by w'^ = max{0, 1 — Since 
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Algorithm 6 SINRGraphLinkSchedule (SGLS) 



Input: Communication graph Q C (V,£ C ), 7 C , iV , P 
Output: A coloring C: £ c — > {1, 2, . . .} 
V'«-£ c 

Construct the directed complete graph Q'(V',£') 



for all e', G 5' do 



if edges i and j have a common vertex in £/ c (-) then 



Wij 



else 



Wij 



''D(U, rj f 



end if 
end for 

for all e', G 5' do 

w'ij <— max{0, 1 - lUj.,} 
end for 

for all v'j G V do 

N{v>) ^ ^D^rtf 
end for 

p - 0; K c <- V 
while V uc ^ do 

p <— p + 1; choose t> ' G V^ c randomly 
<V) - P5 K c - V: c \ {t/}; V; - {t/}; V - 1 
while ^ = 1 and V uc + and max^^ E^'ev^ w iv + «W > do 
for all u' G K c such that J2x>ev w 'x'w + K>x> > do 

£><-l 

for all ^ G do 

g <- 
end if 
end for 

if q= 1 and £„, 6V , > |V;| + A/V) " 1 then 

coo - p; v; <- v; u {«'}; k c <- k c \ {«'} 

else 

g <- 
end if 
end for 
if q = then 

V> <- 
end if 
end while 
end while 
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and Wji represent interferences among links i and j in the STDMA network $(•), w'^ 
and intuitively represent the co-schedulability of links i and j in $(•) (equivalently, 
co-schedulability of vertices z and j in G'(-))- For example, if -1% is greater than or equal 
to 1, then the interference at the receiver of link j from the transmitter of link i is very 
high and these links cannot be scheduled simultaneously. This will result in w'^ being 
equal to indicating that vertices i and j in Q'(-) are not co-schedulable. On the other 
hand, if is slightly greater than (0 < <C 1), w'^ will be slightly less than 1 
indicating that the vertices i and j in G'(-) are co-schedulable. Note that for the SINR 
graph G'(-), the weight of an edge refers to the value of co-schedulability function for 
that edge. 

Next, we determine the normalized noise power at the receiver of each link of $(■). 
This is tantamount to computing the normalized noise power for each edge of Gc(-), i-e., 
at each vertex of G'{)- Note that the normalized noise power function appears as a term 
in the equivalent SINR threshold condition (4.4). 

Our objective is to color the vertices of G'(-) (equivalently, edges of G c {-)) using 
minimum number of colors under the physical interference model, i.e., subject to the 
condition that the SINR at the receiver of every link in $(•) is no less than the commu- 
nication threshold 7 C . Equivalently, for any V' cc C V, the coloring of all vertices v[ € V' cc 
with the same color is defined to be feasible if 

^ HV^ (4.5) 



^0 + E^ c \{<} D(t M f 



In the SINR graph G'(-), this condition translates to the sum of weights of edges incoming 
to a vertex from all co-colored vertices being greater than the sum of the number of 
remaining co-colored vertices and the normalized noise power minus a constant factor 
(unity); this will be proved in Theorem 4.3.1. 

Finally, we color vertices of G'(-), i-e., edges of G c {-), according to the following 
procedure. Let V' uc denote the set of uncolored vertices of G'{-). Initially, V' uc includes 
all vertices of G'{-)- First, we choose a vertex randomly from V' uc . This is assigned a new 
color, say p. Then, we consider every vertex v! from V' uc such that the sum of weights of 
all the edges between u' and the vertices colored with p is positive. Next, for each vertex 
colored with p, we check if the sum of weights of all incoming edges is greater than the 
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sum of the number of vertices colored with p and the normalized noise power at that 
vertex minus a constant factor (unity). If this inequality is satisfied, we further check 
if the sum of weights of all edges incoming to u is greater than the sum of the number 
of vertices colored with p and the normalized noise power at u' minus unity. If this 
inequality is also satisfied, then vertex u' is colored with p. If any of these inequalities 
are not satisfied, vertex v! is colored with a new color. The algorithm exits when all the 
vertices are colored. The pseudocode of the algorithm is provided in Algorithm 6. 



4.2.2 Example 

Consider the STDMA wireless network $(•) whose deployment is shown in Figure 4.1. It 
consists of four labeled nodes whose coordinates (in meters) are 1 = (—40, 5), 2 = (0, 0), 
3 = (95,0) and 4 = (135,0). We use typical values of system parameters in wireless 
networks [42]. These values are shown in Table 4.1, which lead to R c = 100 m. 



l = (-40,5) 



Y 



X 
— - 



2 = (0,0) 



3 = (95,0) 4 = (135,0) 



Figure 4.1: An STDMA wireless network with four nodes. 



Parameter 


Symbol 


Value 


transmission power 


P 


10 mW 


path loss exponent 


P 


4 


noise power spectral density 


N 


-90 dBm 


communication threshold 


7c 


20 dB 



Table 4.1: System parameters for the STDMA network shown in Figure 4.1. 

The communication graph model of the STDMA network is shown in Figure 4.2. 
The communication graph Q C (V, £ c ) consists of four vertices and six directed edges. The 
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vertex and edge sets are given by 

V = {vi,V 2 ,V3,V4}, (4.6) 

£ c = {(1,2), (2,1), (2,3), (3,2), (3,4), (4,3)}. (4.7) 




Figure 4.2: Communication graph model of STDMA network described by Figure 4.1 
and Table 4.1. 

The SINR graph model of the communication graph Q C (V, S c ) is shown in Figure 4.3. 
The SINR graph Q'(V',S') is a complete graph and consists of six vertices and thirty 
directed edges. The vertex set of the SINR graph is given by 

V = {(1,2), (2,1), (2,3), (3,2), (3,4), (4,3)}. (4.8) 

The edge set £' of the SINR graph is enumerated in Table 4.2, along with the interference 
weight function Wij and co-schedulability weight function w'^ for each edge i — > j e <?'(•)• 
The normalized noise powers at vertices of the SINR graph are enumerated in Table 4.3. 

The truncated SINR graph G' t {V,£' t ) is shown in Figure 4.4. The truncated SINR 
graph consists of all vertices of the SINR graph and only those edges whose co-schedulability 
weight function is positive, i.e., S[ = : i,j G £ c and w'^ > 0}. The values of the 

co-schedulability weight functions for all edges and the normalized noise powers at all 
vertices are also shown in the figure. We use the truncated SINR graph to explain the 
SGLS algorithm, since edges having zero weight in the SINR graph do not play any role 
in the SGLS algorithm. Note that, in the truncated SINR graph, the weight of an edge 
refers to the value of the co-schedulability weight function for that edge. 

Initially, the set of uncolored vertices is V uc = {(1,2), (2, 1), (2, 3), (3, 2), (3, 4), (4, 3)}. 
In the first iteration, we randomly choose v' = (1,2) and assign it Color 1 (say, red). So, 
C(l,2) = 1. The set of uncolored vertices is V uc = {(2, 1), (2, 3), (3, 2), (3, 4), (4, 3)} and 
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Figure 4.3: SINR graph model of communication graph shown in Figure 4.2. 



0.0264 0.8145 




0.0256 0.0256 



Figure 4.4: Truncated SINR graph derived from SINR graph shown in Figure 4.3 and 
weight values given in Tables 4.2 and 4.3. 
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Edge % — > j of 

SINR graph Q'{V,£') 



w. 



Edge 2 — > j of 

SINR graph Q'{V,£') 



w 



w 



u 



1,2 



2, 1 



3,2 



1,2 



2,1 



1,2 



0.2725 



0.7275 



3,4 



3.2420 



2,1 



3,2 



2,1 



0.7707 



0.2293 



3,4 



0.7928 



0.2072 



2,1 



3.1430 



4,3 



0.2811 



0.7189 



2,3 



0.7950 



0.2050 



3,4 



0.7686 



0.2314 



2,3 



4,3 



3,2 



1,2 



3,2 







4,3 



3,2 



Table 4.2: Interference and co-schedulability weight functions for edges of SINR graph 
shown in Figure 4.3. 



Vertex v'j of 

SINR graph Q'{V,£') 


a/>;) 


(1,2) 


0.0264 


(2,1) 


0.0264 


(2,3) 


0.8145 


(3,2) 


0.8145 


(3,4) 


0.0256 


(4,3) 


0.0256 



Table 4.3: Normalized noise powers at vertices of SINR graph shown in Figure 4.3. 
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the set of vertices colored 1 is V' C1 = {(1,2)}. From the set of uncolored vertices V^ c , we 
consider every vertex u' such that the sum of weights of edges from the presently colored 
vertex (1, 2) to u' and from u' to the presently colored vertex is positive. From Figure 4.4, 
we obtain two candidates: v! = (3,4) and u' = (4,3). We first examine the candidate 
vertex (3, 4). We check if the weight of the edge from (3, 4) to the presently colored vertex 
(1,2) is no greater than the number of vertices colored with the present color (red) plus 
the normalized noise power at the colored vertex minus unity. Our calculations show 
that inequality holds (0 < 0.0264) and candidate vertex (3,4) cannot be assigned Color 
1. We next examine the candidate vertex (4, 3). We check if the weight of the edge from 
the candidate vertex to (1,2) is no greater than the number of vertices colored with the 
present color plus the normalized noise power at (1,2) minus unity. Our calculations 
show that inequality does not hold (0.2050 ^ 0.0264). Furthermore, we check if the 
weight of the edge from the presently colored vertex (1, 2) to the candidate vertex (4, 3) is 
greater than the number of vertices colored red plus the normalized noise power at (4, 3) 
minus unity. The inequality holds and hence the candidate vertex (4, 3) is assigned Color 
1 (red). So, C(4,3) = 1. The set of uncolored vertices is V' uc = {(2, 1), (2, 3), (3, 2), (3, 4)} 
and the set of vertices colored 1 is V' C1 = {(1, 2), (4, 3)}. Again, from the set of uncolored 
vertices V' uc , we consider every vertex u' such that the sum of weights of edges from the 
presently colored vertices {(1, 2), (4, 3)} to u' and from u' to the presently colored vertices 
is positive. From Figure 4.4, the candidate vertices are (2, 1) and (3,4). Consider the 
candidate vertex (2, 1). For every vertex v' c colored 1, we check if the sum of weights of 
edges from remaining co-colored vertices and the candidate vertex to the colored vertex 
is no greater than the number of co-colored vertices and the normalized noise power at 
the colored vertex minus unity. For the colored vertex (1,2), our calculations show that 
inequality holds (0.2050 ^ 1.0264). So, we discard (2,1), consider the next candidate 
vertex (3,4) and perform an analogous comparison with v' c = (1,2). Since inequality 
holds in this case too (0.2050 < 1.0264), we discard (3,4) and proceed to the next 
iteration. The set of vertices colored so far is shown in Figure 4.5. 

In the second iteration, we randomly choose v' = (2, 3) and assign it Color 2 (say, 
blue). So, C(2,3) = 2. The set of uncolored vertices is V uc = {(2, 1), (3,2), (3,4)} and 
the set of vertices colored 2 is V' C2 = {(2, 3)}. From the set of uncolored vertices V' uc , we 
consider every vertex u' such that the sum of weights of edges from the presently colored 
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0.0264 0.8145 




0.0256 0.0256 

Figure 4.5: Coloring of vertices of truncated SINR graph after first iteration of SGLS 
algorithm. 



vertex (2, 3) to u' and from u' to the presently colored vertex is positive. From Figure 
4.5, no such vertex exists. So, we proceed to the next iteration. The vertices colored so 
far are shown in Figure 4.6. 

In the third iteration, we randomly choose v' = (3,4) and assign it Color 3 (say, 
green). So, C(3, 4) = 3. The set of uncolored vertices is V' uc = {(2, 1), (3, 2)} and the set 
of vertices colored 3 is V' C3 = {(3,4)}. From the set of uncolored vertices V' uc , we consider 
every vertex u' such that the sum of weights of edges from the presently colored vertex 
(3, 4) to v! and from v! to the presently colored vertex is positive. From Figure 4.4, 
we obtain u' = (2, 1) as the only candidate vertex. Next, we check if the weight of the 
edge from the candidate vertex (2, 1) to the presently colored vertex (3,4) is no greater 
than the number of vertices colored with the present color (green) plus the normalized 
noise power at the colored vertex minus unity. Our calculations show that inequality 
does not hold (0.2293 ^ 0.0256). So, we further check if the weight of the edge from 
the presently colored vertex (3,4) to the candidate vertex (2,1) exceeds the number of 
vertices colored with the present color plus the normalized noise power at the candidate 
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0.0264 0.8145 




0.0256 0.0256 

Figure 4.6: Coloring of vertices of truncated SINR graph after second iteration of SGLS 
algorithm. 



vertex minus unity. Since the inequality holds (0.2072 > 0.0264), the candidate vertex 
(2, 1) is assigned Color 3 (green). So, C(2, 1) = 3. The set of uncolored vertices is 
V' uc = {(3,2)} and the set of vertices colored green is V' c = {(3,4), (2, 1)}. Next, from 
the set of uncolored vertices V' uc , we choose that uncolored vertex v! such that the sum 
of weights of edges from v! to the set of presently colored vertices {(3, 4), (2, 1)} and 
from {(3,4), (2,1)} to u' is positive. From Figure 4.6, no such vertex u' exists. So, we 
proceed to the next iteration. Figure 4.7 shows the set of vertices colored so far. 

In the fourth iteration, (3,2) is the only uncolored vertex. So, we choose v' = (3,2) 
and assign it Color 4 (say, pink). The set of vertices colored 4 is V' Ci = {(3,2)} and 
the set of uncolored vertices is V' uc = (p. So, the algorithm ends. The final coloring of 
vertices of the truncated SINR graph by SGLS algorithm is shown in Figure 4.8. 

The output of the SGLS algorithm is enumerated in Table 4.4 and is also shown 
pictorially in Figure 4.9. The resulting link schedule is denoted by \I/(<Si, S2, S3, S4, S5), 
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0.0264 0.8145 




0.0256 0.0256 

Figure 4.7: Coloring of vertices of truncated SINR graph after third iteration of SGLS 
algorithm. 

where 





= {1- 


+ 2, 4- 


-3} 


s 2 


= {2- 


-3}, 




S 3 


= {3- 


+ 4, 2- 


-1} 


<s 4 


= {3- 


-2}. 





Finally, we check if the link schedule enumerated in Table 4.4 is conflict-free, i.e., if 



Time slot 


Color 


Active (transmitter, receiver) pairs 


1 


red 


(1,2), (4,3) 


2 


blue 


(2,3) 


3 


green 


(3,4), (2,1) 


4 


pink 


(3,2) 



Table 4.4: Output of SGLS algorithm for STDMA network described by Figure 4.1 and 
Table 4.1. 
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0.0264 0.8145 




0.0256 0.0256 

Figure 4.8: Coloring of vertices of truncated SINR graph after complete execution of 
SGLS algorithm. 




Figure 4.9: Output of SGLS algorithm for STDMA network described by Figure 4.1 and 
Table 4.1. 
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the SINR threshold condition (2.7) is satisfied at every receiver for the STDMA network 
described by Figure 4.1 and Table 4.1. Only one transmitter-receiver pair is active 
during time slots 2 and 4. Since the receiver is within the communication range of its 
corresponding transmitter for each of these time slots, the SINR threshold condition 
is satisfied trivially for time slots 2 and 4. Two transmitter-receiver pairs are active 
during time slots 1 and 3. In time slot 1, the active transmitter-receiver pairs are (1,2) 
and (4,3). Our computations show that the SINRs at Receivers 2 and 3 are 20.85 dB 
and 21 dB, both of which exceed the communication threshold of 20 dB. In time slot 3, 
the active transmitter-receiver pairs are (2,1) and (3,4). Our computations show that 
the SINRs at Receivers 1 and 4 are 20.87 dB and 20.99 dB respectively, both of which 
exceed the communication threshold. This verifies that the SGLS algorithm yields a 
conflict-free link schedule for the network described by Figure 4.1 and Table 4.1. Note 
that, from (4.1), the spatial reuse of SGLS algorithm for this network is 1.5. 



4.3 Analytical Results 

In this section, we prove the correctness of the SGLS algorithm and derive its running 
time (computational) complexity. We follow the notation of Algorithm 6. 

Theorem 4.3.1. For any V cc CV' if 

E w kv[ > \K\ +M(v[) - 2 V v[ G V cc , (4.9) 
v' 2 ev^\{v[} 

then the coloring of all vertices of V' cc with the same color is feasible. 

Proof. Recall that w',, = or 1 — w v ' v > and that ^ w',, ^ 1. Suppose w',, = for 

v 2 1 2 1 v 2 I ' 3 1 

some v[,v 3 e V cc , v[ ^ v' 3 , then 

v' 2 ev^\{v[} t/ 2 ev^\{«i,t;j} 

= \K\{v[,v' 3 }\, 
= IVLI-2, 

which contradicts the hypothesis since M(v[) > 0. So, an edge connecting any two 
vertices in V' cc must have positive weight. Thus, < w' v i , ^ 1 V v[,v' 2 G V' cc , v[ ^ v' 2 . 



78 



Chapter 4. Point to Point Link Scheduling based on SINR Graph Model 



Equivalently, < 1 — w v 'yi < 1 V v' 1 ,v' 2 G V' cc , v[ ^ v 2 . If two vertices v[,v 2 G V' cc 
(equivalently, edges v[,v 2 G G c {')) have a common vertex in Q c (-), then w^tf = 1, which 
is a contradiction. So, no two vertices in V' cc have a common vertex in G c {-)- From the 
hypothesis, 
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0(4,4 ,»■„#) 



> 7c V 4 G V c . 



^O + E^eV^WI D(t v , ,r vt ' 

2 1 

Therefore, the SINR threshold condition (4.5) is satisfied at the receivers of all vertices 



of V cc . 



With respect to (w.r.t.) the communication graph Q C (V, £ c ), let: 

e = number of edges, 
v = number of vertices. 



Theorem 4.3.2. The running time complexity of SGLS algorithm is 0(e 2 ). 

Proof. \V'\ = |£ c | = e. Since Q'(-) is a directed complete graph, \£'\ = e(e — 1) = 
0(e 2 ). Since the computation of Wij for given edges % and j of Q'(-) takes unit time, 
the computation of interference weight functions for all edges of £?'(•) takes 0(e 2 ) time. 
Similarly, the computation of co-schedulability weight functions for all edges of Q'{-) 
requires 0(e 2 ) time. The computation of normalized noise powers at all vertices of Q'(-) 
takes 0(e) time. 



4.4. Performance Results 
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In G'(-), let C denote the total number of colors used to color all vertices and let N 
denote the number of vertices assigned color z, i.e., Ni = \V' C .\. Since C can never exceed 
the number of vertices in G'(-), i.e., the number of edges in Gc(-), C is 0(e). The time 
required by Lines 20-21 is 0(1), let it be ki, where k\ is a constant. 

With a careful implementation of storing J2v[eV' \{v'}u{u'} w 'v' 1 v' ^ v c e V Cp i Lines 
26-28 take 0(1) time. Thus, Lines 25-29 take 0(\V' Cp \) time, let it be equal to k 2 \V c J, 
where k 2 is a constant. Along similar arguments, Lines 30-34 take 0(1) time, let it be 
equal to k 3 , where k 3 is a constant. The time required by Lines 36-38 is k 4 , where k 4 is 
a constant. Thus, the total running time of the coloring phase is 

C Ni 

r = E(^+E(i v -i(^iv;i + /c 3 )+^)). 

i=l j=l 

Since V' UC ,V' C . C V', it follows that \V' UC \, \V' C .\ ^ \V'\ = e. Furthermore, for any color i, 
V; C UV;, C V. Thus, |V4 C ||V;| ^ ^. Therefore 

C C Ni 2 C Ni C Ni 

i=l i=l j=l i=l j'=l i=l J=l 

e 2 

= A;iO + A; 2 — (e) + k 3 e{e) + k A {e) 
k 

= k x C + k 3 e 2 + -je 3 + k 4 e 
= 0(e 3 ). 

Hence, the total running time complexity of SGLS algorithm is 0(e 3 ). ■ 

4.4 Performance Results 

In this section, we demonstrate the efficacy of SGLS algorithm via simulations. To the 
best of our knowledge, for an STDMA network with constrained transmission power, 
there is no existing work on link scheduling that utilizes an SINR graph representation of 
the network. However, for completeness, we compare the performance of SGLS algorithm 
with the CFLS algorithm proposed in Chapter 3. Note that SGLS is based on SINR 
graph while CFLS is based on communication graph and verifying SINR conditions. 

In the simulation experiments, the location of every node is generated randomly using 
a uniform distribution for its X and Y coordinates. We assume that the deployment 
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area is a circular region of radius R. The values chosen for system parameters P, 7 C , 
[3 and N are prototypical values of system parameters in wireless networks [42]. After 
generating random positions for N nodes, we have complete information of $(•). Once 
the link schedule is computed by every algorithm, a is computed using (4.1). For 
a given set of system parameters, we calculate the average spatial reuse by averaging a 
over 1000 randomly generated networks. Keeping all other parameters fixed, we observe 
the effect of increasing the number of nodes N on the average spatial reuse. 

In the first experiment (Experiment 1), we assume that R = 500 m, P = 10 mW, 
(3 = 4, N = —90 dBm and 7 C = 20 dB. Thus, R c = 100 m. We vary the number 
of nodes from 30 to 110 in steps of 5. Figure 4.10 plots the average spatial reuse vs. 
number of nodes for both the algorithms. 

R = 500 m, P = 10 mW, R = 4, N = -90 dBm, y = 20 dB 
r 'c 




60 70 80 

number of nodes 



100 



Figure 4.10: Spatial reuse vs. number of nodes for Experiment 1. 

In the second experiment (Experiment 2), we assume that R = 700 m, P — 15 mW, 
j3 = 4, iVo = —85 dBm and 7 C = 15 dB. Thus, R c = 110.7 m. We vary the number 
of nodes from 70 to 150 in steps of 10. Figure 4.11 plots the average spatial reuse vs. 
number of nodes for both the algorithms. 

From the figures, we observe that SGLS achieves 5-10% higher spatial reuse than 
CFLS. However, this improvement in performance is obtained at the cost of higher 
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R = 700 m, P = 15 mW, R = 4, N = -85 dBm, y = 15 dB 
r 'c 




CD 
CO 



2 




70 80 90 100 110 120 

number of nodes 



130 



140 



150 



Figure 4.11: Spatial reuse vs. number of nodes for Experiment 2. 
computational complexity. 

4.5 Discussion 

In this chapter, we have proposed a novel point to point link scheduling algorithm based 
on an SINR graph representation of an STDMA wireless network under the physical 
interference model. Our results demonstrate that the spatial reuse for the proposed 
algorithm is higher than that of the ConflictFreeLinkSchedule algorithm. This is due to 
the fact that we have embedded interference conditions between pairs of nodes into the 
edge weights and normalized noise powers at receiver nodes into vertex weights of the 
SINR graph and consequently determined a conflict-free schedule. Our approach has 
the potential to scale with the number of nodes in the network. 



Chapter 



5 



Point to 



Multipoint Link 



Scheduling: A Hybrid Approach 



In this chapter, we investigate point to multipoint link scheduling in STDMA wireless 
networks. We generalize the definition of spatial reuse introduced in Chapter 2 for 
point to multipoint link scheduling. We propose a "hybrid" link scheduling algorithm 
based on a communication graph representation of the network and SINR conditions. 
We demonstrate that the proposed algorithm achieves higher spatial reuse than existing 
algorithms, without any increase in running time complexity. 

The rest of this chapter is organized as follows. In Section 5.1, we describe our 
system model. We describe point to multipoint link scheduling based on the protocol 
interference model in Section 5.2 and describe its limitations in Section 5.3. In Section 
5.4, we introduce spatial reuse as our performance metric and formulate the problem. 
In Section 5.5, we describe the proposed link scheduling algorithm. We evaluate its 
performance in Section 5.6 and derive its computational complexity in Section 5.7. We 
discuss the implications of our work in Section 5.8. 



5.1 System Model 



Our system model and notations are exactly as described in Section 2.1. However, we 
redefine and introduce terms that are applicable to point to multipoint link scheduling. 
If node k is within node j's communication range, then k is defined as a neighbor of 



84 



Chapter 5. Point to Multipoint Link Scheduling: A Hybrid Approach 



j, since k can decode fs packet correctly (subject to Equation 2.6). Note that if node k 
is outside node j's communication range, then it can never decode j's packet correctly 
(from Equation 2.6). The number of neighbors of node j is denoted by 

A point to multipoint link schedule for an STDMA wireless network $(•) is a mapping 
from the set of nodes to time slots. Let C denote the number of time slots in a point to 
multipoint link schedule. For a given time slot i, j th point to multipoint transmission is 
denoted by {t it j — > {r^i, • • • , r i,j,r){t i] )}}, where t it j denotes the index of the node 
which transmits a packet and r^i, • • • , r i,j,r)(tij) denote the indices of neighboring 
nodes (neighbors of ijj) that receive the packet. Note that denotes k th receiver of j th 
transmission in time slot %. Let Mj denote the number of concurrent point to multipoint 
transmissions in time slot i. A point to multipoint link schedule for an STDMA network 
$(•) is denoted by ■ ■ • , B c ), where 

Bi -.= {t itl — > {^1,1,^1,2, . . ■ ,n iljV (t iA )}, ■ ■ ■ ,t i>Mi — > {n,M i ,i,n,M i ,2, • • • ^i^^um^}} 

= set of concurrent point to multipoint transmissions in time slot %. 
Every point to multipoint schedule fi(-) must satisfy the following: 

1. Operational constraints: 

(a) A node cannot transmit and receive in the same time slot, i.e., 

{ti,j}n{ri,k,i,---,ri, k , v (ti,k)} = <f> V * = 1,---,C V j ^ k. (5.1) 

(b) A node cannot receive from multiple transmitters in the same time slot, i.e., 

Kj,!, . . . , r idMUj) } n {r iiM , . . . , r iAv{h k) } = <p V % = 1, . . . , C V j ^ k. 

(5.2) 

2. Range constraint: Every receiver is within the communication range of its intended 
transmitter, i.e., 

lht,.; r r,.; h ,) ^ R c V l = 1,...,C V j = 1,...,M, V k = \ //(/, ; ). (5.3) 

For an example, consider the STDMA wireless network $(•) shown in Figure 5.1(a). 
It consists of six nodes whose coordinates (in meters) are 1 = (—40,5), 2 = (0,0), 
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1 = (-40,5) 



5 = (-75,0) 



■2= (0,0) 



3 = (95,0) 



X 



4= (135,0) 



f 6= (0,-75) 

(a) An STDMA wireless network with six nodes. 

point to muitipointpoint link schedule 
point to multipoint transmissions 



4 -{3} 
1 -» {2,5,6} 


6-»{l,2} 


2 -» {1,3,5,6} 


3 -» {2,4} 


5-»{l,2} 


4 -{3} 
1 -> {2,5,6} 




1 


2 


„ 3 


_4 


5 


1 


time 



time slots 

(b) A point to multipoint link schedule for the network shown in Figure 5.1(a). 
Figure 5.1: Example of STDMA network and point to multipoint link schedule. 
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3 = (95,0), 4 = (135,0), 5 = (-75,0) and 6 = (0,-75). One of the possible point 
to multipoint link schedules for this STDMA network is shown in Figure 5.1(b). The 
schedule length is C = 5 time slots and the schedule is defined by Q(Bi, £> 2 , £> 3 , £> 4 , £> 5 ), 
where 

#i = Ki {^1,1,1,^1,1,2, ri,i, 3 }, *i )2 — Ow}} 
= {1^ {2,5,6}, 4 ^{3}}, 

#2 = {*2,1 ^ {^2,l,l,r 2 ,l, 2 }} 

= {6- {1,2}}, 

#3 = {*3,1 — {7-3,1,1^3,1,2^3,1,3^3,1,4}} 

= {2 -.{1,3,5,6}}, 
Ba = (*4,i — {r 4 ,i,i,r 4 ,i, 2 }} 

= {3 -{2,4}}, 
#5 = {*5,i — {r 5 ,i,i,r 5 ,i,2}} 

= {5- {1,2}}. 

After 5 time slots, the schedule repeats periodically, as shown in Figure 5.1(b). 

A point to multipoint link scheduling algorithm is a set of rules that is used to 
determine a schedule fi(-). Typically, a scheduling algorithm is required to satisfy certain 
objectives. 

Consider k th receiver of j th transmission in time slot i, i.e., receiver r^j^. The power 

received at rj & from its intended transmitter U (signal power) is n U p 7. The power 

received at from its unintended transmitters (interference power) is Y^iti WJt~jy- k ) ■ 
Thus, the SINR at receiver rj Jjfc is given by 

p 



SLNR r = m ' • (5.4) 



N + DiTn- 



According to the physical interference model [15], receiver can successfully de- 
code the packet transmitted by Uj if the SINR at r^-^ is no less than the communication 
threshold j c , i.e., 

SINR,,., > 7c . (5.5) 



5.2. Equivalence of Link Scheduling and Graph Vertex Coloring 
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A link schedule Q(-) is exhaustive if every two nodes j, k who are neighbors of each 
other are included in the schedule exactly twice, once with j being a transmitter and k 
being one of its receivers, and during another time slot with k being a transmitter and 
j being one of its receivers. 

5.2 Equivalence of Link Scheduling and Graph Ver- 
tex Coloring 

In this section, we describe the equivalence between a point to multipoint link schedule 
for an STDMA wireless network and the coloring of vertices of the communication graph 
representation (see Section 2.2.1) of the network. 



Parameter 


Symbol 


Value 


transmission power 


P 


10 mW 


path loss exponent 


P 


4 


noise power spectral density 


N 


-90 dBm 


communication threshold 


7c 


20 dB 


interference threshold 


7i 


10 dB 



Table 5.1: System parameters for STDMA networks shown in Figures 5.1(a) and 5.4. 

Consider the STDMA wireless network <£>(•) whose deployment is shown in Figure 
5.1(a). The system parameters for this network are given in Table 5.1. From (2.4), we 
obtain R c = 100 m. The corresponding communication graph representation Q C (V, S c ) is 
shown in Figure 5.2. The communication graph comprises of 6 vertices and 14 directed 
communication edges. The vertex and communication edge sets are given by 

{v 1 ,v 2 ,v 3 ,v 4 ,v 5 ,v 6 }, (5.6) 
{vi A v 2 , v 2 -^v 1 ,v 1 -^ v 5 , v 5 -^v 1 ,v 1 -^ v 6 , v 6 -^v 1 ,v 2 -^ v 5 , 
v 5 A v 2 , v 2 A v 6 , v e A v 2 , v 2 A v 3 , v 3 A v 2 , v 3 A u 4j v 4 A v 3 }. (5.7) 

Given the above representation of an STDMA network, a point to multipoint link 
schedule Q(-) can be considered as equivalent to assigning a unique color to every vertex 
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Figure 5.3: Vertex coloring of communication graph shown in Figure 5.2 corresponding 
to the link schedule shown in Figure 5.1(b). 
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in the communication graph, such that nodes with the same color transmit simultane- 
ously in a particular time slot. For the example network considered, the link schedule 
shown in Figure 5.1(b) corresponds to the coloring of the vertices of the communication 
graph shown in Figure 5.3. Time slots 1, 2, 3, 4 and 5 in fi(-) correspond to colors red, 
blue, green, magenta and yellow in V respectively. Note that a coloring algorithm that 
uses the least number of colors also minimizes the schedule length. 

Algorithms for assigning nodes to time slots (equivalently, colors) require that two 
vertices Vi, Vj can be colored the same if and only if: 

1. edge Vi A vj ^ £ c and edge Vj A Vi ^ £ c , i.e., there is no primary vertex conflict, 
and 

2. there is no vertex such that t) ; ^ 6 £ c and Vj A t> fe e £ c , i.e., there is no 
secondary vertex conflict. 

These criteria are based on the operational constraints (5.1) and (5.2). 

Algorithms based on the protocol interference model represent the network by a com- 
munication graph and utilize various graph coloring methodologies to devise heuristics 
which yield a minimum length schedule. Hence, such algorithms have the merit of low 
computational complexity. However, recent research suggests that these algorithms yield 
low network throughput. This aspect is elaborated in the following section. 

5.3 Limitations of Algorithms based on Protocol In- 
terference Model 

In this section, we illustrate that algorithms based on the protocol interference model 
can result in schedules that yield low network throughput. Note that the limitations of 
point to multipoint link scheduling algorithm are similar to those of point to point link 
scheduling algorithms described in Section 2.3. 

With the intent of maximizing the throughput of an STDMA network, algorithms 
based on the protocol interference model transform the scheduling problem to a ver- 
tex coloring problem for the communication graph representation of the network. For 
example, the BroadcastSchedule algorithm [16] works in two phases. In Phase 1, the 
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vertices of the communication graph are labeled using the labeler function (Algorithm 
2, Section 3.1). In Phase 2, vertices are considered in increasing order of label. For the 
vertex u under consideration, it discards any color that leads to primary or secondary 
vertex conflicts with u. The least color among the residual set of non-conflicting colors 
is used to color vertex u. If no non-conflicting color exists, vertex u is colored with a 
new color. 

The simplification of the link scheduling problem in a wireless network as a vertex 
coloring problem on the communication graph can result in schedules that violate the 
SINR threshold condition (5.5). Specifically, algorithms based on the protocol interfer- 
ence model do not necessarily maximize the throughput of an STDMA network because: 

1. They can result in high cumulative interference at a receiver, due to hard-thresholding 
based on communication radius. This is because the SINR at receiver de- 
creases with an increase in the number of concurrent transmissions Mj, while R c 
has been defined for a single transmission only. 

Y 



1 = (0,0) 



5 = (200,0) 



6 =(370,0) X 



2 = (-80,0) 



3 = (90,0) 



4 = (280,0) 



Figure 5.4: An STDMA wireless network with six nodes. 



Figure 5.5: Communication graph model of STDMA network described by Figure 5.4 
and Table 5.1. 



For example, consider the STDMA wireless network whose deployment is shown 
in Figure 5.4. The network consists of six nodes whose coordinates (in meters) 
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Figure 5.7: Point to multipoint link scheduling algorithms based on protocol interference 
model can lead to high interference. 



92 



Chapter 5. Point to Multipoint Link Scheduling: A Hybrid Approach 



are 1 = (0,0), 2 = (-80,0), 3 = (90,0), 4 = (280,0), 5 = (200,0) and 6 = 
(370, 0). The system parameters are shown in Table 5.1, which yields R c = 100 m. 
The communication graph model of the STDMA network is shown in Figure 5.5. 
Consider the transmission requests 1 — > {2,3} and 4 — > {5,6}, which correspond 
to vertices v\ and V4 of the graph shown in Figure 5.5. Note that vertices v\ and V4 
do not have primary or secondary vertex conflicts. So, to minimize the number of 
colors, such an algorithm will color these vertices with the same color, as shown in 
Figure 5.6. Equivalently, transmissions 1 — > {2, 3} and 4 — > {5, 6} will be scheduled 
in the same time slot, say time slot i. However, our computations show that the 
SINRs at receivers r^i, r^i^, r ij2 ,i and 7^2,2 are 21.85 dB, 12.42 dB, 15.27 dB and 
19.97 dB respectively. Figure 5.7 shows the nodes of the network along with the 
labeled transmitter-receivers sets, receiver-centric communication zones and SINRs 
at the receivers. From the SINR threshold condition (5.5), transmission t^i — > 
is successful, while transmissions t itl — > r^i^, ti,2 —> r i,2,i and t it2 — > ^2,2 are 
unsuccessful. This leads to low network throughput. 

2. Moreover, these algorithms are not aware of the topology of the network, i.e., they 
determine a link schedule without being cognizant of the exact positions of the 
transmitters and receivers. 

As argued above, point to multipoint link scheduling algorithms based on the pro- 
tocol interference model can result in low network throughput. In essence, algorithms 
that construct an approximate model of an STDMA network (communication graph) 
and concentrate on minimizing the schedule length do not necessarily maximize network 
throughput. This observation is developed into a proposal for an appropriate perfor- 
mance metric in Section 5.4. 

5.4 Problem Formulation 

In this section, we motivate the need for a performance metric that takes into account 
the SINR threshold condition (5.5) as the criterion for successful packet reception. Anal- 
ogous to the notion of spatial reuse, we propose a performance metric for point to multi- 
point link scheduling, which is also termed as spatial reuse. We argue that spatial reuse 
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is directly proportional to the number of successful point to multipoint transmissions. 
Finally, we formulate the scheduling problem from a perspective of maximizing spatial 
reuse. 

Algorithms based on the protocol interference model are inadequate to design efficient 
point to multipoint link schedules. This is because these algorithms are entirely based 
on the communication graph Q C (V, £ c ), which is a crude approximation of $(•), and can 
lead to low network throughput, as argued in Section 5.3. On the other hand, from $(•) 
and Gc{-)i one can exhaustively determine the link schedule Q(-) which yields highest 
network throughput according to the physical interference model. However, this is a 
combinatorial optimization problem of prohibitive complexity (O ( | V| ' v ' ) ) and is thus 
computationally infeasible. 

To overcome these problems, we propose a point to multipoint link scheduling algo- 
rithm for STDMA wireless networks under the physical interference model. Our algo- 
rithm is based on the communication graph model Q C (V,£ C ) as well as SINR computa- 
tions. 

To evaluate the performance of our algorithm and compare it with existing link 
scheduling algorithms, we define the notion of spatial reuse. Consider the point to mul- 
tipoint link schedule Q(-) for the STDMA network $(•). Under the physical interference 
model, transmission tij — * is successful if and only if (5.5) is satisfied. The spatial 
reuse of the link schedule Q(-) is defined as the average number of successful point to 
multipoint transmissions per time slot. Thus 

Spatial Reuse = q := ^ i=1 ^ i=1 ^ , (5.8) 

where 1(A) denote the indicator function for event A, i.e., 1(A) = 1 if event A occurs, 
1(A) = if event A does not occur. Note that in (5.8), the number of nodes that 
successfully receive a transmitted packet is normalized by the number of neighbors of the 
transmitting node. A high value of spatial reuse corresponds to high network throughput. 

The essence of STDMA is to have a reasonably large number of simultaneous and 
successful transmissions. For an STDMA wireless network which is operational for a 
long period of time, say L time slots, the total number of successful point to multipoint 
transmissions is Lq. Thus, a high value of spatial reuse directly translates to higher 
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network throughput and the number of colors C is relatively unimportant. Hence, 
spatial reuse turns out to be a crucial metric for the comparison of various STDMA link 
scheduling algorithms. 

Our goal is to design a low complexity point to multipoint link scheduling algorithm 
that achieves high spatial reuse, where spatial reuse is given by (5.8). We only consider 
link schedules that are feasible and exhaustive. 

5.5 MaxAverageSINRSchedule Algorithm 

Our proposed point to multipoint link scheduling scheduling algorithm under the phys- 
ical interference model is MaxAverageSINRSchedule (MASS), which considers the com- 
munication graph Q C (V,S C ) and is described in Algorithm 7. 

Algorithm 7 MaxAverageSINRSchedule (MASS) 
1: input: STDMA wireless network $(•), communication graph Q C (V,S C ) 

2: output: A coloring C : V — > {1, 2, . . .} 

3: label the vertices of Q c randomly {Phase 1} 

4: for j <— 1 to n do {Phase 2 begins} 

5: let u be such that L(u) = j 

6: C(u) ^- MaxAverageSINRColor(u) 

7: end for{Phase 2 ends} 

In Phase 1, we label all the vertices randomly 1 . Specifically, if Q c (-) has v vertices, 
we perform a random permutation of the sequence (1,2, ... ,v) and assign these labels 
to vertices with indices 1,2, ... ,v respectively. Let L(u) denote the label assigned to 
vertex u. 

In Phase 2, the vertices are examined in increasing order of label 2 and the MaxAver- 
ageSINRColor (MASC) function is used to assign a color to the vertex under considera- 
tion. The MASC function is explained in Algorithm 8. It begins by discarding all colors 
that have a primary or secondary vertex conflict with u, the vertex under consideration. 

Randomized algorithms are known to outperform deterministic algorithms, especially when the 

characteristics of the input are not known a priori [53] . 

2 In essence, the vertices are scanned in a random order, since labeling is random. 



5.6. Performance Results 



95 



Algorithm 8 integer MaxAverageSINRColor(w) 
1: input: STDMA wireless network $(•), communication graph £ C (V, £ c ) 

2: output: A non-conflicting color 

3: C <— set of existing colors 

4: C p <— {C(x) : x is colored and is a neighbor of u} 

5: C s <— {C(x) : x is colored and is two hops away from u} 

6: C nc = C\ {C p U C s } 

7: if C nc 7^ then 

8: r <— color in C„ c which results in maximum average SINR at neighbors of -u 

9: if maximum average SINR ^ 7 C then 
10: return r 
11: end if 
12: end if 
13: return \C\ + 1 



Among the set of non-conflicting colors C nc , it chooses that color for u which results in 
the maximum value of average SINR at the neighbors of u, provided this value exceeds 
the communication threshold. Intuitively, the average SINR is also a measure of the 
average distance of every neighbor of u from all co-colored transmitters. The higher the 
average SINR, the higher is this average distance. We choose that color which results 
in the maximum average SINR at the neighbors of u, so that the additional interference 
at the neighbors of all co-colored transmitters is kept low. If no such color is found, it 
assigns a new color to u. 

5.6 Performance Results 

In this section, we describe our simulation model. We compare the performance of the 
proposed algorithm with existing point to multipoint link scheduling algorithms. 

In our simulation experiments, the location of every node is generated randomly 
in a circular region of radius R. If (Xj,Yj) are the Cartesian coordinates of node j, 
then Xj ~ U[-R,R] and Yj ~ U[-R,R] subject to X] + Yf < R 2 . Equivalently, if 
(Rj,@j) are the polar coordinates of node j, then R? ~ U[0, R 2 } and Qj ~ ?7[0,27r]. 
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Using (2.4) and (2.5), we compute R c and R iy and then map the STDMA network 
$(•) to the communication graph Q(V,£ C ). Once the link schedule f2(-) is computed 
by every algorithm, the spatial reuse q is computed using (5.8). We use two sets of 
prototypical values of system parameters in wireless networks [42]. For a given set of 
system parameters, we calculate the average spatial reuse by averaging q over 1000 
randomly generated networks. Keeping all other parameters fixed, we observe the effect 
of increasing the number of nodes on the average spatial reuse. 

In our experiments, we compare the performance of the following algorithms: 

1. BroadcastSchedule (BS) [16] 

2. MaxAverageSINRSchedule (MASS) 

In our first set of experiments (Experiment 1), we assume that R = 500 m, P = 10 
mW, (3 = 4, N = -90 dBm, 7c = 20 dB and 7i = 10 dB. Thus, R c = 100 m and 
Ri = 177.8 m. We vary the number of nodes from 30 to 110 in steps of 5. Figure 5.8 
plots the average spatial reuse vs. number of nodes for both the algorithms. 
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Figure 5.8: Average spatial reuse vs. number of nodes for Experiment 1. 



In our second set of experiments (Experiment 2), we assume that R = 700 m, P — 15 
mW, (3 = 4, N = -85 dBm, 7c = 15 dB and 7i = 7 dB. Thus, R c = 110.7 m and 
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Ri = 175.4 m. We vary the number of nodes from 70 to 150 in steps of 5. Figure 5.9 
plots the average spatial reuse vs. number of nodes for both the algorithms. 
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Figure 5.9: Average spatial reuse vs. number of nodes for Experiment 2. 



From Figures 5.8 and 5.9, we observe that average spatial reuse increases with the 
number of nodes for both the algorithms. The MASS algorithm consistently yields 
higher average spatial reuse compared to BS. The spatial reuse of MASS is about 15% 
higher than BS in Experiment 1 and 4% higher in Experiment 2. This improvement in 
performance translates to substantially higher network throughput. 

Also, an increase in the number of nodes in a given geographical area leads to an 
increase in the number of vertices having a primary or secondary vertex conflict with 
a given vertex. Hence, the number of non-conflicting colors for a given vertex also de- 
creases. From this reduced set of non-conflicting colors, BroadcastSchedule chooses a 
color randomly, while MaxAverageSINRSchedule chooses a color based on SINR condi- 
tions. Since spatial reuse takes SINR threshold conditions into account, the gap between 
average spatial reuse values increases with number of nodes in Figures 5.8 and 5.9. 
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5.7 Analytical Result 

In this section, we derive an upper bound on the running time (computational) com- 
plexity of the MaxAverageSINRSchedule algorithm. Let v denote the number of vertices 
of the communication graph Q C {V,S C ). 

Theorem 5.7.1. The running time of MaxAverageSINRSchedule is 0(v 2 ). 

Proof. Assuming that an element can be chosen randomly and uniformly from a finite set 
in unit time ([53], Chapter 1), the running time of Phase 1 can be shown to be 0(v). In 
Phase 2, the vertex under consideration is assigned a color using MaxAverageSINRColor. 
The worst-case size of the set of colors to be examined \C nc U C p U C s \ is 0(v). With a 
careful implementation, MaxAverageSINRColor runs in time proportional to \C nc \, i.e., 
0(v). Thus, the running time of Phase 2 is 0(v 2 ). Finally, the overall running time of 
MaxAverageSINRSchedule is 0(v 2 ). ■ 

5.8 Discussion 

In this chapter, we have developed a point to multipoint link scheduling algorithm for 
STDMA wireless networks under the physical interference model, namely MaxAvera- 
geSINRSchedule. The performance of our algorithm is superior to existing algorithms. 
A practical experimental modeling shows that, on an average, our algorithm achieves 
15% higher spatial reuse than the BroadcastSchedule algorithm [16]. Since link sched- 
ules are constructed offline only once and then used by the network for a long period of 
time, this improvement in performance directly translates to higher network throughput. 
The computational complexity of MaxAverageSINRSchedule is also comparable to the 
computational complexity of BroadcastSchedule. Therefore, MaxAverageSINRSchedule 
is a good candidate for efficient STDMA point to multipoint link scheduling algorithms. 



Chapter 6 



A Review of Random Access 
Algorithms for Wireless Networks 

The MAC problem or multipoint to point problem is present in all communication 
networks, both wired and wireless. Multiple nodes (users) can access a single channel 
simultaneously to communicate with each other or a common receiver - the challenge is 
to design efficient channel access algorithms to achieve the desired performance in terms 
of throughput and delay. Several solutions to the MAC problem have been proposed 
depending on source traffic characteristics, channel models and Quality of Service (QoS) 
requirements of the users. 

MAC protocols can be broadly classified into two types: fixed resource allocation 
protocols and random access protocols. Fixed resource allocation protocols such as Time 
Division Multiple Access (TDMA), Frequency Division Multiple Access (FDMA) and 
Code Division Multiple Access (CDMA) assign orthogonal or near-orthogonal channels 
to every user and are mostly implemented in voice-dominant wireless cellular networks. 
These protocols typically require the presence of a central entity (base station) to perform 
channel allocation and admission control, i.e., they are highly centralized. Though fixed 
resource allocation protocols are contention-free and can multiplex users with similar 
traffic characteristics easily, they suffer from low throughput and high channel access 
delay when the traffic is bursty and there are large number of users. On the other hand, in 
random access protocols, users vary their transmission probabilities or transmission times 
based on limited channel feedback, i.e., random access protocols are highly distributed. 
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Random access protocols are more suitable for scenarios wherein many users with varied 
traffic requirements have to be multiplexed, i.e., the traffic is bursty. 

Random access algorithms for satellite communications, multidrop telephone lines 
and multitap bus ("traditional random access algorithms") have been well studied for 
the past four decades. These algorithms can be broadly classified into three categories: 
ALOHA [17], [59], Carrier Sense Multiple Access [60] and tree (or stack or splitting) 
algorithms [18]. Traditional random access algorithms have been implemented in practi- 
cal systems. For example, ALOHA is used in most cellular networks to request channel 
access and also in satellite communication networks. Carrier Sense Multiple Access with 
Collision Detection (CSMA/CD) is used to resolve contentions in Local Area Networks 
(LANs). 

On the other hand, random access algorithms that incorporate physical layer char- 
acteristics such as SINR and channel variations have only been studied recently. These 
algorithms, which have been primarily proposed for wireless networks, can be broadly 
classified into three categories: algorithms based on signal processing and diversity tech- 
niques, channel-aware ALOHA algorithms based on adapting the retransmission prob- 
abilities of contending users and "tree-like" algorithms based on adapting the set of 
contending users. Existing random access algorithms, such as Carrier Sense Multiple 
Access with Collision Avoidance (CSMA/CA ), are not channel-aware and can lead to 
low throughput. Thus, the design of physical layer aware random access algorithms can 
be a potential step towards achieving higher data rates in future wireless networks. 

The organization of this chapter is as follows. Section 6.1 provides a summary of tra- 
ditional random access algorithms along with the canonical system model, performance 
metrics and well-known random access techniques such as ALOHA and tree algorithms. 
This helps us understand channel-aware generalizations of these algorithms. In Section 
6.2, we review research papers which employ signal processing and diversity techniques 
to correctly decode packets in random access wireless networks. We critically review 
some of the research which focus on channel-aware ALOHA and tree-like algorithms 
for wireless networks in Sections 6.3 and 6.4 respectively. Finally, we motivate the use 
of variable transmission power to increase the throughput of random access wireless 
networks in Section 6.5. 
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6.1 Traditional Random Access Algorithms 

In this section, we describe the idealized slotted multiaccess model, which can be used to 
represent various multiaccess media such as satellite channels, multidrop telephone lines 
and multitap bus. We explain traditional random access algorithms such as ALOHA 
and tree algorithms. We also describe the performance metrics used to analyze and 
evaluate random access algorithms, namely, throughput, delay and stability. 

Consider an idealized slotted multiaccess system with m transmitting nodes and one 
receiver. The assumptions of the model are [22]: 

1. Slotted system: All transmitted packets have the same length and each packet 
requires one time unit, called a slot, for transmission. 

2. One of the following is usually assumed: 

(a) Poisson arrivals: Packets arrive at each of the m nodes according to an inde- 
pendent Poisson process. Let A be the overall arrival rate to the system and 
let — be the arrival rate at each transmitting node. 

(b) Backlogged model: Every node always has a packet to transmit. Once a 
node transmits a packet successfully, a new packet is generated and awaits 
transmission. 

3. Collision or perfect reception: If two or more nodes transmit a packet in a given 
slot, then there is a collision and the receiver obtains no information about the 
contents or the sources of transmitted packets. If only one node transmits a packet 
in a given slot, the packet is correctly received. 

4. {0, 1, e} immediate feedback: At the end of each slot, every node obtains feedback 
from the receiver specifying whether packet, 1 packet or more than one packet 
(e denotes error) were transmitted in that slot. 

5. Retransmission of collisions: Each packet involved in a collision must be retrans- 
mitted in some later slot, with further such retransmissions until the packet is 
successfully received. A node with a packet that must be retransmitted is said to 
be backlogged. 
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6. Only one of the following is assumed: 

(a) No buffering: If one packet at a node is currently waiting for transmission or 
colliding with another packet during transmission, new arrivals at that node 
are discarded and never transmitted. 

(b) Infinite set of nodes: The system has an infinite set of nodes and each new 
packet arrives at a new node. 

For the analysis and performance evaluation of random access algorithms, the metrics 
of interest are: 

1. Delay: Index packets as 1, 2, 3, . . . according to their arrival instants. Let Dj denote 
the delay experienced by j th packet. Then the average packet delay is defined as 



V = lim E 



m ^ J 

3=1 



(6.1) 



2. Throughput: The following are the two most common definitions of throughput: 

(a) Throughput is the supremum of input packet arrival rates A such that the 
packet delay remains bounded, i.e., 

7i = sup A. (6.2) 

D<oo 

(b) Let n(t) denote the number of packets successfully transmitted in [0, t\. Define 



T(A) = 



Throughput is then defined as 




if V < oo, 
otherwise. 



T 2 = supT(A). (6.3) 

A 

3. Stability: A random access algorithm is stable if the throughput T > and 
unstable if T = 0. 

The research of random access algorithms began with the unslotted ALOHA (pure 
ALOHA) algorithm proposed by Abramson [17]. Each node, upon receiving a packet, 
transmits it immediately rather than waiting for a slot boundary. If a packet is involved 
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in a collision, it is retransmitted after a random delay. It can be shown that unslotted 
ALOHA achieves a maximum throughput of ^ ^ 0.1839 [22]. An advantage of unslotted 
ALOHA is that it can be used with variable-length packets. 

Slotted ALOHA is a variation by Roberts [59] of the original unslotted ALOHA pro- 
tocol proposed by Abramson. Each node simply transmits a newly arriving packet in 
the first slot after the packet arrival. When a collision occurs, every node sending a 
colliding packet discovers the collision at the end of the slot and becomes backlogged. 
Backlogged nodes wait for a random number of slots before retransmitting. The maxi- 
mum throughput of slotted ALOHA can be shown to be ± « 0.3678 [22]. Drift-analytic 1 
methods reveal that slotted ALOHA is unstable. To stabilize ALOHA, some techniques 
estimate n or p r , so as to maintain the attempt rate G(n) at 1, resulting in a maximum 
stable throughput of \ [61], [62]. Unlike unslotted ALOHA, slotted ALOHA cannot be 
easily used with variable-sized packets. In slotted ALOHA, long packets must be broken 
up to fit into slots and short packets must be padded out to fill up slots. 

Keeping the random access spirit of the ALOHA protocol, researchers attempted to 
design more efficient protocols. A highly successful approach consists of improving the 
control of the channel by carrier sensing, i.e., the Carrier Sense Multiple Access (CSMA) 
technique. In [22], the authors show that CSMA outperforms ALOHA. Research has 
shown that CSMA based protocols can achieve a throughput close to 0.9 [63]. The 
Ethernet protocol, which is used to connect computers on a wired LAN, utilizes Carrier 
Sense Multiple Access with Collision Detection (CSMA/CD). 

In splitting algorithms, the set of colliding nodes splits into subsets, one of which 
transmits in the next slot. For a given colliding node, the choice of its subset depends 
on a pre-determined rule such as, the outcome of tossing an unbiased coin, a function 
of its arrival time or a function of its node identifier. If the collision is not resolved, a 
further splitting into subsets takes place. The algorithm proceeds recursively until all 
collisions are resolved. 

In the Basic Tree Algorithm (BTA) [18], when a collision occurs, say in k th slot, all 
nodes not involved in the collision go into a waiting mode, and all those involved in the 
collision split into two subsets, according to the pre-determined rule. The first subset 

1 Drift in state n is defined as the expected change in backlog over one time-slot, starting in state n. 
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Figure 6.1: Basic Tree Algorithm for three nodes a, b and c. 
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Table 6.1: Transmitting and waiting sets for basic tree algorithm shown in Figure 6.1. 
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Figure 6.2: Stack representation of transmitting and waiting nodes for basic tree algo- 
rithm shown in Figure 6.1. 
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("left" subset) transmits in slot k + 1, and if that slot is idle or successful, the second 
subset ("right" subset) transmits in slot k + 2. Alternatively, if another collision occurs 
in slot k + 1, the first of the two subsets splits again, and the second subset waits for the 
resolution of that collision. Figure 6.1 exemplifies the operation of BTA for three nodes 
a, b and c. Observe the binary tree structure of the sets of transmitting and waiting 
nodes in the figure. The transmitting and waiting sets in terms of subtrees of this binary 
tree are shown in Table 6.1, where U = {a, b, c} denotes the set of all nodes that were 
involved in the initial collision. The labeling of the subtrees is recursive; for example, LR 
denotes the right subtree of the left subtree of the original binary tree. The transmission 
order corresponds to that of a stack, as shown in Figure 6.2. In each slot, the stack is 
popped and all the nodes that were at the top of the stack transmit their packets. In 
case of a collision, the stack is pushed with nodes that join the right subset and then 
pushed again with nodes that join the left subset. In case of a success or idle, no push 
operations are performed on the stack. A Collision Resolution Period (CRP) is defined 
to be completed when a success or idle occurs and there are no remaining elements on 
the stack. In Figure 6.1, the length of the CRP is 9 slots. 

During the operation of BTA, many new packets might arrive while a collision is 
being resolved. To solve this problem, at the end of a CRP, the set of nodes with new 
arrivals is immediately split into j subsets, where j is chosen so that the expected number 
of packets per subset is slightly greater than 1. The maximum throughput, optimized 
over the choice of j as a function of expected number of waiting packets, is 0.43 packets 
per slot [18]. 
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Figure 6.3: Modified Tree Algorithm for three nodes a, b and c. 
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Table 6.2: Transmitting and waiting sets for modified tree algorithm shown in Figure 



abc 



abc 



be 



be 



be 



10 



time 



Figure 6.4: Stack representation of transmitting and waiting nodes for modified tree 
algorithm shown in Figure 6.3. 
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There exist variants of BTA that yield higher throughput. For example, in Modified 
Tree Algorithm (MTA), if a collision in slot k is followed by an idle in slot k + 1, then 
nodes which collided in slot k refrain from transmitting in slot k + 2. Instead, they 
further split into two subsets, one of which transmits in slot k + 2. As an example, the 
operation of MTA for three nodes a, b and c is shown in Figure 6.3. Observe that the 
length of the CRP is 8 slots. For this example, the transmitting and waiting sets of 
subtrees are shown in Table 6.2, while the corresponding stack representation is shown 
in Figure 6.4. If an idle occurs in the current slot and a collision occurred in the previous 
slot (see Slot 5 in Figure 6.4), then the stack is popped a second time but the nodes at 
the top of the stack are not transmitted. Instead, these nodes split into right and left 
subsets and these subsets are pushed on the stack. This leads to fewer collisions and 
higher throughput compared to BTA. The maximum stable throughput of MTA is 0.46 
packets per slot [64]. 

In First Come First Serve (FCFS) splitting algorithm [22], nodes involved in a col- 
lision split into two subsets based on the arrival times of collided packets. Using this 
approach, each subset consists of all packets that arrived in some given interval, and 
when a collision occurs, that interval will be split into two smaller intervals. By always 
transmitting packets that arrived in the earlier interval first, the algorithm transmits suc- 
cessful packets in the order of their arrival. The FCFS algorithm is stable for A < 0.4871 
[22] . Conflict resolution protocols based on tree algorithms have provable stability prop- 
erties [65]. 

We should point out that the random access algorithm proposed in Chapter 7 has a 
"tree structure" analogous to that of MTA. The detailed explanations of BTA and MTA 
provide a basic background to understand the dynamics of the proposed algorithm. 

So far, we have summarized the methodology of traditional random access algo- 
rithms. In subsequent sections, we will focus on random access algorithms that are 
better suited for wireless networks such as WLANs and Wireless Metropolitan Area 
Networks (WMANs). 
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6.2 Signal Processing in Random Access 

The collision model (Section 6.1, Assumption 4) is simple in that the onus of scheduling 
packets is left entirely to the MAC layer. On the contrary, physical layer techniques 
like multipacket reception, capture and network-assisted diversity are able to correctly 
decode packets from collisions by means such as coding and signal processing. These 
techniques are potential steps towards alleviating the burden of decoding packets from 
the MAC layer to the physical layer [66]. In this section, we review representative re- 
search work which exploits signal processing and diversity techniques to correctly decode 
the received packets in random access wireless networks. 

With the advent of multiaccess techniques such as CDMA and Multiuser Detection 
[67], the first fundamental change in the collision model has been propounded in [68]. The 
authors offer the generalization that, in the presence of simultaneous transmissions, the 
reception can be described by conditional probabilities instead of deterministic failure. 
They propose the MultiPacket Reception (MPR) model defined by the matrix 



C = 



Cio Cn 
C20 C21 C22 



(6.4) 



where CV,- is the conditional probability that, given i users transmit, j out of i transmis- 
sions are successful. Given k users transmit at the same time, the average number of 
successfully received packets is given by 

k 

C k = ^2jC kj . (6.5) 
3=0 

They show that ALOHA under MPR achieves stable throughput lim^oo Ck assuming 
that the limit exists. The stability and delay of finite-user slotted ALOHA with multi- 
packet reception has been analyzed in [69]. 

In [70], the authors analyze the probability of capture in a multipoint to point wire- 
less network. Analogous to the physical interference model, the capture model assumes 
that if a user's SINR exceeds a threshold 7, then that user's packet will be successfully 
received. They consider a realistic multiplicative propagation model in which the re- 
ceived power is obtained by multiplying the transmitted power by independent random 
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variables representing fading, shadowing and path loss effects. To model the near-far 
effect, they assume that the distance r of a mobile station from the base station is a 
random variable with distribution function Fn(r). They show that, under broad condi- 
tions, the roll-off parameter 5 of the distribution of power received from a mobile station 
is determined by the path loss exponent and Fr(t). Additionally, 5 is insensitive to 
other effects such as Rayleigh or Rician fading and log-normal shadowing. Finally, they 
show that in the limit of a large number of transmitters, the probability of capture is 
determined by the power capture threshold 7 and 5. Though the analysis provided in 
[70] is mathematically robust, the authors do not describe any multiple access algorithm 
which achieves high throughput in wireless networks under the capture model, i.e., their 
result is more existential than constructive. 

In [71], the authors propose Network- Assisted Diversity Multiple Access (NDMA), 
a technique for resolving collisions in wireless networks. They consider a wireless slot- 
ted random access network with Rayleigh fading and Additive White Gaussian Noise 
(AWGN). In NDMA, if k users collide in a given slot, they repeat their transmissions 
k — 1 times so that k copies of the collided packets are received. Using signal separation 
principles, the receiver resolves a k x k source mixing problem to extract the signals of 
individual users, without incurring any penalty in throughput. The protocol has been 
extended to blind user detection [72] and has provable stability [73]. A good review of 
NDMA protocols is given in [74]. 

An alternative to employing signal processing techniques in random access wireless 
networks is to appropriately model the wireless channel and modify the well-researched 
ALOHA protocol. We review such research work in the next section. 

6.3 Channel-Aware ALOHA Algorithms 

In this section, we review representative research work whose central theme is to adapt 
the retransmission probabilities of users in random access wireless networks. In other 
words, we review research work which develops channel-aware ALOHA algorithms for 
wireless networks. 

In [75], the authors develop a channel-aware ALOHA protocol for wireless networks. 
They assume a slotted system, block fading, {0, l,e} feedback and a backlogged model 
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(Section 6.1, Assumption 2b). They develop a distributed random access protocol in 
which each node only has knowledge of its own channel gain and nodes have long- 
term power constraints. A node transmits only if its channel gain exceeds Hq. For 
a system with n nodes, the authors show that the optimum transmission probability 
is where a(n) G (0, 1] and a(n) — > 1 as n — > oo. Asymptotically, the ratio of 
the throughput of channel-aware ALOHA to the throughput of a centralized scheduler 
(which has knowledge of channel gains of all nodes) is shown to be \. 

Opportunistic ALOHA algorithms for wireless networks have been studied in [76]. 
The authors consider a general reception model which encompasses {0,1, e} feedback, 
capture as well as multipacket reception. Under the assumption that the Channel State 
Information (CSI) is known to each user, they propose a variant of slotted ALOHA, 
where the transmission probability is allowed to be a function of the CSI. The maximum 
throughput for the finite-user infinite-buffer model is derived. Finally, the theory is 
applied to CDMA networks with Linear Minimum Mean Square Error (LMMSE) receiver 
and matched filters. 

The performance of slotted ALOHA in a wireless network with multiple destinations 
under the physical interference model is evaluated in [19]. A packet is successful only if 
it is captured at the receiver of its intended destination. The authors assume Poisson 
packet generation, {0,1, e} feedback and circularly symmetric Gaussian distribution of 
users around each destination. They use a modified version of Rivest's pseudo-Bayesian 
estimator [22] to estimate the backlog. Their simulation results demonstrate the effect of 
arrival rate, capture threshold, variance of user distribution and number of destinations 
on the throughput and energy efficiency per destination. 

In [77], the author analyzes the throughput of slotted ALOHA in a multipoint to 
point wireless ad hoc network under the physical interference model. The cluster head 
employs reverse link power control, similar to IS-95 CDMA systems [78], to ensure that 
equal power is received from all nodes who attempt transmission in a time slot. The 
wireless channel is modeled as a multipacket reception channel. Assuming that one new 
packet arrives at each node in every time slot, the state of the system is characterized 
by a discrete time Markov chain with a steady state distribution. Finally, the author 
describes a technique to compute the network throughput. 

In [79] , the authors introduce spatial reuse slotted ALOHA, a random access proto- 
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col for random homogeneous mobile wireless networks. The occurrence of a collision is 
determined by the SINR at a receiver, i.e., the authors employ the physical interference 
model. They assume that nodes are randomly placed in a two-dimensional plane accord- 
ing to a Poisson point process and each node chooses a random destination at some finite 
distance. The powers at which stations can transmit are assumed to be independent and 
identically distributed (i.i.d.) and the wireless channel is characterized by its propaga- 
tion path loss. Nodes move according to the random waypoint mobility model [80]. The 
authors characterize the interference process using tools from stochastic geometry. Sub- 
sequently, they determine the probability of channel access that maximizes the expected 
projected distance traversed per hop towards a destination, termed as "spatial density 
of progress". Under the assumption that there is some non-degenerate node mobility, 
the authors show that the spatial density of progress is proportional to the square-root 
of the density of the nodes. Though the authors present a distributed ALOHA protocol 
and address certain implementation issues, their model does not represent real-world 
scenarios. Practical deployments of wireless nodes are better modeled by a uniform 
distribution in a finite plane rather than a Poisson point process in an infinite plane. 
Also, most of their results do not hold for static wireless networks (say, wireless mesh 
networks) since ergodicity assumptions no longer hold. Finally, their proposed routing 
protocol requires every node to have knowledge of locations and MAC states (receiver or 
transmitter) of all other nodes, which requires a lot of message passing between nodes 
(especially with mobile nodes) and is thus not scalable. 

Instead of adapting the transmission probabilities of users in random access wireless 
networks, one can also adapt the transmission times of users based on the channel 
state and feedback from the receiver. Such techniques, which can be broadly termed 
as splitting algorithms or tree-like algorithms for wireless networks, are reviewed in the 
next section. 

6.4 Splitting Algorithms 

In this section, we review representative research work on random access algorithms 
whose main idea is to adapt the set of contending users based on feedback from the 
channel or the common receiver. In such work, the authors develop and analyze splitting 
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(or tree or stack) algorithms for various models of the wireless channel and evaluate the 
performance of their algorithms via simulations. 

In [81], the authors propose an opportunistic splitting algorithm for a multipoint to 
point wireless network. They assume a slotted system, block fading channel and {0, 1, e} 
feedback. Assuming that each user only knows its own channel gain and the number of 
backlogged users, the authors propose a distributed splitting algorithm to determine the 
user with the best channel gain over a sequence of mini-slots. The algorithm determines 
a lower threshold Hi and a higher threshold H^ for each mini-slot, such that only users 
whose channel gains lie between between Hi and Hh are allowed to transmit their packets. 
Based on results from "partitioning a sample with binary type questions" [82], they 
show that the average number of mini-slots required to determine the user with the best 
channel is 2.5, independent of the number of users and the fading distribution. However, 
their algorithm is impractical because it assumes that every user can accurately estimate 
the number of backlogged users. 

In [83], the authors consider a random access network with infinite users, Poisson 
arrivals and {0, k, e} immediate feedback, where k is any positive integer. In contrast 
to standard tree algorithms (BTA, MTA, FCFS) that discard collided packets (Section 
6.1, Assumption 4), they propose an algorithm that stores collided packets. The receiver 
extracts information from the collided packets by relying on successive interference can- 
cellation techniques ([67], Chapter 7) and the tree structure of a collision resolution 
algorithm. Though their algorithm achieves a stable throughput of 0.693, it requires 
infinite storage and increased input voltage range at the receiver, which are not feasible 
in practical systems. 

In [84] , the author considers a multipoint to point wireless channel with and without 
capture and MPR. The channel provides Empty(E)/Non-Empty(NE) feedback to all 
active users and 'success' feedback to successful users only. The users do not need to 
know the starting times and ending times of collision resolution periods. For such a 
channel with E/NE binary feedback, the author proposes and analyzes a stack multiple 
access algorithm that is limited sensing and does not require any frame synchronization. 
The author considers two models for capture, namely Rayleigh fading with incoherent 
and coherent combining of joint interference power. For MPR, the author assumes a 
maximum of two successes during a collision. The maximum throughput of the algorithm 
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is numerically evaluated be to 0.6548 when capture and MPR are present, and 0.2891 
when both effects are absent. Though a novel splitting algorithm is proposed in [84], 
the author does not take into account throughput gains possible by varying transmission 
powers of users. 

So far, we have reviewed research papers that either utilize signal processing tech- 
niques or adapt transmission probabilities or transmission times to increase the through- 
put in random access wireless networks. The throughput can be further increased by 
allowing users to use variable transmission powers. We review research papers which 
employ this idea in the next section. 

6.5 Towards Power Controlled Random Access 

In this section, we review representative research papers which focus on power control 

techniques in random access wireless networks. We then motivate the use of variable 

transmission power to increase the throughput in random access wireless networks. 

In [85] , the author considers a time-slotted CDMA-based wireless network wherein a 

finite number of nodes communicate with a common receiver. The author formulates the 

problem of determining the set of nodes that can transmit in each slot along with their 

corresponding transmission powers, subject to constraints on maximum transmission 

power and the SINRs of all transmissions exceeding the communication threshold. Due 

to its NP-hard nature, the problem is relaxed to a case wherein a node transmits with 

a certain probability in each slot. Equivalently, the problem of joint power control 

and link scheduling is transformed to a problem of power controlled random access, 

wherein the objective is to determine the probability of transmission Aj and transmission 

power Pi for each node i, subject to constraints on maximum transmission power and 

the "expected SINR" exceeding the communication threshold. The author seeks to 

minimize a weighted sum of the maximum transmission power and maximum reciprocal 

probability, i.e., minimize (max« Pj + Amaxj ^-). This convex optimization problem is 

solved using techniques from geometric programming [86]. Finally, the author derives 

the probability of outage 2 and delay distribution of buffered packets and demonstrates 

2 An outage occurs on a link if the received (actual) SINR on the link is less than the communication 
threshold. 
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the efficacy of the schemes via simulations. 

In [87], the authors investigate transmission power control and rate adaptation in 
random access wireless networks using game theoretic techniques. They consider mul- 
tiple transmitters sharing a time-slotted channel to communicate equal-length packets 
with a common receiver. A user's packet is successfully received if the SINR at the re- 
ceiver is no less than the communication threshold, i.e., the authors employ the physical 
interference model. The random access problem is formulated as a game wherein each 
user selects its strategy (transmit or wait) at each stage of the game in a non-cooperative 
(independent) or cooperative manner. The authors evaluate equilibrium strategies for 
non-cooperative and cooperative symmetric random access games. Finally, the authors 
describe distributed power control and rate adaptation games for non-cooperative users 
for a collision channel with power-based capture. Their numerical results demonstrate 
improved expected user utilities when power control and rate adaptation are incorpo- 
rated, at the expense of increased computational complexity. Though the authors pro- 
pose a distributed random access algorithm based on game theoretic techniques, their 
algorithm is impractical because it assumes that every user knows n, the number of 
backlogged users, in each slot. However, in practice, n can only be estimated using 
techniques such as Rivest's pseudo-Bayesian algorithm [62]. 

Though researchers have addressed the problem of random access in wireless networks 
by considering various channel models, different types of feedback and realistic criteria for 
successful packet reception, only few of them exploit the idea that throughput gains are 
achievable in a random access wireless network by varying transmission powers of users. 
In general, varying the transmission powers of users leads to higher long-term average 
power. However, there exist wireless networks whose users do not have stringent energy 
requirements. For such scenarios, it would be useful to investigate the throughput gains 
achievable in the network by varying the transmission powers of users. 

We envisage developing a power controlled random access algorithm for wireless net- 
works under the physical interference model. We seek an algorithm that yields higher 
throughput than traditional random access algorithms. In cognizance of these require- 
ments, we propose a power controlled splitting algorithm for wireless networks in Chapter 
7. The algorithm is so designed that successful packets are transmitted in the order of 
their arrivals, i.e., in an FCFS manner. 
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In the system model considered in Chapter 7, if multiple transmissions occur, the 
receiver can decode a certain user's packet correctly only if the received SINR exceeds a 
threshold, i.e., we consider a channel with power-based capture. The notion of capture 
has been addressed previously, though in different contexts [19], [84], [88]. However, in 
Chapter 7, we motivate the idea that a user can transmit at variable power levels to 
increase the chances of capture. Moreover, unlike [19], [84], [88], we assume {0,1, c, e} 
feedback, where 0, 1 and e denote idle, success and error respectively (Section 6.1, 
Assumption 4), and c denotes capture in the presence of multiple transmissions. Note 
that the system model considered in Chapter 7 is different from those considered in 
existing works on splitting algorithms for wireless networks. For example, in [84], the 
author proposes a novel splitting algorithm, but does not take into account throughput 
gains possible by varying the transmission power. Though the authors of [81] propose 
a splitting algorithm to determine the user with the best channel gain, their algorithm 
is impractical because it assumes that every user can accurately estimate the number of 
backlogged users. 

To the best of our knowledge, there is no existing work on variable power splitting al- 
gorithms for a wireless network under the physical interference model. The specification 
of the proposed algorithm along with its performance analysis and evaluation constitute 
the subject matter of the next chapter. 



Chapter 7 

Power Controlled FCFS Splitting 
Algorithm for Wireless Networks 

In this chapter, we propound a random access algorithm that incorporates variable trans- 
mission powers in a multipoint to point wireless network. Specifically, we investigate 
random access in wireless networks under the physical interference model wherein the 
receiver is capable of power-based capture, i.e., a packet can be decoded correctly in 
the presence of multiple transmissions if the received SINR exceeds the communication 
threshold. We propose an interval splitting algorithm that varies the transmission pow- 
ers of users based on channel feedback. We derive the maximum stable throughput of 
the proposed algorithm and demonstrate that it achieves better performance than the 
FCFS splitting algorithm [22] with uniform transmission power. 

The rest of the chapter is organized as follows. We describe our system model in 
Section 7.1 and motivate variable control of transmission powers of contending users 
in Section 7.2. We describe the proposed random access algorithm and provide two 
illustrative examples in Section 7.3. We model the algorithm dynamics by a Markov 
chain and derive its maximum stable throughput in Section 7.4. The performance of the 
proposed algorithm is evaluated in Section 7.5. We conclude in Section 7.6. 

7.1 System Model 

Consider a multipoint to point wireless network. We assume the following: 



118 



Chapter 7. Power Controlled FCFS Splitting Algorithm for Wireless Networks 



1. Slotted system: Users (nodes) transmit fixed- length packets to a common receiver 
over a time-slotted channel. All users are synchronized such that the reception of 
a packet starts at an integer time and ends before the next integer time. 

2. Poisson arrivals: The packet arrival process is Poisson distributed with overall rate 
A, and each packet arrives to a new user that has never been assigned a packet 
before. After a user successfully transmits its packet, that user ceases to exist and 
does not contend for channel access in future time slots. 

3. Channel model: The wireless channel is modeled by the path loss propagation 
model. The received signal power at a distance D from the transmitter is given by 

where P is the transmission power and is the path loss factor. We do not 
consider fading and shadowing effects. 

4. Power-based capture: According to the physical interference model [15], a packet 
transmission from transmitter tij to receiver r in i th time slot is successful if and 
only if the SINR at receiver r is greater than or equal to the communication 
threshold 7c 1 , i.e., 



where 

Mi 

P 

No 

5. {0, 1, c, e} immediate feedback: By the end of each slot, users are informed of the 
feedback from the receiver immediately and without any error. The feedback is 
one of: 

1 In literature, 7 C is also referred to as capture ratio [84], capture threshold [19] and power ratio 
threshold [70]. 



p. 

gSg > 7 ( 7 n 



number of concurrent transmitters in i th time slot, 
j th transmitter in i th time slot (j — 1,2, ... , Mi), 
Euclidean distance between Uj and r, 
transmission power of Uj, 
thermal noise power spectral density. 
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(a) idle (0): when no packet transmission occurs, 

(b) perfect reception (1): when one packet transmission occurs and is received 
successfully, 

(c) capture (c): when multiple packet transmissions occur and only one packet 
is received successfully, or 

(d) collision (e): when multiple packet transmissions occur and no packet recep- 
tion is successful. 

The receiver can distinguish between 1 and c by using energy detectors [83], [89]. 
Thus, by the end of every slot, only two bits are required to provide feedback from 
the receiver to all users. Note that two bits are required to provide feedback even 
for the classical {0, 1, e} feedback model. Thus, our {0, 1, c, e} immediate feedback 
assumption does not increase the number of bits required for feedback.. 

6. Gated Channel Access Algorithm (CAA): New packets are transmitted in the first 
available slot after previous conflicts are resolved. The time interval from the slot 
where an initial collision occurs up to and including the slot in which all users 
recognize that all packets involved in the collision have been successfully received, 
is called a Collision Resolution Period (CRP). Thus, new arrivals are inhibited 
from transmission during the CRP. 

7. Equal distances: We assume that each user is at the same distance D from the 
common receiver. 

7.2 Motivation and Problem Formulation 

The maximum stable throughput of the well-known FCFS splitting algorithm is 0.4871 
[22] , which is the highest throughput amongst a wide class of random access algorithms 
for wired networks. However, in a wireless network, transmission power of a node pro- 
vides an extra degree of freedom, and higher throughputs are achievable. 

Consider a scenario wherein all contending nodes transmit with equal power P in a 
given time slot. When only one node transmits, its packet is successfully received if the 
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SINR threshold condition (7.1) is satisfied, i.e., 

P > j c N D?. (7.2) 

When M nodes transmit concurrently with equal power P, where M ^ 2, the SINR 
corresponding to i th transmission is given by 

p 

SINR* = -5-, (7.3) 

a quantity which is always less than 1. Since j c > 1 for all practical narrowband 
communication receivers [70], SINRj < j c V % and all M transmissions are unsuccessful 2 . 
Thus, when multiple nodes transmit with equal power, a collision occurs irrespective of 
the transmission power P. 

However, the above situation can be circumvented by varying transmission powers of 
users in some special cases. With relatively small attempt rates, when a collision occurs, 
it is most likely between only two packets [22]. In this case, if the receiver is capable 
of power-based capture, a collision between two nodes can be avoided by using different 
transmission powers. Specifically, one of the nodes, say Aq, transmits with minimum 
power Pi such that, if it were the only node transmitting in that time slot, then its 
packet transmission will be successful. From (7.1), the required nominal power is 

P 1 = lc N Q D?. (7.4) 

The other node, say A^, transmits with minimum power P2 such that if there is exactly 
one other node transmitting at nominal power Pi, then the packet transmitted by N 2 
will be successful. From (7.1) and (7.4), we obtain 

AT -L p i 1 ' 

P 2 = lc (N D? + Pi), 
:.Pi = 70(1 + 70)^0^- (7.5) 

Note that ^ = 1 + 7 C . We do not consider more than two power levels for the following 
reasons: 



2 For a spread spectrum CDMA system with processing gain L, (7.3) gets modified to SINR; = 

p 

jv +1(m-i) p t^]- For such a wideband system, j c < 1, and more than one packet can be decoded 
correctly in the presence of multiple transmissions. However, in this thesis, we consider narrowband 
systems only. 
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1. it complicates the power-control algorithm, and 

2. most mobile wireless devices have constraints on peak transmission power. 

Note that the above power control technique converts some collisions into "captures". 
Thus, it has the potential of increasing the throughput of random access algorithms 
employing uniform transmission power. 

We seek to design a distributed algorithm incorporating this power control technique, 
while still ensuring that the algorithm transmits successful packets in the order of their 
arrival, i.e., in an FCFS manner 3 . 

7.3 PCFCFS Interval Splitting Algorithm 

In this section, we present an algorithmic description of the proposed Power Controlled 
First Come First Serve (PCFCFS) splitting algorithm. We also explain the behavior of 
the proposed algorithm by providing two illustrative examples. 

7.3.1 Description 

We first describe the notation. Slot k is defined to be the time interval [k, k+1). At each 

integer time k (k ^ 1), the algorithm specifies the packets to be transmitted in slot k to 

be the set of all packets that arrived in an earlier interval [T(k), T(k) + 4>(k)), which is 

defined as the allocation interval for slot k. The maximum size of the allocation interval is 

denoted by <p , a parameter which will be optimized for maximum throughput in Section 

7.4. Packets are indexed as 1, 2, ... in the order of their arrival. Since the arrival times 

are Poisson distributed with rate A, the inter-arrival times are exponentially distributed 

with mean j-. Let <Zj denote the arrival time of i th packet. Using the memoryless 

property of the exponential distribution (and without loss of generality), we assume 

that ai = 0. The transmission power of i th packet in slot k is denoted by Pi(k), where 

Pi(k) G {0, Pi, P 2 }. Note that, if Pj(fc) = 0, then i th packet is not transmitted in slot k. 

3 Since successful packets are transmitted in an FCFS manner, the delay experienced by a packet 
will not be significantly higher than the average packet delay. Thus, from a QoS perspective, FCFS 
transmission of packets not only guarantees average delay bounds, but also ensures fairness of user 
packets. 
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Algorithm 9 describes the proposed Power Controlled First Come First Serve (PCFCFS) 
splitting algorithm, which is the set of rules by which the users compute allocation in- 
terval parameters {T(k + l),0(fc + l),cr(k + 1)} and transmission power Pi(k + 1) for 
slot k + 1 in terms of the feedback and allocation interval parameters for slot k. In 
our algorithm, every allocation interval is tagged as a "left" (£) or "right" (1Z) inter- 
val. a(k) denotes the tag (C or TV) of allocation interval [T(k),T(k) + <ft(k)) in slot k. 
Moreover, whenever an allocation interval is split, it is always split into two equal-sized 
subintervals, and these subintervals (£,7Z) are said to correspond to each other. 

In Phase 1 of the algorithm, we initialize various quantities, r denotes the number 
of slots for which the algorithm operates; ideally r — > oo. By convention, the initial 
allocation interval is [0, mm(</» , 1)), which is a right interval (71). The initial channel 
feedback is assumed to be idle (0). 

In Phase 2 of the algorithm, we determine power levels, obtain channel feedback and 
compute allocation interval parameters for each successive slot k. In Phase 2a, all users 
whose arrival times lie in the left half of the current allocation interval transmit with 
higher power P 2 , while all users whose arrival times lie in the right half of the current 
allocation interval transmit with nominal power P\. However, if a capture occurred 
in the previous slot k — 1, all users in the current allocation interval transmit with 
nominal power P ± . Therefore, our algorithm always transmits successful packets in an 
FCFS manner. In Phase 2b, the allocation interval parameters are modulated based 
on the channel feedback. More specifically, if a collision occurs, then the left half of 
the current allocation interval becomes the new allocation interval. If a capture occurs, 
then the right half of the current allocation interval becomes the new allocation interval. 
If a success occurs and the current allocation interval is tagged as a left interval, then 
the corresponding right interval becomes the new allocation interval. If an idle occurs 
and the current allocation interval is tagged as a left interval, then the left half of the 
corresponding right interval becomes the new allocation interval. Otherwise, if a success 
or an idle occurs and the current allocation interval is tagged as a right interval, the 
waiting interval truncated to length O becomes the new allocation interval, and a new 
Collision Resolution Period (CRP) begins in the next time slot k + 1. Note that the 
transmit power levels in PCFCFS are variable and based on channel feedback, i.e., they 
are adaptive. 
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input: O , Pi, P 2 , arrivals a lt a 2 , a 3 , . . . in [0, r) {Phase 1 begins} 

T(l) <- 

0(1) <- min(0 o , 1) 

cr(i) = ?e 

feedback = {Phase 1 ends} 
for k <— 1 to t do {Phase 2 begins} 
if feedback ^ c then {Phase 2a begins} 
for all i such that T(k)^ ai <T(k)+^ do 

^ (A;) = P 2 
end for 

for all % such that T(fc)+^<ai<T(k)+4>(fc) do 

Pi(k) = P l 
end for 
end if {Phase 2a ends} 

transmit packets whose arrivals times lie in [T(k),T(k) + (p(k)) and obtain channel 
feedback {Phase 2b begins} 
if feedback = e then 
T(k + 1) <— T(k) 

a(k + l) <- C 
else if feedback = c then 
+ 1) <- T(fc) + ^ 



a(ife + l) <- ft 
else if feedback = 1 and a(k) = £ then 

r(fc + 1) <- T(fc) + 4>{k) 
4>(k + i) <- 0(fc) 

a(ife + l) <- ft 
else if feedback = and cr(fc) = £ then 

r(fc + 1) <— + 4>{k) 

a(k + l) <- £ 
else 

+ 1) <- + 0(fc) 
0(fc + l) = min(0 o , k - T(k)) 
a(k + l) <- ft 
end if {Phase 2b ends} 
end for{Phase 2 ends} 
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7.3.2 Examples 
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Figure 7.1: PCFCFS splitting algorithm illustrating a collision followed by an idle. 



To illustrate the rules of the PCFCFS splitting algorithm for a single CRP, consider 
the example shown in Figure 7.1. In slot k, the allocation interval has no node in its 
left half and two nodes in its right half. Both these nodes transmit with nominal power 
Pi and a collision occurs. So, the allocation interval is split, with the left interval L 
being the allocation interval for slot k + 1. An idle occurs in slot k + 1. Next, the right 
subinterval R is further split, with RL being the allocation interval for slot k + 2. The 
left node in RL transmits with higher power P 2 , while the right node in RL transmits 
with nominal power P ± , resulting in a capture of the packet transmitted by the left node. 
The allocation interval is further split, with RLR forming the allocation interval for slot 
k + 3. Since a capture occurred in RL in slot k + 2, the corresponding right subinterval 
RR is returned to the waiting interval in slot A; + 3. Post-capture, the lone node in RLR 
transmits with nominal power Pi, resulting in a success and completing the CRP. For 
the same sequence of arrival times, the behavior of the FCFS algorithm with uniform 
transmission power P is shown in Figure 7.2, where a(k) denotes the length of the 
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Figure 7.2: FCFS splitting algorithm illustrating a collision followed by an idle. 



126 



Chapter 7. Power Controlled FCFS Splitting Algorithm for Wireless Networks 



allocation interval in slot fc. Note that the FCFS algorithm requires 5 slots to resolve 
the collisions, while the proposed PCFCFS algorithm requires only 4 slots. 
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Figure 7.3: PCFCFS splitting algorithm illustrating a collision followed by another 
collision. 



To further illustrate the rules of the PCFCFS splitting algorithm for a single CRP, 
consider the example shown in Figure 7.3. In slot k, the allocation interval has three 
nodes in its left half and one node in its right half. All 'left half nodes transmit with 
higher power P2, while the 'right half node transmits with nominal power Pi, leading to 
a collision. So, the allocation interval is split, with the left interval L being the allocation 
interval for slot fc + 1. In slot fc + 1, the allocation interval has one node in its left half, 
which transmits with higher power P 2 , and two nodes in its right half, which transmit 
with nominal power Pi. Hence, a collision occurs, and the allocation interval L is split 
into two equal sized subintervals LL and LR, with LL being the allocation interval for 
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Figure 7.4: FCFS splitting algorithm illustrating a collision followed by another collision. 
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slot k+2. Since a collision is followed by another collision, the right interval R is returned 
to the waiting interval in slot k + 2. In slot k + 2, there is only one node in the allocation 
interval. Since this lone node lies in the right half of the allocation interval, it transmits 
with nominal power Pi, leading to a success. Thus, LR becomes the allocation interval 
for slot k + 3. For this allocation interval, the node in the left half transmits with higher 
power F 2 and the node in the right half transmits with nominal power Pi, resulting in 
a capture of the packet transmitted by the former node. Consequently, LRR becomes 
the new allocation interval for slot k + 4. Finally, in slot k + 4, the lone node transmits 
with nominal power Pi, leading to a deterministic success and completing the CRP. For 
the same sequence of arrival times, the behavior of the FCFS algorithm with uniform 
transmission power P is shown in Figure 7.4. Note that the FCFS algorithm requires 
6 slots to resolve the collisions, while the proposed PCFCFS algorithm requires only 5 
slots. 



7.4 Throughput Analysis 

The evolution of a CRP can be represented by the Discrete Time Markov Chain (DTMC) 
shown in Figure 7.5. Every state in the DTMC is a pair (cr,i), where a is the status 
{L, V , R, R' , C} and i is the number of times the original allocation interval (of length O ) 
has been split. State (R, 0) corresponds to the initial slot of a CRP. If an idle or a success 
occurs, the CRP ends immediately and a new CRP begins in the next slot. If a capture 
occurs, a transition occurs to state (C, 1), where C indicates that capture has occurred 
in the allocation interval. If a collision occurs in (R, 0), a transition occurs to state 
(L, 1). Each subsequent idle in a left allocation interval generates one additional split 
with a smaller left allocation interval, corresponding to a transition to (L',i + 1), where 
V indicates that the current left allocation interval has been reached after a collision 
(in some time slot) followed by one or more idles. A collision in an allocation interval 
generates one additional split with a smaller left allocation interval, corresponding to a 
transition to (L,i + 1), where L indicates that the current left allocation interval has 
been reached just after a collision. A capture in an allocation interval generates an 
additional split with a smaller right allocation interval and corresponds to a transition 
to (C,i + 1). This is followed by a success from (C,i + 1) to (R,0), thus ending the 
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Figure 7.5: Discrete Time Markov Chain representing a CRP of PCFCFS splitting 
algorithm. 
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CRP. A success in a left allocation interval leads to the corresponding right allocation 
interval with no additional split, which causes a transition from (L, i) to (R, i), or (L', i) 
to (R', i). A success in (R', i) causes a transition to (R, 0), thus ending the CRP. It can 
be easily verified that the states and transitions in Figure 7.5 constitute a Markov chain, 
i.e., each transition from every state is independent of the path used to reach the given 
state. 

We now analyze a single CRP. Assume that the size of the initial allocation interval 
is 0o (corresponding to state (R, 0)). Each splitting of the allocation interval halves 
this, so that states (L,i), (L',i), (R,i), (R',i) and (C, i) in Figure 7.5 correspond to 
allocation intervals of size 2~ l O - Since the arrival process is Poisson with rate A, the 
number of packets in the original allocation interval is a Poisson random variable (r.v.) 
with mean X(/) . Consequently, the a priori distributions on the number of packets in 
disjoint subintervals are independent and Poisson. Define Gi as the expected number of 
packets in an interval that has been split i times. Thus 

d = 2-^00 = 2-^0 V i > 0, (7.6) 
:.Gi = \g^ V i^l. (7.7) 

We view (R, 0) as the starting state as well as the final state. For brevity in notation, 
the transition probability from state (A,i) to state {B,j) is denoted by where 
A,B e {L, L', R, R', C} and i, j e {0} UZ + (see Figure 7.5). For example, the transition 
probability from (L, 1) to (C, 2) is denoted by Pl 1: c 2 - 

L R 

(x L packets) ( yXR packets) 

I 1 1 1 1 

LL LR RL RR 

(xll packets) {xlr packets) (xrl packets) (%rr packets) 

Figure 7.6: Notation for number of packets in left and right subintervals of the original 
allocation interval. 

Pr ,Ro 1S the probability of an idle or success in the first slot of the CRP. Since 
the number of packets in the initial allocation interval is Poisson with mean Go, the 
probability of or 1 packet is 

Pro.Ro = (l + G )e~ Go . (7.8) 
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-Pr ,Ci is the probability of capture in the first slot of a CRP. Let xl and xr denote 
the number of packets in the left and right halves of the original allocation interval 
respectively, as shown in Figure 7.6. Capture occurs if and only if xl — 1 and xr = 1. 
xl and xr are independent Poisson r.v.s of mean G\ each. Thus 

Pro,Ci = Pr(x L = l,x R = 1), 

= Pr(x L = l)Vi(x R = 1), 

= G^ 2G \ 

•••^oA = fe^'. (7.9) 

State (L, 1) is entered after collision in state (R, 0). Using (7.8) and (7.9), this occurs 
with probability 



PRoM — 1 - Pr ,R - Pr ,Ci, 

.■.Prom = 1- (l + G + ^)e- Go . (7.10) 
Since a capture is always followed by a deterministic success, 



1 V % > 1. (7.11) 



Lemma 7.4.1. The outgoing transition probabilities from (L,i), where i ^ 1, are given 
by 

(1 - e~ G > - Gie- Gi )Gie- Gi 
p Li,Ri - t ^r^; — - — , (7-12) 



Lp i 

4 C 



1- (l + Gi-i + ^le-^-i 



1 - (1 + ^ + ^)e" G ' 

-Gi-i 



1- (l + Gi_i + %)e- 



Proof. Refer to Figure 7.5. For z = 1, (L, i) is entered only via a collision in (i?, z — 1). 
For i = 2, (L, i) is entered only via a collision in (L, i — 1) or (R,i — 1). For i ^ 3, (L, z) is 
entered only via a collision in (L', z — 1), (L, z — 1), (i?, z — 1) or (i?', z — 1). In every case, 
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a subinterval Y is split into YL and YR, and YL becomes the new allocation interval. 
Let xyl and xyr denote the number of packets in YL and YR respectively. A priori, 
xyl and xyr are independent Poisson r.v.s of mean Gi each. The event that a collision 
occurred in the previous state is {xyl + xyr ^ 2} fl {xyl = %yr = 1} C =: Cy. Note that 
xyl + xyr = xy is a Poisson r.v. of mean Gi-\. From (7.7), Gi = \Gi-\ V i ^ 1. The 
probability of success in (L,i) is the probability that xyl = 1 conditional on Cy, i.e., 

P LijRi = Pl(xy L = l\Cy), 

Pt{C y \x yl = l)Pr(xy L = 1) 
Pr(Cy) 

Pr({x YR = 1} n {x YR = 1} C ) Pr(x YL = 1) 
Pr(Cy) 
Pr(xy fi > 2) Pr(xy L = 1) 

Pr({xy L + x yi? > 2} n {xy L = £y R = 1} C ) ' 

Pr(x YR > 2) Pt(x yl = 1) 
Pr(xy ^ 2) - Pr(a;y L = 1) Pr(x YR = 1) ' 

(1 _ e -Gi _ G ie - Gi )Gie- Gi 
1 _ g-G*-! _ G^ie-^-i - G?e- 2G * ' 
(1 - e" G < - Gie- Gi )Gie- Gi 
1- (l + Gi_i + %i)e-Gi-i 

The probability of idle in (L,i) is the probability that xyl = conditional on Cy, 

i.e., 

P LiA+1 = Pt(x yl = 0\C y ), 

Pr(Cy|xy L = 0) Pr(xy L = 0) 
Pr (Cy) 

Pr({x YR > 2} n {x YR = 1} C ) Pr(x yL = 0) 

Pr(Cy) 

Pr(x YR > 2) Pr(xy L = 0) 

Pr({x YL + x Y r > 2} n {x YL = x YR = 1} C ) ' 

Pr(xy R ^ 2) Pr(x YL = 0) 

Pr(xy ^ 2) - Pr(xy L = 1) Pr(x YR = 1) ' 

(l-e~ G ' - G,e- G ')e- G ' 

1 _ e - G -i - Gi_ ie - G '-i - Gfe- 2G * ' 

(1 _ e-Gi _ G ie - Gi )e- Gi 
• ■ Li,L' i+1 - 7 ^r^; • (/-L'J 

Let xyll and £y_L_R denote the number of packets in YLL and YLR respectively. 
xyll and xylr are independent Poisson r.v.s of mean G i+ i each, and xyll+xylr = ^yl- 



7.4. Throughput Analysis 



133 



The probability of capture in (L,i) is the probability that xyll = 1 and xylr = 1 
conditional on Cy, i.e., 

PLi,c i+1 = Pt(xyll = 1,xylr = 1\C y ), 

Pt(C y \x Y ll = I, xylr = 1) Pt(x Y ll = l,x YLR = 1) 

Pt({x YL + X Y R > 2} n {X YL = X YR = 1} C ) 

Pr(C y \x YL = 2) Pt(xyll = 1) ^(xylr = 1) 

Pr({xy L + X Y R > 2} n {X YL = X Y R = 1} C ) ' 

Pr(x rfl ^ 0) Pr(x yLL = 1) Pt(x YL r = 1) 
Pr(x Y >2)- Pr(x YL = 1) Pr(x YR = 1) ' 
l.G? +1 e^ 

1 - e- G -i - Gi-ie- *-! - G\e-^ ' 

.:P Li ,c i+1 = 7 4 " s v • (7.18) 



l-(l + Gi_i + %)e-<* 
From (7.16), (7.17) and (7.18), we obtain 



l 



Phi,L i+l 1 Phi,Ri Phi,L' i+1 PLi,Ci+ii 

1- (i + G i + f) e -G< 
^,L i+1 = H> • (7.19) 



1- + + 



e - G »-i 



Lemma 7.4.2. The outgoing transition probabilities from (R,i) are given by 

= 1 - (f+ fr)r^ V ^ (7 ' 20) 
1 _ d + (7. + ^ e ~ Gl 

p <«- = i-(i + G<) ^. v 1 > L (7 ' 21) 



Proof. Refer to Figure 7.5. For i ^ 1, (i?, z) is entered only via a success in (L, z). Recall 
that (L, z) was entered only via a collision from a previous state. We use the notation 
introduced in the proof of Lemma 7.4.1. Define the event 

S YL ■= C Y H {x YL = 1}, 

= {x Y l + x Y r > 2} n {x YL = x YR = 1} C H {x YL = 1}, 

= {x Y R>i}n{x Y R = i} c n{xYL = i}, 

:.S YL = {x Y r >2}n{x YL = 1}. (7.22) 
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Let xyrl and xyrr denote the number of packets in YRL and YRR respectively. 
xyrl and xyrr are independent Poisson r.v.s of mean G i+i each. Since xyr ^ 2, a 
success or an idle can never occur in state (R, i). Note that xyr = xyrl + xyrr- The 
probability of capture in state (R, i) is the probability that xyrl = 1 and xyrr = 1 
conditional on Syl, i-e., 

Pr i: C 1+1 = Pr(x Y RL = l,XYRR = l\xYR>2,XYL = l), 

= Pt(x Y rl = 1,x Y rr = 1\x Y r ^ 2), 

Pr(xy R ^ 2|x yflL = 1, x Y rr = 1) Pr(o:y RL = 1, x YRR = 1) 

Pr(x YR > 2) 

Pr(^yfiL + x Y rr > 2|x yRL = l,x yflfi = 1) Pr(s yfl£ , = l,x yflR = 1) 

Pr(x YR > 2) 
1.Pt(x Y rl = 1)Pt(xyrr = 1) 
Pr(x YR > 2) 
G 2 -e -20 ^ 1 



1 - e~ G * - Git~ Gi ' 



g?„ -a 



• p ^a +1 = r^iTG^- (7 - 23) 

From (7.23), we obtain 



p _ i p 



1 



G' 2 r 

4 e 



l_(l + G .) e -G- 

G?- 



1- (l + Gi + ^)e- Gi 



Lemma 7.4.3. The outgoing transition probabilities from (L', i) are given by 

P t ,, L , +i = (1 ,- e "-- G -VT G - V,>2, (7.26) 

= 1 - (1 fCe^ V ' * 2 ' <7 ' 27> 

2 



*w - u^,:^.- ^ 2 - « 7 - 28 > 
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Proof. Refer to Figure 7.5. For i = 2, (£/, z) is entered only by an idle in (L, i — 1). For 
z ^ 3, state (L', z) is entered by an idle in (Z/, z — 1) or an idle in (L, i — 1). In every case, 
a residual right subinterval, say Z, is split into ZL and Zi?, and ZL becomes the new 
allocation interval. Note that (L', i) can be entered if and only if there is a collision (in 
some time slot) followed by one or more idles. Therefore, Z must contain at least two 
packets. Let xzl and xzr denote the number of packets in ZL and ZR respectively. A 
priori, xzl and x Z r are independent Poisson r.v.s of mean Gi each. Let xz denote the 
number of packets in Z. Thus xz = x Z l + %zr, %z is a Poisson r.v. of mean and 
x z >2. 

The probability of success in (L',i) is the probability that xzl = 1 conditional on 
x z > 2, i.e., 

P L[A = Pr{x ZL = l\xz>2), 

Pr(x z > 2\x ZL = l)Pr{x ZL = 1) 

Pr(x z > 2) 
Pt(x Z l + x Z r > 2\x Z l = 1) P?(x Z l = 1) 

Pr(x z > 2) 
Pr(x ZR > 1) Pr(x ZL = 1) 
Pr(x z > 2) 

The probability of idle in (Z/, z) is the probability that x Z l = conditional on x z ^ 2, 

i.e., 

P L5)L5+1 = Pr(x ZL = 0\x z > 2), 

PrQr z >2|q; ZL = 0)Pr(x ZL = 0) 

Pr(x z ^ 2) 
PrQr ZL + x Z R > 2\x ZL = 0) Pr(x ZL = 0) 

Pr(x z > 2) 
PrQr Zi ^2)Pr(x ZL = 0) 
Pr(x z ^ 2) 

(1 - e" Gl - G i e~ G '0e~ G ' 1 ,„ N 

•■■^«« = \-(i + G J,)U-, • < 7 - 30 ' 

Let xzll and xzlk denote the number of packets in ZLL and ZZi? respectively. A 
priori, xzll and xzlr are independent Poisson r.v.s of mean G i+ \ each. The probability 
of capture in (£/, z) is the probability that xzll = 1 and rc^Li? = 1 conditional on xz ^ 2, 
i.e., 
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Pl' v c i+1 = Pt(x Z ll = 1,x Z lr = l\x z > 2), 

Pr(x z > 2\x Z ll = 1, xzlr = 1) Pt(x Z ll = 1, xzlr = 1) 

Pr(x z > 2) 
l.Pr(x ZLL = l)Pr(x ZL R = l) 
Pr(x z > 2) 

Q2 i6 -2G i+1 



1 - e^'- 1 - G , j_ie~ Gi - 1 ' 



G 2 r 

•'• /Vr = i-a+G^)^-' (7 ' 31) 

From (7.29), (7.30) and (7.31), we obtain 

Pv.,L l+1 = 1 - Pl>.,r>. - Pl>,l>- Pl>.,c>, (7-32) 



1 - (1 + G t + 



e 



Ph '^ ~ 1 - (1 + G^e-^ ■ (7 ' 33) 



Lemma 7.4.4. TTie outgoing transition probabilities from (R 1 , i) are given by 

Pr>,r = ^zg7 V i> 2, (7.34) 
= j^za- V^2, (7.35) 



l-(l + G, + %)e- 
1 - e" 



= V ; V < > 2. (7.36) 



Proof. Refer to Figure 7.5. For i ^ 2, state (i?', i) is entered if and only if a success 
occurs in state (L', i). When (L f , i) was entered, a residual right subinterval Z was split 
into ZL and Zf?, and ZL became the new allocation interval. Recall that x z ^ 2, since 
(I/, i) can only be entered after a collision followed by one or more idles. A success in 
(L 1 , i) implies xzl = 1- Hence, (R', i) is entered if and only if both these events occurs, 
i.e., x z ^ 2 and x ZL = 1. Therefore, (i?', i) can be entered if and only if x Z r ^ 1. Note 
that there can never be an idle from (it!', i). 

The probability of success in (R', i) is the probability that xzr = 1 conditional on 
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x Z r ^ 1, i.e., 

Pr^Ro = ^(x ZR = 1\x Z r > 1), 

Pt(x zr >1\x zr =1)Pt(x zr =1) 
Pr(x Zi? ^ 1) 

■■•^o = Y^al- (7-37) 

Let xzrl and denote the number of packets in ZRL and ZRR respectively. 

Note that xzr = xzrl + %zrr- xzrl and xzrr are independent Poisson r.v.s of mean 
G i+ i each. The probability of capture in state (it!', i) is the probability that xzrl = 1 
and x Z rr = 1 conditional on Xzr ^ 1, i.e., 

PR' v C l+1 = P?(XZRL = ^-,X Z RR = 1\X Z R > 1), 

Pr(x Z R ^ 1|x Zj rl = I, x Z rr = 1)Pt(x Z rl = l,x ZRR = 1) 

Pr(x ZR > 1) 
l.PrQr^ L = l)Pr(x z ^ = l) 
Pt(x zr > 1) 



1 — e _Gi 



From (7.37) and (7.38), we obtain 



Pr^M+i - 1 - PR'i,Ro - PR'i,C i+1 , 

G? 



l-(l + C, + %)e- 
1 - e" 



P ^,L i+1 = V , ' 4 ' • (7-39) 



In summary, Figure 7.5 is a DTMC and the transition probabilities are given by 
(7.8), (7.9), (7.10) and (7.11), and Lemmas 7.4.1, 7.4.2, 7.4.3 and 7.4.4. 

We now analyze the DTMC in Figure 7.5. Observe that no state can be entered 
more than once before the return to (R, 0). Let Qx t denote the probability that state 
(X, i) is entered before returning to (R, 0), where X £ {L, L', R, R', C} and % £ Z + . In 
other words, Q Xi denotes the probability of hitting (X, i) in a CRP given that we start 
from (R, 0). Note that Qci = Pr ,c\ and = Pr ,l!- The probabilities Qx, can be 
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calculated iteratively from the initial state (R, 0) as follows: 

Q<h = ^fe- G °, (7.40) 



2, 



Qin = l-(l + G + ^)e- G °, (7.41) 

Qc 2 = Ql.Pl^c+Qr.PruC (7.42) 

Ql' 2 = QiaPi*^, (7.43) 

Ql 2 = Q Ll P Ll ,L 2 +QRi P RuL 2 , (7.44) 

Qu = Qv^Pv^Li + Ql^Pl^uL, + Qr^Pr^uL, 

+ Qr>_PrU^ V 2>3, (7.45) 

Ql> % = Qlu^U^+^-^WJ V <>3, (7.46) 

Q Rl = QlA,* Vi>l, (7.47) 

= Ql^Pl'^ V i > 2, (7.48) 

Qc, = Qlu^UA+^^A + ^-i^-iA + Q^^A 3.(7.49) 

Let random variable K denote the number of slots in a CRP. Thus, K equals the 
number of states visited in the Markov chain, including the initial state (R, 0), before 
the return to (R, 0). Thus 

oo 

E[K] = 1 + J2(Ql 1 + Ql>+Qr 1 +Qr>+Qc 1 ), (7.50) 
i=\ 

where we assume Ql> = Qr> = 0. 

We evaluate the change in T(k) from one CRP to the next, i.e., we evaluate the 
difference in left endpoints of initial allocation intervals of successive CRPs. For the 
assumed initial interval of size 0o, this change is at most 4>q. However, if left allocation 
intervals have collisions or captures (e.g., RL in Figure 7.1), then the corresponding 
right allocation intervals (e.g., RR in Figure 7.1) are returned to the waiting interval, 
and the change is less than </> . Let random variable F denote the fraction of <po returned 
in this manner over a CRP, so that <f> (l — F) is the change in T(k). We distinguish 
between two cases: 

1. If a left allocation interval of type (L,i) has a collision or a capture, then the 
corresponding right allocation interval (R, i) is returned to the waiting interval. 
Let denote the probability that (L, i) has a collision or a capture. Hence, 
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denotes the probability that (L,i) has two or more packets. Thus, = 

P^m+i + P Li,c i+1 - Using (7.14) and (7.15), we obtain 

1 _ (i + G-)e~ Gi 
U Lt = / { +G f 2 S V^l. (7.51) 

2. If a left allocation interval of type (I/, i) has a collision or a capture, then the 
corresponding right allocation interval (R', i) is returned to the waiting interval. 
Let Uu denote the probability that (L',i) has a collision or a capture. Hence, 
Ul> denotes the probability that (L',i) has two or more packets. Thus, U L '. = 
Pv v L i+ i + p L' v c i+1 - Using (7.27) and (7.28), we obtain 

1-(1 + Gi)e- Gi 

In either case, the fraction of the original allocation interval returned on such a 
collision or a capture is 2~\ Therefore, the expected value of F is 

oo 

E[F] = ^2(Q Li U Li + Q L >U L {)2- i , (7.53) 
i=i 

where we assume JJu =0. 

From (7.6), (7.50) and (7.53), we observe that E[K] and E[F\ are functions only 
of the product A</> - Note that as i — > oo, Gj = 2~ i A0 o — > 0. Using the Taylor series 
expansion for e x or L'Hopital's Rule, we can easily prove that: 



1. 



&«w = \' (7 - 54) 

= \< < 7 - 55 > 

1™Pl;, C)+ , = i, (7.56) 

lim^,L l+1 = \, (7.57) 



limP^Ho = 1, (7.58) 
limP^ A+1 = 0, (7.59) 
limPflj^ = 0, (7.60) 
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3. 



limPz^ = 0, (7.61) 

i — >oo 

KmP LiA+l = i (7.62) 

limPL iA+1 = i (7-63) 

lim^W = j, (7.64) 



4. 



limP^A^ = \, (7.65) 
\imP Ri:Li+1 = \. (7.66) 

The proofs of these results are given in Appendix A. Hence, Ql v Q L '., Q R r and Q Ci tend 
to zero with increasing i as 2~ l , while Q Ri tends to zero with increasing i as 4~\ Thus, 
E[K] and -E^-F] can be easily evaluated numerically as functions of \(f) . 

Define the time backlog to be the difference between the current time and the left 
endpoint of the allocation interval, i.e., k — T(k). Note that all packets that arrived in 
the interval T(k), k have not yet been successfully transmitted, i.e., they are backlogged. 
Moreover, we define the drift D to be the expected change in time backlog, k — T(k), over 
a CRP, assuming an initial allocation interval of O - Thus, D is the expected number of 
slots in a CRP less the expected change in T(k), and is given by 

D = E[K] - 0o (1 - E[F\). (7.67) 

The drift is negative if E[K] < O (1 — E[F\). Equivalently, the drift is negative if 

A < ^yauc. 

The right hand side of (7.68), £, is a function of A0o and is plotted in Figure 7.7. We 
observe that £ takes its maximum value at A0o = 1-4. More precisely, £ has a numerically 
evaluated maximum of 0.5518 at A0 O = 1-4. If 0o is chosen to be ^ 5 4 lg = 2.54, then 
(7.68) is satisfied for all A < 0.5518. Thus, the expected time backlog decreases whenever 
it is initially larger than </>o, and we infer that the algorithm is stable for A < 0.5518. 
We have therefore proved the following result. 

Proposition 7.4.5. The maximum stable throughput of the PCFCFS algorithm is 0.5518. 
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Figure 7.7: Plot of £ versus A0r> 

7.5 Numerical Results 

In our numerical experiments, we use values of system parameters that are commonly 
encountered in wireless networks [42]. We compare the performance of the following 
algorithms: 

1. FCFS with uniform power P 1; 

2. PCFCFS. 

For each algorithm, the value of the initial allocation interval is chosen so as to achieve 
maximum stable throughput. For FCFS, maximum stable throughput occurs when its 
initial allocation interval, a = 2.6 [22]. From Section 7.4, the maximum throughput of 
PCFCFS occurs at 4> = 2.54. Let n suc denote the number of successful packets in [0, r) 
and di denote the departure time of i th packet. 
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For a given set of system parameters, we compute the following performance metrics: 



Throughput = 
Average Delay = 

Average Power = 



££r(*-ai) 



n s 



En S uc sr^ii ~p OA 
i=l 2^k=\a,i] r i\^) 



(7.69) 
(7.70) 

(7.71) 



Keeping all other parameters fixed, we observe the effect of increasing the arrival rate 
on the throughput, average delay and average power. 



Parameter 


Symbol 


Value 


communication threshold 


7c 


7 dB 


noise power spectral density 


N 


-90 dBm 


path loss exponent 


P 


4 


transmitter-receiver distance 


D 


100 m 


initial allocation interval of FCFS 


a 


2.6 s 


initial allocation interval of PCFCFS 


<Po 


2.54 s 


algorithm operation time 


T 


3 x 10 5 s 



Table 7.1: System parameters for performance evaluation of PCFCFS and FCFS algo- 
rithms. 

The system parameters for our numerical experiments are shown in Table 7.1. From 
(7.4) and (7.5), we obtain Pi = 0.50 mW and F 2 = 3.01 mW. We vary the arrival rate 
A from 0.40 to 0.60 packets/s in steps of 0.01. Figure 7.8 plots the throughput versus 
arrival rate for the PCFCFS and FCFS algorithms. Figure 7.9 plots the average delay 
per successful packet versus arrival rate for both the algorithms. Finally, Figure 7.10 
plots the average power per successful packet versus arrival rate for both the algorithms. 

For arrival rates exceeding 0.56, the throughput of PCFCFS is less than the arrival 
rate (Figure 7.8) and the average delay of PCFCFS increases rapidly (Figure 7.9), which 
leads to a substantial increase in the number of backlogged packets and system insta- 
bility. Hence, the maximum stable throughput of PCFCFS is between 0.55 and 0.56. 
Thus, Figures 7.8 and 7.9 corroborate our result that the maximum stable throughput 
of PCFCFS is 0.5518 (see Section 7.4). 
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Figure 7.8: Throughput versus arrival rate for PCFCFS and FCFS algorithms. 
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Figure 7.9: Average delay versus arrival rate for PCFCFS and FCFS algorithms. 
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Figure 7.10: Average power versus arrival rate for PCFCFS and FCFS algorithms. 

For both PCFCFS and FCFS, the departure rate (throughput) equals the arrival rate 
for all arrival rates up to 0.487 (Figure 7.8). Hence, both these algorithms are stable for 
arrival rates below 0.487. For arrival rates exceeding 0.487, the departure rate of FCFS 
is strictly lower than its arrival rate, leading to packet backlog and system instability. 
On the other hand, for PCFCFS, the departure rate still equals its arrival rate for arrival 
rates between 0.487 and 0.5518. In other words, the PCFCFS algorithm is stable for 
a higher range of arrival rates compared to FCFS algorithm. However, the PCFCFS 
algorithm becomes unstable for arrival rates exceeding 0.5518. 

The PCFCFS algorithm achieves higher throughput and lower average delay than 
the FCFS algorithm, albeit at the cost of expending higher average power. For example, 
at A = 0.55, PCFCFS achieves 13.3% higher throughput and 96.7% lower average delay 
than FCFS, at the cost of 170% higher power. 

7.6 Conclusions 

In this chapter, we have considered random access in wireless networks under the phys- 
ical interference model. By recognizing that the receiver can successfully decode the 
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strongest packet in presence of multiple transmissions, we have proposed PCFCFS, a 
splitting algorithm that modulates transmission powers of users based on observed chan- 
nel feedback. PCFCFS achieves higher throughput and substantially lower delay than 
those of the well known FCFS algorithm with uniform transmission power. We show 
that the maximum stable throughput of PCFCFS is 0.5518. PCFCFS can be imple- 
mented in those scenarios where users are willing to trade some power for a substantial 
gain in throughput. Moreover, if users can estimate the arrival rate of packets, then 
they can employ FCFS algorithm for arrival rates up to 0.4871 and PCFCFS algorithm 
for higher arrival rates, thus leading to further reduction in average transmission power. 



Chapter 8 

Flow Control: An Information 
Theory Viewpoint 



This thesis has so far explored various aspects of link scheduling in wireless networks. 
An equally interesting problem is to analyze flow control. We formulate the problem of 
controlling the rate of packets at the ingress of a packet network (possibly a wireless link) 
so as to maximize the mutual information between a source and a destination. We discuss 
various nuances of the problem and describe related work. We then derive the maximum 
entropy of a packet level flow that conforms to linearly bounded traffic constraints, by 
taking into account the covert information present in the randomness of packet lengths. 
Our results provide insights towards the design of flow control mechanisms employed by 
an Internet Service Provider (ISP). 

The rest of the chapter is organized as follows. In Section 8.1, we define the problem 
of information theoretic analysis of flow control in a packet network. In Section 8.2, we 
introduce a Generalized Token Bucket Regulator (GTBR) as our flow control mechanism. 
The concepts of flow entropy and information utility are defined in Section 8.3. We 
formulate the problem of determining the GTBR with maximum information utility in 
Section 8.4. In Section 8.5, we derive a necessary condition for the optimal GTBR and 
compute its parameters. We explain the results from an information theoretic viewpoint 
in Section 8.6 and discuss the implications of our work in Section 8.7. 
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8.1 System Model 
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Figure 8.1: Flow control of a source's packets over a packet network. 

Our system model is shown in Figure 8.1, wherein a source sends packets to a des- 
tination over a packet-switched network (possibly a wireless network). The packets 
transmitted by the source are regulated (policed) by a flow control mechanism at the 
ingress of the network. We are interested in the packet probability distribution that 
maximizes the mutual information between the source and the destination. In other 
words, given the description of the flow control mechanism and the stochastic charac- 
terization of the packet network, we seek the maximum amount of information (in the 
Shannon sense) that can be transmitted from the source to the destination. 

The problem can be stated as: 

max/(X;y), (8.1) 

Px 

where 

X = random variable representing randomness in packet contents, lengths 

and timings at the source, 
y = random variable representing randomness in packet contents, lengths 
and timings at the destination, 
px = probability distribution of X. 

(8.1) can be simplified to: 

max (H(X) - H(X\Y)) . (8.2) 

PX 

Thus, to maximize the information transfer from the source to the destination, we not 
only have to characterize the entropy of the source's packets H(X), but also the con- 
ditional entropy of the source's packets given the packets received at the destination 
H(X\Y). 

We state the following remarks about our problem formulation: 
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1. It is well-known that, in a packet-switched network, information can be transmit- 
ted not only by the contents, but also by the lengths and timings of packets. [20] 
is perhaps the first work to recognize this fact. Information transmitted by the 
lengths and timings of packets is referred to as covert information or side infor- 
mation. The channel that is used to convey covert information is called covert 
channel. Covert channels have been investigated in [20], [21], [90]. 

2. By flow control, we mean a rate control mechanism that regulates the packets 
transmitted by a source (subscriber) at the ingress of a network. Note that we 
do not consider end-to-end flow control mechanisms such as Transmission Con- 
trol Protocol (TCP). For simplicity, we consider a flow control mechanism that is 
described by a linearly bounded service curve 1 [91]. 

3. In the packet network, packets can be received incorrectly at the destination due 
to fluctuations in the channel, like that in a wireless channel. We assume the 
existence of link layer mechanisms such as Forward Error Correction (FEC) which 
ensure that all packets are correctly received at the destination. 

The packet network shown in Figure 8.1 only guarantees that the contents and lengths 
of the packets transmitted by the source are the same as those at the destination. 
However, the network can arbitrarily vary the timings between packets. Equivalently, 
the network can highly distort the covert timing information carried by the packets. 

Taking a cue from this, we only take into account information that is carried by 
the contents and lengths of the packets. Consequently, the probability distribution of 
packet contents and lengths at the destination is the same as that at the source. Hence, 
H(X\Y) = and (8.2) simplifies to 

maxH(X). (8.3) 

Px 

In other words, we seek the probability distribution of packet contents and lengths that 
maximize the source entropy H(X). 

1 Consider a flow through a system S with input and output functions A(t) and B(t) respectively. 
S offers to the flow a service curve $(t) if and only if is a wide sense increasing function, with 
i?(0) = 0, and B{t) ^ inf s ^ t {^(s) + d{t - s)} for all t ^ 0. 



150 



Chapter 8. Flow Control: An Information Theory Viewpoint 



Typically, the entity that owns a network, say an Internet Service Provider (ISP), 
implements certain mechanisms to ensure that packets transmitted by a subscriber are 
not lost in the network. However, to allocate network resources efficiently and guarantee 
zero loss of packets, the entity also mandates that the aggregate traffic of a subscriber 
be upper bounded by an envelope or a service curve. For example, the entity can 
mandate that the aggregate traffic of the subscriber be linearly bounded. A linearly 
bounded service curve can be implemented by a class of regulators known as token 
bucket regulators. 




Destination 



Source 



Figure 8.2: Token bucket regulation of a source's packets over a packet network. 

The system model that we analyze incorporates a Token Bucket Regulator (TBR) 
and is shown in Figure 8.2. A source transmits packets to a destination over a network, 
where every packet consists of an integer number of bits. The packets transmitted by the 
source are regulated by a TBR or leaky bucket regulator [92]. Intuitively, the regulator 
collects tokens in a bucket of depth B, which fills up at a certain rate r. Each token 
corresponds to the permission to transmit one bit into the network. The packets to be 
transmitted by the source accumulate in its data buffer over time. If there is a packet 
of length n bits in the data buffer at a given time, then it can be sent into the network 
only if n ^ B + r. If the packet is transmitted, then n tokens are depleted from the 
token bucket. 

A TBR can be used to smoothen the bursty nature of a subscriber's traffic. We as- 
sume that the network is owned by an ISP. From a Quality of Service (QoS) perspective, 
a TBR can be considered to be a part of the Service Level Agreement (SLA) between a 
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subscriber and an ISP. The SLA mandates that the ISP should provide end-to-end loss 
and delay guarantees to a subscriber's packets, provided the traffic profile of the sub- 
scriber adheres to certain TBR constraints. Specifically, the onus of the ISP is to ensure 
that every packet of a conforming source successfully reaches its destination within a 
certain permissible delay. 

The Standard Token Bucket Regulator (STBR), as defined by the Internet Engi- 
neering Task Force (IETF) and shown in Figure 8.2, enforces linear-boundedness on the 
flow. An STBR is characterized by its token increment rate r and bucket depth B. We 
will be more general and consider a TBR in which the token increment rate and bucket 
depth (maximum burst size) can vary from slot to slot. Such a TBR, which we define 
as a Generalized Token Bucket Regulator (GTBR), can be used to regulate Variable 
Bit Rate (VBR) traffic 2 from a source [93]. The continuous-time analogue of a GTBR 
is the time-varying leaky bucket shaper [94] in which the token rate and bucket depth 
parameters can change at specified time instants. 

The idea is to develop the notion of information utility of a GTBR. Specifically, 
we derive the maximum information that a GTBR-conforming traffic flow can convey 
in a finite time interval, by taking into account the additional information present in 
the randomness of packet lengths. These aspects are further elucidated in subsequent 
sections. 

8.2 Generalized Token Bucket Regulator 

In this section, we mathematically describe our system model and define a GTBR. We 
also explain the differences between our system model and those considered in existing 
literature. 

Consider a system in which time is divided into slots and a source which has to com- 
plete its data transmission within S slots. In our discrete-time model, we will evaluate 
the system at time instants 0, 1, . . . , S — l,S. Slot k is defined to be the time interval 
[k, k + 1), i.e., data transmission commences with slot and terminates with slot (S — l). 



2 For example, a pre-recorded video stream. 
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Figure 8.3: Relative time instants of parameters defined in (8.4). 
The traffic from the source is regulated by a GTBR. We define: 

rfc = token increment for slot k, 

B k = bucket depth for slot (k + l), 

£k = length of packet (in bits) transmitted in slot k, 

Uk = residual tokens at start of slot k. (8.4) 

rfc, B k , £k and u^, whose relative time instants are shown in Figure 8.3, are all non- 
negative integers. Let r := (r ,ri, . . . ,r S -i) denote the token increment sequence and 
B := (B , Bi, . . . , B s _ 2 ) denote the bucket depth sequence. The system starts with zero 
tokens. So, u — 0. A GTBR with the above parameters is denoted as 1Z g (S, r, B). 
The constraints imposed by 7Z g (S, r, B) on the packet lengths is 

ii^Ui + n V i = 0,l,...,S- 1. (8.5) 

If (8.5) is satisfied, then £ = (£o,£±, . . . ,£s-i) is a conforming packet length vector and 
the number of residual tokens will evolve according to 

u = 0, 

u i+ i = mm(ui + ri-£i,Bi) V i = Q,l,...,S -2, 
u s = us-i + rs-i- £s-i- (8.6) 

(8.6) is referred to as the token evolution equation. 

Note that if n = r V i = 0, 1, . . . , S - 1 and Bi = B V % = 0, 1, . . . , S - 2, then the 
GTBR K g (S, r, B) degenerates to the STBR K S (S, r, B). 

We should point out that our system model is similar to that of [95] . However, unlike 
[95], our traffic regulator is a deterministic mapping of an input sequence to an output 
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sequence. Also, the rate of our regulator is defined by the average token increment rate 
and not by the peak rate. 

The system model encompasses that of [90], wherein the authors have derived the 
information utility of an STBR and suggested a pricing viewpoint for its application. 
Our interest, however, is more theoretical. Specifically, we consider an STBR as a special 
case of a GTBR and describe a framework for their information-theoretic comparison. 
The main objective is to investigate whether a GTBR can achieve higher flow entropy 
than an STBR and explain the properties of entropy-maximizing GTBRs. These aspects 
are addressed in the following sections. 

8.3 Notion of Information Utility 

In this section, we introduce the concept of information utility of a GTBR. We de- 
rive the entropy of a flow that is regulated by a GTBR by considering the information 
present in the contents and lengths of the packets. We formulate the problem of com- 
puting the maximum flow entropy and subsequently describe a technique to compute 
the information utility of the GTBR. 

Consider a source which has a large amount of data to send and whose traffic is reg- 
ulated by a GTBR. We seek to maximize the information that the source can convey to 
the destination in the given time interval or the entropy present in the source traffic flow 
in an information-theoretic sense. For a given transmission interval S, token increment 
sequence r and bucket depth sequence B, the maximum entropy achievable by any flow 
which is constrained by the GTBR 7Z g (S, r,B) is termed as the information utility of 
the GTBR K g (-). 

The source can send information to its destination via two channels: 

i. Overt channel: The contents of each packet. Let £i denote the length of a packet in 
bits. The value of each bit is or 1 with equal probability and is independent of the 
values taken by the preceding and succeeding bits. Thus, this packet contributes 
li bits of information. 

ii. Covert channel: We consider the length of a packet as an event and associate a 
probability with it. Thus, side information is transmitted by the randomness in 
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the packet lengths. 

The joint entropy of i. and ii. is the sum of their entropies. 

During any slot k, the only method by which past transmissions can constrain the 
rest of the flow is by the residual number of tokens u k . So, Uk captures the state of 
the system. The key observation is that the future entropy depends only on the token 
bucket level u k in slot k. Hence, entropy is a function of system state u k and is denoted 
by H k (u k ). 

During slot S, the source signals the termination of the current flow by transmitting 
a special string of bits (flag). The information transmitted by this fixed sequence of bits 
is zero. Thus 

H s (u s ) = 0Vn s . (8.7) 

For a given state u k of the system, if a packet of length £ k bits is transmitted with 
probability p£ k (u k ), then 

1. The overt information transmitted is l k bits, 

2. As the event occurs with probability pe k (u k ), the covert information transmitted 
is (-log 2 p ek (u k )) bits, 

3. Since i k is random, u k+ \ is also random (from (8.6)). Thus, H k+ i{u k+ i) is also a 
random variable. 

Adding all of the above and averaging it over all conforming packet lengths, we obtain 
the entropy in the current slot (stage) 

H k (u k ) = Y Pt k (u k )\ik -log 2 (Pe k (u k )) + H k+1 (mm(u k + r k -£ k ,B k ))j 

4=0 

Vfc = 0,...,S'-l. (8.8) 

The equation above, which will be referred to as the flow entropy equation, intuitively 
states that the flow entropy of the current state is given by the sum of the entropy of the 
packet contents, the entropy of the packet lengths and the flow entropy of possible future 
states in the next slot. Note that (8.8) is similar to the backward recursion equation 
from dynamic programming [96]. Finally, the packet length probabilities must satisfy 

Pt k (u k ) = 1 Vfc = 0,...,5-1. (8.9) 
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Let p k (u k ) = (po(uk) , Pi(uk) , ■ ■ ■ ,p Uk+ r k (uk)) denote the vector of packet length proba- 
bilities for slot k with u k residual tokens. The dependence of pi h and on u k is assumed 
to be understood and is not always stated explicitly. So, p k = (po,pi, ■ ■ ■ ,p Uk +r k )- 

Our objective is to determine the sequence of probability mass functions (pg_i, Ps-2> 
" ' ' ) Po) which maximizes the flow entropy H o (0) for a given GTBR TZ g (S, r, B). From 
(8.7) 

H*(u s ) = 0. (8.10) 

From (8.8) 

H k {u k ) = ^ Pik (4 - ^g 2 {p tk ) + H k+i ( min(w fc + r k - 4, B k ))j 

4=0 

VJfe = 0,l,...S-l. (8.11) 

Given H£ +1 (u k+1 ) V Uk+i, there exists an optimum probability vector p£ — (pq,p\, ■ ■ ■ , Pu k+rk ) 
which maximizes the flow entropy H k (u k ), i.e., 

u k +r k 

H *k( u k) = Pi ~ lo ^(Pi k ) + K+1 ( min(w fc + r k - 4, B k ))j 

4=0 

V A; = 0,1,..., 5-1. (8.12) 

Thus, the problem of computing the entire sequence of probability vectors (Ps_!, 
P<y_2> ' ' ' j Po) decouples into a sequence of subproblems. The subproblem for slot 
k is: Given the function Hl +1 (u k +i) V u k +i, determine the probability vector p^ = 
(Po,Pi, • • • ,Pu k +r k ) so as to 

u k +r k 

maximize ^ p lk [£ k - log 2 (p 4 ) + H* k+1 ( min(w fe + r k - £ k , B k ))j , 

4=0 

u k +r k 

subject to ^ p£ k — 1. (8.13) 

4=0 

(8.13) is an equality-constrained optimization problem and can be solved using the 
technique of Lagrange multipliers [97]. Define the Lagrangian 

u k +r k 

C(p k , Afc) = ^ pe k [i k ~ log 2 (P4) + H k+i{ mm («fc + r k - 4, B k ))j 

4=0 

' u k +r k 



4=o J 
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At the optimal point (p* k , \* k ), we must have 

DC 



d Pik 
dC 



Solving (8.16) yields 



= V ^ 4 < u k + r k , (8.15) 



= 0. (8.16) 



Uk+T k 

J2pIM = !. ( 8 - 17 ) 

4=0 

which is (8.9) for the case of optimal probabilities. Solving (8.15), we obtain 

p* tk (u k ) = 2 4 ~ log2e+ ^+i (min(ufc+rfc - 4 ' Bfc))+A ^ Mfc) V < 4 < u k + r k . (8.18) 
From (8.17) and (8.18), the optimal Lagrange multiplier is given by 

2 tk+H *+^ n{ - Uk+rk ~ tk > Bk)) J . (8.19) 



From (8.18) and (8.19), the optimum packet length probability is given by 

24+fl'fc + i(min(u fc +r fc -4,B fc )) 

Pt k ( U k) = v « fc +r fe 9 a k +H*(mm(u k +r k -a k ,B k )) ' ( 8 - 20 ) 
2^a k =0 Z 

From (8.12) and (8.20), we finally obtain 

( Uk+T k \ 
2 tk+H ^ nin{ - Uk+rk -^ Bk)) \. (8.21) 
4=0 ' 

(8.21) will be referred to as the optimal flow entropy equation. 

The information utility of the GTBR 7Z g (S, r, B) is defined to be Hq(0), the maxi- 
mum flow entropy. Hq(0) is computed by starting with Hg(u s ) = 0, and using (8.21) 
to compute the optimal flow entropy Hl(u k ) for all u k and then proceeding backward 
recursively for k — S — 1, S — 2, . . . , 0. 



8.4 Problem Formulation 

Having developed a method to compute the information utility of a GTBR in Section 
8.3, we seek answers to the following questions: 
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a. Can a GTBR achieve higher information utility than that of an STBR? 

b. If yes, what is the increase in information utility? 

For the information-theoretic comparison of a GTBR TZ g (S, r,B) and an STBR 
1Z S (S', r, B), we impose the following conditions: 

1. TZg(-) and 1Z S (-) must operate over the same number of slots, i.e., 

S = S'. (8.22) 

2. The aggregate tokens of TZ g (-) and 1Z S (-) must be equal, i.e., 

5-1 

^Tr, = Sr. (8.23) 

3. The aggregate bucket depth of 1Z g (-) must not exceed that of 7Z S (-) 3 , i.e., 

5-2 

J2 B i < (S-1)B. (8.24) 

4. The bucket depth of TZ S (-) cannot be very high compared to its token increment 
rate. To quantify this, we mandate 4 

2r < B < 5r. (8.25) 



5. The token increment rate of 1Z g (-) in every slot must not exceed the bucket depth 
of TZ S (-), i.e., 

Ti ^ B. (8.26) 



3 Equality is present in (8.23) because every additional token directly translates to the permission 

to transmit one more bit, leading to increase in information utility. As this is not necessarily true for 

bucket depth, we permit inequality in (8.24). 

4 This assumption is practically justifiable. For example, in [94], the authors use r = 6 Mbps and 

B = 12 Mbps for their simulations. 
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If Conditions 1, 2, 3, 4 and 5 are satisfied, then GTBR 1Z g (-) and STBR 1Z S (-) are said 
to be comparable to each other. 

The optimal GTBR problem is formally stated as: 
Given an STBR 1Z S (S, r, B), determine the token increment sequence r and bucket depth 
sequence B of a comparable GTBR TZ g (S, r, B) so as to 

maximize Hq(0), 
subject to J2t =o r 'i = Sr > ( 8 - 27 ) 
Ef=o 2 ^ (8.28) 

Note that we are maximizing a real-valued function over two finite sequences of non- 
negative integers. 

8.5 Results 

In this section, we derive a necessary condition for the optimal GTBR in terms of 
aggregate bucket depth. We also compute the parameters of the optimal GTBR for 
some representative cases. 

8.5.1 Analytical Result 

Proposition 8.5.1. For an optimal GTBR, equality must hold in (8.28), except when 
S is small. In other words, ifB* is the bucket depth sequence of an optimal GTBR, it 
must satisfy 

S-2 

J^B* = (S-l)B. (8.29) 

i=0 

Proof. We prove by contradiction. Define gk{u) = 2 H k ( - u \ Since H%(u) ^ 0, guiu) ^ 1. 
From (8.21) 

u+r k 

g k {u) = ^2Vi(min(« + r,-f,B,)). (8.30) 



8.5. Results 



159 



gs-i(u) = 2" +rs - 1+1 — 1 is an increasing sequence in u. Using (8.30), we can show that 
g k {u) is an increasing sequence in u V k = 0, . . . , S — 1. Let /ij = maximum number of 
tokens possible in slot i. Thus 

Ho = 0, (8.31) 
fii = min(/ij_i + r i _i,S i _ 1 ) Vz = 1, . . . ,S - 1. (8.32) 

If Mj ^ /ij, then we say that state Ui is reachable in slot z, otherwise it is unreachable. 

Let 1Z g {S, r, B) be an optimal GTBR, for which equality does not hold in (8.28). 
Then Y. S i=S Bi^(S- l)B - 1. Consider another GTBR K' g (S, r', B') with r' = r and 
B' = (Bq, . . . , Bk-i, B k + 1, B k+1 , • • • , Bs-2) for some k. Let H' k *(u) denote the optimal 
flow entropy of H' g (') in slot k with u residual tokens. Define g' k (u) = 2 H 'k ("). From 
(8.21) 

u+r k 

g' k {u) = ^2^ +1 (min( W + ^-£,^)). (8.33) 

B' satisfies (8.28). g'^u) = gi(u) V i — k+1, . . . , S and V u. Since min(-u+rfc— £, B k + 1) ^ 
mm(u+r k —£ : B k ), it follows that g k (mm(u+r k —£, B k + l)) ^ g k (mm(u+r k —£ : B k )) ^ I. 
If we determine a reachable state u such that <7^(m) > gk( u ), then <7q(0) > fi'o(O), since 
the flow entropy in slot is computed slot-by-slot as a linear sum of future possible flow 
entropies with positive weights. Thus, the problem now reduces to determining a slot k 
and a reachable state u such that g' k (u) > gk(u). One of the following must hold: 

1. There exists ani e {1, . . . , S — 1} such that = B^i < + rj_i, or 

2. There is no z such that /ij = < + r^-i. 

Case 1: Consider the smallest z such that /i« = B^_i < + rj_i. Substituting 
A; = z — 1 in (8.30), we obtain 

U+Ti-l 

9i-i(u) = ^2 2^ i (miii(B i _i,« + r i _i - £)), 

u+ri_i— 1 "tt+r-i-i 

= ^ 2^(^-0 + ^ 2^ i (« + r i _i-^). (8.34) 

^=0 £=M+r 4 _i-Bi_i 
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Substituting k = i — 1 in (8.33), we obtain 

u+n-i 

9i-i( u ) = 2 £ ^(min(Si_i + l,M + ri_i -£)), 

u+rj_i— 1 u+r-i-i 

■■g'i-M = 2< ft (B i _ 1 + l)+ 2^ i (« + r i _ 1 -^).(8.35) 

(8.34) and (8.35) hold only if 

u + ri-i-Bi-i-l^O. (8.36) 

u = is a state which is reachable in the original system as well as in the primed 
system and satisfies (8.36). Since gi(u) is an increasing sequence in u, (8.34) and (8.35) 
imply 3 t '_i(/Ji-i) > gi-i(lJ,i-i). Consequently, ^(0) > g (0). 

Case 2: If no such i exists, then we must have Bi ^ r + h r; V z = 0, . . . , S — 2. 

Adding these (S — 1) inequalities and using n ^ B (from (8.26)), 

5-2 

^Bi ^ (5r - r s _i) + (Sr - r S -i - r 5 _ 2 ) + (Sr - r 5 _i - r 5 _ 2 - r 5 _ 3 ) H , 

i=0 

^ (Sr -B) + (Sr - 2B) + (Sr - 3B) + ■ ■ ■ , (8.37) 
= S(S-l)r-aB. (8.38) 

We cannot have r { = B V % (from (8.25), (8.26) and (8.27)). Thus, a cannot be of the 
order of S 2 . Thus, the lower bound on J2i=o Bi given by (8.37) and (8.38) is a loose 
lower bound. From (8.25), (8.28) and (8.38), E 4 = 5 i 

grows as S 2 and is upper-bounded 
by 5(S — l)r, which is impossible (except when S is small). So, we discard Case 2. 

From the result of Case 1, Hq*(0) > Hq(0). So, our assumption that 1Z g (-) is an 
optimal GTBR is incorrect. Therefore, equality must hold in (8.28) for every optimal 
GTBR. ■ 



8.5.2 Numerical Results 

For a given data transmission time S, token increment sequence r and bucket depth se- 
quence B, we determine the optimal GTBR by exhaustive search over the reduced search 
space obtained from Proposition 8.5.1. Our computation results are shown in Table 8.1. 
H s and H* denote the information utility of the STBR 1Z s (S,r,B) and the optimal 
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STBR 

parameters 
(S,r,B) 


optimal token increment 
sequence of GTBR 
r* 


optimal bucket depth 
sequence of GTBR 
B* 


information utility 


H s 
(bits) 


h; 

(bits) 


percentage 
increase 


(4,3,6) 


(6 3 3 0) 


(6 6 6) 


20.04 


20.92 


4.4% 


(4,3,7) 


(6 4 2 0) 


(6 8 7) 


20.08 


21.16 


5.4% 


(4,3,8) 


(7 3 2 0) 
(8 3 10) 


(7 9 8) 
(8 9 7) 


20.10 


21.32 


6.1% 


(4,3,9) 


(8 3 10) 
(9 2 10) 


(8 10 9) 
(9 10 8) 


20.10 


21.44 


6.7% 


(4,3,10) 


(9 3 0) 


(9 12 9) 


20.10 


21.51 


7.0% 


(4,3,11) 


(10 2 0) 
(11 1 0) 


(10 12 11) 
(11 12 10) 


20.10 


21.54 


7.2% 


(4,3,12) 


(12 0) 


(12 12 12) 


20.10 


21.56 


7.2% 


(4,3,13) 


(12 0) 


(13 13 13) 


20.10 


21.56 


7.2% 


(4,4,8) 


(8 4 4 0) 


(8 8 8) 


25.08 


26.04 


3.8% 


(4,4,9) 


(8 5 3 0) 
(9 4 3 0) 


(8 10 9) 
(9 10 8) 


25.11 


26.24 


4.5% 


(4,4,10) 


(9 5 2 0) 


(9 12 9) 


25.13 


26.39 


5.0% 


(4,4,12) 


(11 4 1 0) 


(11 14 11) 


25.14 


26.59 


5.8% 


(4,4,16) 


(16 0) 


(16 16 16) 


25.14 


26.70 


6.2% 


(4,5,10) 


(10 5 5 0) 


(10 10 10) 


29.91 


30.92 


3.4% 


(4,5,12) 


(11 6 3 0) 


(11 14 11) 


29.96 


31.24 


4.3% 


(4,6,12) 


(11 7 6 0) 
(12 7 5 0) 


(11 13 12) 
(12 13 11) 


34.60 


35.66 


3.1% 


(5,3,6) 


(6 3 3 3 0) 


(6 6 6 6) 


25.68 


26.57 


3.5% 


(5,3,9) 


(8 3 3 1 0) 


(8 10 10 8) 


25.88 


27.33 


5.6% 


(5,3,12) 


(11 2 2 0) 


(11 13 13 11) 


25.90 


27.59 


6.5% 


(5,3,15) 


(15 0) 


(15 15 15 15) 


25.90 


27.64 


6.7% 


(6,2,4) 


(4 2 2 2 2 0) 


(4 4 4 4 4) 


23.00 


23.77 


3.4% 


(6,3,6) 


(6 3 3 3 3 0) 


(6 6 6 6 6) 


31.33 


32.23 


2.9% 



Table 8.1: Entropy-maximizing GTBR for given data transmission time, token rate and 
bucket depth of a comparable STBR. 
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21.6 




bucket depth B of comparable STBR 



Figure 8.4: Information utility of GTBR vs. bucket depth of comparable STBR. 
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Figure 8.5: Information utility of GTBR vs. token increment rate of comparable STBR. 



8.5. Results 



163 



GTBR 1Z g (S, r*, B*) respectively. We also observe the variation in information utility 
of the optimal GTBR with important parameters of the comparable STBR, namely its 
bucket depth B and token increment rate r. For a data transmission time of 4 slots and 
token increment rate of 3 bits, Figure 8.4 shows the variation of information utility of 
the GTBR versus the bucket depth of the comparable STBR. For a data transmission 
time of 4 slots and bucket depth of 15 bits, Figure 8.5 shows the variation of information 
utility of the GTBR versus the token increment rate of the comparable STBR. 
Based on our computations, we draw the following inferences: 

1. A generalized token bucket regulator can achieve higher information utility than 
that of a standard token bucket regulator. The increase in information utility is 
significant (up to 7.2%), especially for higher values of B. 

2. The optimal bucket depth sequence B* is uniform 5 or near-uniform (the standard 
deviation is very small compared to the mean). 

3. The optimal token increment sequence r* is a decreasing sequence and is not 
uniform. 

4. For a fixed data transmission time S and token increment rate r: 

(a) If B = 2r, B* is always uniform and r* is uniform except for the terminal 
values. 

(b) As B increases from 2r to min(5, S)r, the variance of r* increases rapidly 
with a concentration of tokens in first few stages, the variance of B* increases 
slowly, while H* initially increases and then saturates at some final value. H* 
is an increasing and concave sequence 6 in B (see Figure 8.4). 

5. For a fixed data transmission time S and bucket depth B, H* an increasing, highly 
linear and slightly concave sequence in r (see Figure 8.5). 

±f — ±S 1 — ■ ■■ — D s _ 2 . 

6 The sequence of first-order differences {B\ — B$ , B 2 — B\, ■ ■ ■ , B* s _ 2 — B* s _ 3 ) is a decreasing and 
non-negative sequence. 
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8.6 Information-Theoretic Interpretation 

In this section, we provide explanations for empirical results in Section 8.5. The expla- 
nations are intuitive and rely on basic results from information theory. 

Consider a system with n states, where Pi denotes the probability of state i and 
^2™ =1 Pi = 1. From classical information theory, system entropy H increases with de- 
creasing Kullback-Leibler distance between the given probability mass function (pmf) 
and the uniform pmf [98]. H is maximized only if pi — ■ ■ ■ — p n — ^. Also, H* in- 
creases with n. Analogously, a GTBR can achieve higher information utility than that 
of an STBR because the pmfs of the packet lengths at each stage are more uniform and 
have a larger support. Recall that, for given r and B, information utility is computed 
recursively using (8.6) and (8.21). 

We argue that the optimal bucket depth sequence B* must be uniform or near- 
uniform for maximum information utility. If B* is neither uniform nor near uniform, 
then Bj = min^ Bi is much smaller than B. This restricts the range of values taken by 
Uj + i and £j + i (from (8.5) and (8.6)). The support of packet length pmfs at stage j + 1 is 
reduced, leading to lower flow entropy at stage j + 1 and consequently lower information 
utility. Thus, B* must be uniform or near-uniform to maximize the minimum support 
of packet length pmfs at each stage. In Table 8.1, the observation that min^ B * — B — 1 
or minj B* = B throughout corroborates our claim that B* is near-uniform. 

We argue that for maximum information utility, the optimal token increment se- 
quence r* must be a decreasing sequence, subject to ^ Bi for every i. If r; L > Bi for 
any i, then a packet of length zero cannot be transmitted in slot i (from (8.6)) and will 
have zero probability. This decreases the support of the packet length pmfs in slot i and 
leads to lower information utility. More importantly, from (8.8), 

n> 

H* (0) = J2pl(°){ i o-log 2 (pl(0))+H* 1 (min(r -e ,B ))y (8.39) 

^o=0 

The major contribution to information utility Hq(0) is from the support of the packet 
lengths [0,r ] and the pmf of the packet lengths (pq(0)), while the contribution from 
H*(-) is relatively smaller. So, to maximize Hq(0), r should be allowed to take its 
maximum possible value, subject to r ^ -Bo, and the pmf of the packet lengths should 
be close to the uniform pmf. The observation that r = B consistently in Table 8.1 
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corroborates this. Also, a high value of r leads to larger supports of packet length pmfs 
at intermediate and later stages. Similarly, the first few elements of r* tend to take large 
values till the aggregate tokens are exhausted. However, their contribution to Hq(0) is 
not as pronounced and equality may not hold in r« ^ B{. Thus, r* must be a decreasing 
sequence and the first few elements of r* tend to take their maximum possible values, 
subject to Ti ^ Bi, to achieve uniformity and larger supports of packet length pmfs at 
intermediate and later stages. 

This "greedy" nature of r* is evident when S and r are kept constant and B increases 
(Result 4b). A similar argument is applicable when S and B are kept constant and r 
increases (Result 5). The only difference is that a unit increase in r will necessarily 
increase H* by at least S bits (S bits are contributed by the packet contents alone, 
which also explains the dominant linear variation in Figure 8.5), while a unit increase in 
B will increase H* only by an amount equal to the difference in covert information. The 
increase in covert information is positive only if the optimal token increment and bucket 
depth sequences (r* and B*) result in larger support and more uniformity for the packet 
length pmfs. Indeed, when B increases beyond the maximum number of tokens possible 
at any stage (maxjjyUj}), clamping the residual number of tokens at every stage becomes 
ineffective and the system behaves as if bucket depth constraints were not imposed at 
all (Figure 8.4). 

8.7 Discussion 

In this chapter, we have studied linearly bounded flows over a packet network. We con- 
sidered a source whose traffic is regulated by a generalized token bucket regulator and 
which seeks to maximize the entropy of the resulting flow. Recognizing that the random- 
ness in packet lengths acts as a covert channel in the network, the source can achieve 
maximum entropy by sizing its packets appropriately. We have formulated the problem 
of computing the GTBR with maximum information utility in terms of constrained token 
increment and bucket depth sequences. A GTBR can achieve higher information utility 
than that of a standard IETF token bucket regulator. Finally, we have information- 
theoretically interpreted the observation that an entropy-maximizing GTBR always has 
a near-uniform bucket depth sequence and a decreasing token increment sequence. 
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Conclusions 



The recent revolution in wireless communications has motivated researchers and engi- 
neers alike to design ever better wireless networks that deliver high data rates to users. 
The joint design of physical and MAC layers is the key to breaking the "bandwidth 
bottleneck" of wireless networks, which has been the primary inspiration for this thesis. 

This thesis has focused on link scheduling in wireless mesh networks by taking into 
account physical layer characteristics. The assumption made throughout this thesis is 
that a packet is received successfully only if the SINR at the receiver exceeds a cer- 
tain threshold, termed as communication threshold. The thesis has also discussed the 
complementary problem of flow control. 

The first part of this thesis has considered link scheduling in STDMA wireless net- 
works. The network is modeled by a finite set of store-and-forward nodes that commu- 
nicate over a wireless channel characterized by propagation path loss. We have consid- 
ered two nuances of the scheduling problem: point to point link scheduling wherein a 
transmitted packet is intended for a single neighbor only, and point to multipoint link 
scheduling wherein a transmitted packet is intended for all neighbors in the vicinity. 

Specifically, in Chapter 2, we have introduced the system model of an STDMA 
wireless network. We have discussed two prevalent models for specifying the criteria for 
successful packet reception: the protocol interference model which mandates a "silence 
zone" around a receiver and is better suited to represent WLANs, and the physical 
interference model which mandates that the SINR at a receiver be no less than the 
communication threshold and is more appropriate to represent mesh networks. We have 
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described the equivalence between a link schedule and the coloring of edges of a certain 
graph representation of the network, termed as communication graph. We have argued 
that STDMA link scheduling algorithms can be broadly categorized into three classes: 
those based entirely on a communication graph representation of the network, those 
based on communication graph and SINR threshold conditions and those based entirely 
on an SINR graph representation of the network. We have reviewed representative 
research papers from each of these classes. We have described limitations of algorithms 
that are based only on the communication graph. Subsequently, we have introduced 
spatial reuse as a performance metric that corresponds to aggregate network throughput. 

Next, in Chapter 3, we have critically examined ArboricalLinkSchedule, a point to 
point link scheduling algorithm proposed in [16] . While this is one of the earlier works on 
link scheduling with nice theoretical properties, it could yield a higher schedule length 
in practice. Specifically, the methodology employed by ArboricalLinkSchedule is to 
represent the network by a communication graph, partition the graph into minimum 
number of subgraphs and color each subgraph in a greedy manner. We have modified 
the algorithm to reuse colors while coloring successive subgraphs of the communication 
graph. We have shown that the modified algorithm yields lower schedule length in 
practice, albeit at a cost of slightly higher running time complexity. Subsequently, 
we have proposed the ConflictFreeLinkSchedule algorithm that not only utilizes the 
communication graph, but also verifies SINR threshold conditions at receivers. We have 
demonstrated that the proposed algorithm achieves higher spatial reuse than existing 
algorithms, even under fading and shadowing channel conditions. We have argued that 
the running time complexity of the proposed algorithm is only marginally higher than 
those of existing algorithms. 

Taking a step ahead, in Chapter 4, we have provided a somewhat different per- 
spective on point to point link scheduling. For an STDMA network, we recognize that 
interferences between pairs of links can be embedded into edge weights and normalized 
noise powers at receivers of links can be embedded into vertex weights of a certain graph 
representation of the network, termed as SINR graph. We have then proposed SINR- 
GraphLinkSchedule, a novel link scheduling algorithm that is based on the SINR graph. 
We have proved the correctness of the algorithm and shown that it has polynomial run- 
ning time complexity. We have demonstrated that the proposed algorithm achieves high 



169 



spatial reuse compared to algorithms which utilize a communication graph model of the 
network, including ConflictFreeLinkSchedule algorithm. 

In Chapter 5, we have considered point to multipoint link scheduling and generalized 
the definition of spatial reuse for this scenario. We have proposed a scheduling algorithm 
based on a communication graph representation of the network and "neighbor-average" 
SINR threshold conditions. Moreover, we have demonstrated that the proposed algo- 
rithm achieves higher spatial reuse than existing algorithms, without any increase in 
running time complexity. 

Overall, we have observed the tradeoff between accuracy of the network representa- 
tion, spatial reuse and algorithm running time complexity in our successive results. For 
a more accurate network representation, higher spatial reuse is achieved, but at a cost of 
higher running time complexity. For example, since the SINR graph representation of an 
STDMA network is more accurate than the communication graph representation, SINR- 
GraphLinkSchedule achieves higher spatial reuse than that of ConflictFreeLinkSchedule, 
but at a cost of increased running time complexity. 

A summary of existing and proposed link scheduling algorithms investigated in the 
first part of the thesis is provided in Table 9.1. 



Type of link 
scheduling 


Wireless network 
model 


Existing 
algorithms 


Proposed 
algorithm 


Point 

to 
point 


communication graph 


ArboricalLinkSchedule [16] 


ALSReuseColors 
(Chapter 3) 


communication graph 
and 

SINR conditions 


GreedyPhysical [27] 
TGSA [32] 


ConflictFreeLinkSchedule 
(Chapter 3) 


SINR graph 




SINRGraphLinkSchedule 
(Chapter 4) 


Point 
to 

multipoint 


communication graph 


BroadcastSchedule [16] 




communication graph 
and 

SINR conditions 




MaxAverageSINRSchedule 
(Chapter 5) 



Table 9.1: Link scheduling algorithms investigated in Chapters 3, 4 and 5. 



The second part of this thesis has considered link scheduling in random access wireless 
networks. Specifically, it has focused on random access algorithms for wireless networks 
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that take into account channel effects and SINR conditions at the receiver. 

In Chapter 6, we have reviewed representative research papers on such random access 
techniques. We have also motivated the use of variable transmission power in random 
access wireless networks. 

Subsequently, in Chapter 7, we have investigated a random access scenario wherein 
multiple transmitters (users) attempt to communicate with a single receiver over a wire- 
less channel characterized by propagation path loss. We have assumed that the receiver 
is capable of power based capture and proposed an interval splitting algorithm that 
varies transmission powers of users based on their arrival times and quaternary chan- 
nel feedback. We have modeled the algorithm dynamics by a Discrete Time Markov 
Chain and consequently shown that its maximum stable throughput is 0.5518. We have 
demonstrated that the proposed algorithm has higher throughput and lower delay than 
the FCFS interval splitting algorithm with uniform transmission power. 

The third and final part of this thesis has considered information-theoretic analysis 
of flow control in packet networks. We have defined the problem of maximizing the 
information carried by packets from a source to a destination, subject to a flow control 
mechanism at the ingress of the network. We have considered a linearly bounded flow and 
focused on the information carried by the randomness in packet contents and lengths. 
Consequently, we have formulated the problem of maximizing the entropy of a packet 
level flow that is shaped by a generalized token bucket regulator. We have demonstrated 
that the optimal regulator has a decreasing token increment sequence and a near-uniform 
bucket depth sequence. Finally, we have provided information theoretic interpretations 
for these observations. 

To sum it up, in this thesis, we have investigated both fixed and random access 
flavors of link scheduling problems in wireless networks from a physical layer viewpoint. 
Finally, we have discussed a flow control problem in packet networks. 

Various avenues for further research have emerged from our investigations. We out- 
line some possible directions for future work. 

I. It would be interesting to derive approximation bounds of ConflictFreeLinkSched- 
ule and SINRGraphLinkSchedule algorithms under reasonable assumptions on 
node deployment and interference regions. The assumptions and approximation 
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techniques employed in [27] may provide some pointers in this direction. 

2. Though distributed link scheduling algorithms for STDMA wireless networks un- 
der the protocol interference model have been proposed in [31], [99], the design 
of distributed link scheduling algorithms under the physical interference model 
remains a challenging problem. 

3. Various generalizations of the PCFCFS algorithm are worth investigating. For 
example: 

(a) Design a variable power splitting algorithm under the assumption that users 
are at unequal distances from the receiver and can adjust their minimum 
transmission powers accordingly. 

(b) Design a splitting algorithm for the case when the receiver is capable of de- 
coding more than one packet correctly (as in wideband systems) and the users 
can employ n transmission power levels, where n > 2. 

(c) Analyze the throughput improvement in CSMA/CA based WLANs when 
power control is employed in conjunction with binary exponential backoff. 
The work done in [100] can be a useful starting point. 

4. A challenging task would be to analyze the expected delay of the PCFCFS algo- 
rithm. A useful starting would be [101], [102], which have employed techniques to 
obtain upper and lowers bounds on the expected delay of the FCFS algorithm. 

5. Our results in Chapter 8 show the existence of upper bounds on the entropy of 
regulated flows. It would be interesting to construct source codes which come close 
to this bound. Furthermore, it would be insightful to develop a rate-distortion 
framework for a generalized token bucket regulator, perhaps using the techniques 
employed in [95]. 



Appendix A 



Proofs of Limiting Transition 
Probabilities 



According to L'Hopital's Rule, if \im x ^ c f (x) and \im x ^ c g(x) are both zero or are both 
±00 and, if linx^c has a finite value or if the limit is ±00, then 



hm — - = hm— — . (A.l) 

x^c g[x) x^c g\x) 



We will employ L'Hopital's Rule to prove (7.54) - (7.66) 

In this appendix, we will only provide the proofs of (7.54), (7.55), (7.56) and (7.57). 
The proofs of (7.58) - (7.66) are similar to those of (7.54) - (7.57) and are omitted. 



A.l Proof of (7.54) 

Proof. In (7.25), substitute G, t = x. From (7.7), = 2G, t = 2x. As i -> 00, 

Gi = 2~' l \(J) — > 0. Thus, using L'Hopital's Rule successively, we obtain 
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lim P L > w. 



lim P L , m 



lim 

x^O 



(1 — e x )xe 



— e )xe 
(l + 2x)e~ 2x ' 



lim 



lim 

x^O 

lim 

x^O 



xe 2x ) 



dx 



■(1 - e- 2x - 2xe~ 2x Y 
e~ x + xe~ x — e" 2x 

Axe~ 2x 
^(e~ x + xe~ x — e~ 2x ) 



lim 

x^O 
1 

2' 



2e 
4e" 23; 



dx 
-2x 



l -(Axe~ 2x ) 



- xe x 



■ 8xe 2x ' 



A. 2 Proof of (7.55) 

Proof. In (7.26), substitute Gi = x. Thus, using = 2x and applying L'Hopital's 
Rule successively, we obtain 

lim Pr / L i = lim ; — - — , 

i^oo L i> L i+i x ^ 1 - (1 + 2x)e~ 2x 

4-{e~ x — e~ 2x — xe~ 2x ) 
= lim# -, 

x^o A(i _ e -2x _ 2xe- 2x ) 

e -2x _|_ 2xe~ 2x — e~ x 
= lim , 

x->o Axe 

-f(e~ 2x + 2xe~ 2x - e~ x ) 
= lim^ -= -, 

£(Axe~ 2x ) 

e -x _ 4 X g-2x 

= lim 7: 7T- , 

x-+o Ae~ 2x - 8xe~ 2x 

.'. lim P L i ti — —. 



A.3 Proof of (7.56) 

Proof. In (7.27), substitute Gi = x. Thus, using = 2x and applying L'Hopital's 
Rule successively, we obtain 



A.4. Proof of (7.57) 
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lim P L 'c i+1 = I™ 



.-. lim P L , C 



+1 



x^o 1 - (1 + 2x)e~ 2xl 

l im dxU^ ) 

x^o £(1 - e ~ 2x - 2xe~ 2x ) 

lim 

4e~ 2x 

1 

8' 



5 



A.4 Proof of (7.57) 



Proof. In (7.28), substitute Gi = x. Thus, using = 2x and applying L'Hopital's 
Rule successively, we obtain 



1 - (l + x + ^)e~ x 
lim Pt/ r,.,, = lim ; r — - — , 

woo x^o 1 - (1 + 2x)e~ 2x ' 

lim — 6 ^ — xe ~ :r ~ T e ~ a ) 
x->o A(i _ e -2^ _ 2xe - 2;c ) ' 

1 _ X 

= lim , 

x^o Ae~ x 

lim P L 'l,. 4 
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